Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 715,749 T→C 8.0% V113A (GTT→GCT)  ybfP → lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913715,7490TC 8.0% 210.0 / 7.6 75V113A (GTT→GCT) ybfPlipoprotein
Reads supporting (aligned to +/- strand):  ref base T (39/30);  new base C (3/3);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

CGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATA  >  NC_000913/715601‑715894
                                                                                                                                                    |                                                                                                                                                 
cGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGt                                                                                                                                                   >  2:603119/1‑149 (MQ=255)
  tggtgATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATTCGGGttt                                                                                                                                                 >  2:522945/1‑149 (MQ=255)
     tgATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGcc                                                                                                                                              <  2:777120/149‑1 (MQ=255)
      gATAACTTACCTGATAATACCCGCCGCTTGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCa                                                                                                                                             <  2:47825/149‑1 (MQ=255)
      gATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCa                                                                                                                                             <  1:986005/149‑1 (MQ=255)
      gATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCa                                                                                                                                             >  1:47825/1‑149 (MQ=255)
         aaCTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAggg                                                                                                                                          <  2:263500/149‑1 (MQ=255)
           cTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCa                                                                                                                                        >  1:110142/1‑149 (MQ=255)
              aCCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGTGCAAGc                                                                                                                                     >  2:1124722/1‑149 (MQ=255)
                cTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGccc                                                                                                                                   >  1:327555/1‑149 (MQ=255)
                 tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAACAGAACTTACTTATCTTCATATTATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCg                                                                                                                                  <  1:63150/149‑1 (MQ=255)
                 tGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCg                                                                                                                                  >  2:733781/1‑149 (MQ=255)
                   ataataCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGCTTGCCCTGGCAAACCCgca                                                                                                                                <  2:323941/149‑3 (MQ=255)
                        aCCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCg                                                                                                                           >  1:689931/1‑149 (MQ=255)
                         cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGa                                                                                                                          <  2:1000511/149‑1 (MQ=255)
                         cccGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGa                                                                                                                          >  1:1131040/1‑149 (MQ=255)
                          ccgccgCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCCCAAGAAAAGGTCGGTTTAAATCCCGGATGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCTGCTCTCGAc                                                                                                                         >  1:767041/1‑149 (MQ=255)
                          ccgccgCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGAc                                                                                                                         >  2:4678/1‑149 (MQ=255)
                                   gCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCa                                                                                                                >  1:322099/1‑149 (MQ=255)
                                     aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAATAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTGTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAt                                                                                                              <  2:737019/149‑1 (MQ=255)
                                     aGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTAAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGCGTTTGCCAGGGCAAGCCCGGCAATCGACCGAAAGTCaaa                                                                                                              >  1:1057916/1‑148 (MQ=255)
                                      gTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGGTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCCAGCACGGCACTCGACCGATAGTCAAtt                                                                                                             >  1:481513/1‑148 (MQ=255)
                                         gTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAgg                                                                                                          >  2:137625/1‑149 (MQ=255)
                                            aGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGAtt                                                                                                       >  2:151248/1‑149 (MQ=255)
                                              gCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTAt                                                                                                     >  2:45950/1‑149 (MQ=255)
                                               cGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATg                                                                                                    <  1:527168/149‑1 (MQ=255)
                                                      gAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTAAACCTCATTTTTTATGCGGGTTTGCCATTGCAAGCACGGAACTCGACCGATATTCAATGAGGATTATGGCAGCGa                                                                                             >  2:840288/1‑149 (MQ=255)
                                                        aaatGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGaaa                                                                                           <  1:817085/149‑1 (MQ=255)
                                                        aaaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGaaa                                                                                           <  1:1055267/149‑1 (MQ=255)
                                                         aaaGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAg                                                                                          >  2:243046/1‑149 (MQ=255)
                                                          aaGAACTTACTTATCATCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGc                                                                                         <  1:21475/149‑1 (MQ=255)
                                                                    ttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGTTTTGCCAGGGCAAGCCCGGCCCTCGACCGATAGTCAATGAGGATTATGGAAGCGAAAGCGGAATTTTTg                                                                               >  2:58676/1‑149 (MQ=255)
                                                                    ttaTCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCAGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTg                                                                               >  2:739823/1‑149 (MQ=255)
                                                                          tCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                         <  1:1170364/149‑1 (MQ=255)
                                                                          tCATAATATGTCAGCACAAGAAAAGGGCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCga                                                                         <  2:462297/149‑1 (MQ=255)
                                                                            ataataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTAAGCGAAAGCGGATTTTTTGGCTCgagt                                                                       >  1:573629/1‑147 (MQ=255)
                                                                             taataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatga                                                                      >  1:200721/1‑149 (MQ=255)
                                                                              aataTGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCgatgat                                                                     <  1:733781/149‑1 (MQ=255)
                                                                               atatGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTTGCTCGatgatg                                                                    >  2:1098103/1‑149 (MQ=255)
                                                                                  tGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCAt                                                                 >  1:411821/1‑149 (MQ=255)
                                                                                    tCAGCTCAAGAAAAGGTCGGTTTAAATCCAGGTTGGACATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGTCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGt                                                               <  1:207893/149‑1 (MQ=255)
                                                                                       gCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCgg                                                            >  2:1015340/1‑149 (MQ=255)
                                                                                          caAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGat                                                         <  1:341577/149‑1 (MQ=255)
                                                                                          caAGAAAACGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTATGCCAGCGCAAGCCAGGCACTCGACCGAAAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGAAGATGCATGTCGGGat                                                         >  1:352477/1‑149 (MQ=255)
                                                                                             gAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGattat                                                      >  1:1122896/1‑149 (MQ=255)
                                                                                             gAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGAAGCATGTCGGGattat                                                      >  1:348554/1‑149 (MQ=255)
                                                                                               aaaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTg                                                    <  1:603119/149‑1 (MQ=255)
                                                                                                aaGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGa                                                   >  1:532620/1‑149 (MQ=255)
                                                                                                           aaaTCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCt                                        >  2:526975/1‑149 (MQ=255)
                                                                                                               cccGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCa                                    >  2:77356/1‑149 (MQ=255)
                                                                                                                 cGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                  <  1:1124722/149‑1 (MQ=255)
                                                                                                                 cGGTTGGCAATGCTACACCTCATTTTATATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGt                                  <  2:406789/149‑1 (MQ=255)
                                                                                                                   gTTGGCAATGCTACACCTCATTTTTTATTCGGGTATGCCAGGGAAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTTTTGAGTTACGACGCTGCCAGTCa                                <  1:694336/149‑1 (MQ=255)
                                                                                                                    ttGGCAATGCTACACCTCATTTTTTATGCGGGTTTGACATGGCAAGCCCGGGACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAg                               <  1:155191/149‑1 (MQ=255)
                                                                                                                         aaTGCTACACCTCATTTTTTATGCGGGTTTGACATTGCAAGCCCGGCACTCGACAGATAGACAATGATGATTATGTAAGCTAAAGCGGATTTTTTGGATCGATGATGCATGTCGGGAATATTGAGTTACGACGCTGCCAGTCAGAAAAt                          >  2:697981/1‑149 (MQ=255)
                                                                                                                          aTGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAAtt                         <  1:85508/149‑1 (MQ=255)
                                                                                                                            gCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAACGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGc                       >  2:1058852/1‑149 (MQ=255)
                                                                                                                              tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                     >  2:302053/1‑149 (MQ=255)
                                                                                                                              tACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCa                     <  2:1131040/149‑1 (MQ=255)
                                                                                                                                  ccTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCaa                 <  1:137625/149‑1 (MQ=255)
                                                                                                                                    tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGCCCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa               >  1:1008889/1‑149 (MQ=255)
                                                                                                                                    tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa               <  2:689931/149‑1 (MQ=255)
                                                                                                                                    tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa               >  1:491935/1‑149 (MQ=255)
                                                                                                                                    tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGa               >  2:980782/1‑149 (MQ=255)
                                                                                                                                    tCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAAGGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGGCGGGAGTTTTGGGTTACGCCGCGGCCAGTCAGCAAATTGccggaatca               >  2:217884/1‑142 (MQ=255)
                                                                                                                                     cATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCTATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa              <  2:926595/149‑1 (MQ=255)
                                                                                                                                     cATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGaa              <  1:77356/149‑1 (MQ=255)
                                                                                                                                        tttttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTg           <  1:45950/149‑1 (MQ=255)
                                                                                                                                          tttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaa         <  2:322099/149‑1 (MQ=255)
                                                                                                                                           ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa        <  2:532620/149‑1 (MQ=255)
                                                                                                                                           ttATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGaaa        >  1:1005352/1‑149 (MQ=255)
                                                                                                                                             tgccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGCATATTTTGGCACGATGATGCATGTCGGGATTATTGAGTTACGACGCTGACAGTCAGAAAAATGCCAGCAAGAATTGAAAGc      >  2:490918/4‑149 (MQ=255)
                                                                                                                                              tGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCg     <  1:753870/149‑1 (MQ=255)
                                                                                                                                               ccGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGa    >  1:481823/2‑149 (MQ=255)
                                                                                                                                                cGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt            >  1:940374/1‑140 (MQ=255)
                                                                                                                                                cGGGCTTGCCCTGGCAAACCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAAtt            <  2:940374/140‑1 (MQ=255)
                                                                                                                                                 gggTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATa  >  1:271669/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
CGTGGTGATAACTTACCTGATAATACCCGCCGCTGGCAGTTGTCAGGCGGGAACGAAAAAGAACTTACTTATCTTCATAATATGTCAGCACAAGAAAAGGTCGGTTTAAATCCCGGTTGGCAATGCTACACCTCATTTTTTATGCGGGTTTGCCAGGGCAAGCCCGGCACTCGACCGATAGTCAATGAGGATTATGTCAGCGAAAGCGGATTTTTTGGCTCGATGATGCATGTCGGGATTATTGAGTTACGACGCTGCCAGTCAGAAAATTGCCAGCAAGAATTGAAAGCGATA  >  NC_000913/715601‑715894

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: