Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,644 G→A 5.0% T327M (ACG→ATG)  cobN ← aerobic cobaltochelatase subunit CobN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,6440GA 5.0% 1222.6 / 31.6 377T327M (ACG→ATG) cobNaerobic cobaltochelatase subunit CobN
Reads supporting (aligned to +/- strand):  ref base G (169/189);  new base A (6/13);  total (175/202)
Fisher's exact test for biased strand distribution p-value = 2.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 383 total aligned reads displayed.

TCTGGGCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAG  >  NC_002947/3984495‑3984792
                                                                                                                                                     |                                                                                                                                                    
tCTGGGCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGcgcg                                                                                                                                                      <  1:8156105/150‑1 (MQ=255)
   gggCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGACGCGTGa                                                                                                                                                   <  2:7778018/150‑1 (MQ=255)
    ggCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGAt                                                                                                                                                  >  1:1907020/1‑150 (MQ=255)
    ggCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGAt                                                                                                                                                  >  2:3915151/1‑150 (MQ=255)
     gCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATa                                                                                                                                                 >  2:5429391/1‑150 (MQ=255)
         tGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgc                                                                                                                                             <  2:7290208/150‑1 (MQ=255)
         tGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgc                                                                                                                                             <  1:3785647/150‑1 (MQ=255)
         tGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgc                                                                                                                                             <  2:7231750/150‑1 (MQ=255)
           ggCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgcgc                                                                                                                                           >  2:8545600/1‑150 (MQ=255)
           ggCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATgcgc                                                                                                                                           <  2:432897/150‑1 (MQ=255)
              aGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCg                                                                                                                                        >  1:6347464/1‑150 (MQ=255)
              aGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCg                                                                                                                                        >  2:2813633/1‑150 (MQ=255)
               ggCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGt                                                                                                                                       <  1:1510176/150‑1 (MQ=255)
                 cGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTcc                                                                                                                                     <  2:2046469/150‑1 (MQ=255)
                   ggCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCATTATCACGCGCCCggtcag                                                                                                                                   <  2:5591184/150‑6 (MQ=255)
                    gCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGc                                                                                                                                  >  1:3887860/1‑150 (MQ=255)
                       aGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGGACGGGCGCATTATCACGCGCCCGGTCAGCttc                                                                                                                               <  2:8342546/150‑3 (MQ=38)
                         cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGGAGCTGGACGGGCGCATTATCACGCGCCCGGTCAGCttcaa                                                                                                                             <  1:7193612/150‑5 (MQ=37)
                         cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCgg                                                                                                                             <  1:7489816/150‑1 (MQ=255)
                         cTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCgg                                                                                                                             <  1:801098/150‑1 (MQ=255)
                          tCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgc                                                                                                                            >  2:3674183/1‑150 (MQ=255)
                             aCCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaat                                                                                                                         >  2:7665030/1‑150 (MQ=255)
                               ccAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaaggc                                                                                                                       >  2:6339374/1‑150 (MQ=255)
                                cAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgca                                                                                                                      >  1:4800435/1‑150 (MQ=255)
                                 aGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgcaa                                                                                                                     <  2:2510694/150‑1 (MQ=255)
                                 aGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgcaa                                                                                                                     >  2:2279779/1‑150 (MQ=255)
                                 aGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGgcaatgcaa                                                                                                                     >  2:4039038/1‑150 (MQ=255)
                                    ggcggGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAAtgt                                                                                                                  <  1:3647179/150‑1 (MQ=255)
                                      cggGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGc                                                                                                                <  2:4776337/150‑1 (MQ=255)
                                       gggCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCa                                                                                                               >  1:6244650/1‑150 (MQ=255)
                                             gCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCa                                                                                                         <  2:7027937/150‑1 (MQ=255)
                                             gCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAAGGTGCATGGCCa                                                                                                         >  2:416726/1‑150 (MQ=255)
                                                   cGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgc                                                                                                   <  1:2792803/150‑1 (MQ=255)
                                                    gACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcg                                                                                                  >  1:4448883/1‑150 (MQ=255)
                                                     aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCTGGgcgcgc                                                                                                 <  2:8077407/150‑1 (MQ=255)
                                                     aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGgcgcgc                                                                                                 <  2:2299510/150‑1 (MQ=255)
                                                     aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgc                                                                                                 <  2:5218329/150‑1 (MQ=255)
                                                     aCAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgc                                                                                                 <  1:688813/150‑1 (MQ=255)
                                                      cAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGTAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTcgcgcg                                                                                                <  1:6973752/150‑1 (MQ=255)
                                                              aTGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAgg                                                                                        >  1:7453220/1‑150 (MQ=255)
                                                                  gcTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCt                                                                                    <  1:2008511/150‑1 (MQ=255)
                                                                   cTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCTGGGCGCGCGCGACCTGGCCATg                                                                                   <  2:6809164/150‑1 (MQ=255)
                                                                     cGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCTGGGCGCGCGCGACCTGGCCATGca                                                                                 <  2:1602011/150‑2 (MQ=255)
                                                                      ggggTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCACGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTc                                                                                <  2:8714043/150‑1 (MQ=255)
                                                                         gTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCt                                                                             <  2:3447690/150‑1 (MQ=255)
                                                                          tGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTg                                                                            <  2:5884484/150‑1 (MQ=255)
                                                                             gcgcGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCt                                                                         <  1:8586221/150‑1 (MQ=255)
                                                                                cgGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTcc                                                                      <  1:3197387/150‑1 (MQ=255)
                                                                                      cAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                <  1:476967/150‑1 (MQ=255)
                                                                                      cAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagcca                                                                <  1:2498634/150‑1 (MQ=255)
                                                                                       aGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCACTGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCccagccag                                                               >  2:4638486/1‑150 (MQ=255)
                                                                                        gACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAgg                                                              <  1:4241079/150‑1 (MQ=255)
                                                                                          cAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCt                                                            <  1:1042642/150‑1 (MQ=255)
                                                                                           aaCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTg                                                           >  1:5063443/1‑150 (MQ=255)
                                                                                            aCGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTgg                                                          <  2:1759168/150‑1 (MQ=255)
                                                                                              gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                        <  1:560355/150‑1 (MQ=255)
                                                                                              gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                        >  2:2760564/1‑150 (MQ=255)
                                                                                              gTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGtt                                                        >  2:7535824/1‑150 (MQ=255)
                                                                                                 ggACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTc                                                     <  1:2438969/150‑1 (MQ=255)
                                                                                                 ggACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTc                                                     <  1:2198187/150‑1 (MQ=255)
                                                                                                 ggACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTc                                                     <  2:7881537/150‑1 (MQ=255)
                                                                                                    cTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTg                                                  >  1:7868539/1‑150 (MQ=255)
                                                                                                     tGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgc                                                 <  1:6652409/150‑1 (MQ=255)
                                                                                                      ggCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTgcg                                                <  1:8456527/150‑1 (MQ=255)
                                                                                                        cTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGca                                              >  1:7236669/1‑150 (MQ=255)
                                                                                                          gcgcTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACAGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCAAGCCAGGCTGGTTGTCCTGCGcaca                                            >  2:2989417/1‑150 (MQ=255)
                                                                                                           cgctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAg                                           >  1:5937161/1‑150 (MQ=255)
                                                                                                            gctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGa                                          <  2:657517/150‑1 (MQ=255)
                                                                                                            gctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGa                                          >  1:3444963/1‑150 (MQ=255)
                                                                                                             ctcgctGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGAt                                         <  1:3091802/150‑1 (MQ=255)
                                                                                                                  tGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCt                                    >  1:7103032/1‑150 (MQ=255)
                                                                                                                  tGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCt                                    >  1:6764997/1‑150 (MQ=255)
                                                                                                                  tGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCt                                    >  2:6996983/1‑150 (MQ=255)
                                                                                                                     gcCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCGgca                                 <  1:206181/150‑1 (MQ=255)
                                                                                                                         ggCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCAc                             >  2:187933/1‑150 (MQ=255)
                                                                                                                           ccATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCg                           >  1:2720988/1‑150 (MQ=255)
                                                                                                                            cATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCgg                          <  2:3120517/150‑1 (MQ=255)
                                                                                                                              tGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGa                        <  2:4807063/150‑1 (MQ=255)
                                                                                                                                  cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                    >  2:6741476/1‑150 (MQ=255)
                                                                                                                                  cTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTc                    >  2:3242555/1‑150 (MQ=255)
                                                                                                                                   ttGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                   >  1:3334111/1‑150 (MQ=255)
                                                                                                                                   ttGAAGCTGACCGGCCTCGTGATAATGCGCCCGTACAGCACCGGCAAATCAATGTGCATGGCCAGCTAGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCACAGCCAGGCAGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcg                   >  1:5272088/1‑150 (MQ=255)
                                                                                                                                     gAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTcgcg                 <  1:990048/150‑1 (MQ=255)
                                                                                                                                      aaGCTGAGCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGTCAATGCATTGTGCATGGCCAGGTCGCTCGCGCCCTGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcgg                <  2:4919111/150‑1 (MQ=255)
                                                                                                                                      aaGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcgg                <  1:1134800/150‑1 (MQ=255)
                                                                                                                                      aaGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcgg                <  2:5048806/150‑1 (MQ=255)
                                                                                                                                       aGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc               <  2:746695/150‑1 (MQ=255)
                                                                                                                                       aGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc               <  1:6019791/150‑1 (MQ=255)
                                                                                                                                       aGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggc               <  1:2150125/150‑1 (MQ=255)
                                                                                                                                        gCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCgcggcg              >  1:6459721/1‑150 (MQ=255)
                                                                                                                                        gCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGAGGTCgcggcg              >  2:7392316/1‑150 (MQ=255)
                                                                                                                                         cTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGa             <  2:5985509/150‑1 (MQ=255)
                                                                                                                                          tggacGGGCGCATTATCACGCGCCCGGTCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGaa            >  2:1256482/5‑150 (MQ=255)
                                                                                                                                           gACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGaaa           >  2:612861/1‑150 (MQ=255)
                                                                                                                                             acGGGCGCATTATCACGCGCCCGGTCAGCTTCGGCAATGCAATGTGCATGGCCAGGTCGCGCCCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACg         >  1:8653658/2‑150 (MQ=255)
                                                                                                                                               gggCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTgca                                 <  1:7411310/124‑1 (MQ=255)
                                                                                                                                               gggCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACggg       <  1:4159121/150‑1 (MQ=255)
                                                                                                                                                   gcgTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCa   >  2:6480549/1‑150 (MQ=255)
                                                                                                                                                   gcgTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCa   >  2:35241/1‑150 (MQ=255)
                                                                                                                                                    cgTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAg  <  1:1793748/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCTGGGCATTGGGCAGGCGGGCCAGCTCGACCCAGCGGCGGGCCAGCTCGGCGACAAAGTCCATGCGCTCGGGGTGGGCGCGGTAGCAGACAACGTCGGACTGGCTGCGCTCGCTGCGCCAGGCCATGTCCTTGAAGCTGACCGGGCGCGTGATAATGCGCCCGTCCAGCTCCGGCAATGCAATGTGCATGGCCAGGTCGCGCGCGCCCAGGCCCTGCTCGCTGGCTTCCCAGCCAGGCTGGTTGTCCTGCGCACAGATGGCCTGCAGCACCGGGATGTCGCGGCGAAACGGGCGCAG  >  NC_002947/3984495‑3984792

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