Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,421,938 T→C 6.6% W76R (TGG→CGG)  moeA → molybdenum::molybdopterin ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,421,9380TC 6.6% 176.0 / 2.3 60W76R (TGG→CGG) moeAmolybdenum::molybdopterin ligase
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base C (2/2);  total (27/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 = 8.25e-01

GGCCGAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAG  >  NC_002947/2421790‑2422085
                                                                                                                                                     |                                                                                                                                                   
ggcCGAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACt                                                                                                                                                     >  8:17184/1‑149 (MQ=255)
 gcCGAAGCTACACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTg                                                                                                                                                    >  5:290326/1‑149 (MQ=255)
   cGAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTgg                                                                                                                                                   >  1:225516/1‑148 (MQ=255)
   cGAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGc                                                                                                                                                  <  2:104/149‑1 (MQ=255)
    gAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTgg                                                                                                                                                   >  3:170107/1‑147 (MQ=255)
     aaGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAg                                                                                                                                                <  4:170107/149‑1 (MQ=255)
        cTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGc                                                                                                                                             <  5:54029/149‑1 (MQ=255)
           cACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAg                                                                                                                                          >  6:333535/1‑149 (MQ=255)
            aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                         <  8:89155/149‑1 (MQ=255)
            aCCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGc                                                                                                                                         >  4:39383/1‑149 (MQ=255)
             ccGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGcc                                                                                                                                        >  6:126558/1‑149 (MQ=255)
                   cgcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                   <  8:157596/148‑1 (MQ=255)
                    gcgAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCg                                                                                                                                 <  4:166988/149‑1 (MQ=255)
                      gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                               <  6:290326/149‑1 (MQ=255)
                      gagaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                               <  3:39383/149‑1 (MQ=255)
                        taCCGCAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACCGCGCCATGGACGGGTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                             <  8:269481/148‑1 (MQ=255)
                        gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                             <  7:91924/149‑1 (MQ=255)
                        gaCCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTca                                                                                                                             <  5:243257/149‑1 (MQ=255)
                             aaaaCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCa                                                                                                                                           <  8:239567/130‑1 (MQ=255)
                             aaaaCCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCa                                                                                                                                           >  7:239567/1‑130 (MQ=255)
                                   tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                               >  6:5052/1‑136 (MQ=255)
                                   tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGt                                                                                                                               <  5:5052/136‑1 (MQ=255)
                                   tGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtt                                                                                                                  >  3:110140/1‑149 (MQ=255)
                                     ccGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTAtttt                                                                                                                >  3:64865/1‑149 (MQ=255)
                                           gctgAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCgg                                                                                                          <  5:252924/146‑1 (MQ=255)
                                                 gCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATgc                                                                                                    >  2:155077/1‑149 (MQ=255)
                                                   cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                              <  5:287703/121‑1 (MQ=38)
                                                   cGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTc                                                                                                                              >  6:287703/1‑121 (MQ=38)
                                                       ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGTGCCCCa                                                                                              <  2:110557/149‑1 (MQ=255)
                                                       ggCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCa                                                                                              <  4:64865/149‑1 (MQ=255)
                                                         ccGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCTGGGCGGAACACCCCAATGGACGGTTATGCCTTGCGCATGGCTGCCTGGCAGGCCCAGCCGCTGCCGGTCAGCCAGCGTATTTACGCCGGCCATGCGCCCCAg                                                                                             >  1:79733/1‑149 (MQ=255)
                                                             gTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCt                                                                                        <  1:155077/149‑1 (MQ=255)
                                                                 tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCaa                                                                                    >  3:72314/1‑149 (MQ=255)
                                                                 tGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCaa                                                                                    >  7:46488/1‑149 (MQ=255)
                                                                   gCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGTCGGCCATGCGCCCCAGCCGCTGCAAc                                                                                   <  4:72314/148‑1 (MQ=255)
                                                                                    ctcgctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGa                                                                 >  2:286994/1‑149 (MQ=255)
                                                                                     tcgctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCa                                                                                                                                                 >  1:8303/1‑68 (MQ=255)
                                                                                     tcgctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCa                                                                                                                                                 <  2:8303/68‑1 (MQ=255)
                                                                                      cgctcgACCTGCCGCCCT‑GGCCGAACAGCACCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATc                                                               <  2:79733/149‑1 (MQ=255)
                                                                                        ctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                   >  7:147867/1‑79 (MQ=255)
                                                                                        ctcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGc                                                                                                                                   <  8:147867/79‑1 (MQ=255)
                                                                                        ctccACCTCCGGCCCG‑GGCCCCACCGGGCCTCGGCCGGTGAGGCCTTGCGCCTGGCGGCCTGGCAGGGCCAGCCGCTGCCGCGCAGCCAGCGGATTTTCCCCGGCCATGCGCCCGAGCCGCGGCATCCGGGTCCCTGCGCCCTGATCtt                                                             <  4:266797/149‑1 (MQ=255)
                                                                                         tcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                             <  7:116815/84‑1 (MQ=21)
                                                                                         tcgACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACCGGCAGCGGCTGGCCCTGCCAGTca                                                                                                                             >  8:116815/1‑84 (MQ=21)
                                                                                           gACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTc                                                            <  7:17184/147‑1 (MQ=255)
                                                                                           gACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCAc                                                          <  5:26370/149‑1 (MQ=255)
                                                                                                    cccT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGcc                                                 >  3:248065/1‑149 (MQ=255)
                                                                                                           ccGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGcc                                           >  4:36050/1‑149 (MQ=255)
                                                                                                            cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                          >  2:232683/1‑149 (MQ=255)
                                                                                                            cGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCACTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCg                                          <  5:333535/149‑1 (MQ=255)
                                                                                                                    gcCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGc                                  <  4:248065/149‑1 (MQ=255)
                                                                                                                             ggTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGt                         <  1:286994/149‑1 (MQ=255)
                                                                                                                             ggTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGt                         <  2:225516/149‑1 (MQ=255)
                                                                                                                              gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                        <  1:78478/149‑1 (MQ=255)
                                                                                                                              gTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTc                        <  5:126558/149‑1 (MQ=255)
                                                                                                                                        gcgcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCa                <  4:110140/147‑1 (MQ=255)
                                                                                                                                          gcCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGa             >  5:110252/1‑148 (MQ=255)
                                                                                                                                               gCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAacac       >  7:3558/1‑149 (MQ=255)
                                                                                                                                                 tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc      >  3:212625/1‑148 (MQ=255)
                                                                                                                                                 tGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACAcc      >  2:77328/1‑148 (MQ=255)
                                                                                                                                                   aCTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCg     >  4:158398/1‑147 (MQ=255)
                                                                                                                                                     tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTTTTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg  >  3:208538/1‑148 (MQ=255)
                                                                                                                                                     tGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAg  >  6:246213/1‑148 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCCGAAGCTGCACCGATCCGCGAGACCGAAAACCTGCCGCTGGCCGAAGCCGAGGGCCGCGTGCTGGCCACCGACCTGGTGGCCTCGCTCGACCTGCCGCCCT‑GGCCGAACAGCGCCATGGACGGTTATGCCTTGCGCCTGGCTGACTGGCAGGGCCAGCCGCTGCCGGTCAGCCAGCGTATTTTCGCCGGCCATGCGCCCCAGCCGCTGCAACCGGGTACCTGCGCGCGGATCTTCACCGGCGCGCCTGTGCCGGAGGGGGCCGACTGCGTCGAGATGCAGGAAAACACCGAAG  >  NC_002947/2421790‑2422085

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