Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,630,056 G→C 5.1% intergenic (+20/‑120) PP_3198 → / → PP_3199 2Fe‑2S ferredoxin/monooxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,630,0560GC 5.1% 792.4 / 22.9 235intergenic (+20/‑120)PP_3198/PP_31992Fe‑2S ferredoxin/monooxygenase
Reads supporting (aligned to +/- strand):  ref base G (122/101);  new base C (4/8);  total (126/109)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.26e-01

Only 100 of 241 total aligned reads displayed.

CGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTG  >  NC_002947/3629907‑3630205
                                                                                                                                                     |                                                                                                                                                      
cGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCGGCGCGATAg                                                                                                                                                        <  2:3549348/150‑1 (MQ=255)
  gcgcgcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGt                                                                                                                                                      <  1:2045308/150‑1 (MQ=255)
     cgcgcTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGc                                                                                                                                                   >  2:3007873/1‑150 (MQ=255)
          ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCCCAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                              >  1:895274/1‑150 (MQ=255)
          ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCa                                                                                                                                                  <  2:3766080/146‑1 (MQ=255)
          ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCa                                                                                                                                                  >  1:3766080/1‑146 (MQ=255)
          ttCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATccc                                                                                                                                              >  1:3265094/1‑150 (MQ=255)
           tCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGATGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                                             >  2:596564/1‑150 (MQ=255)
           tCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGATGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                                             >  2:5397362/1‑150 (MQ=255)
           tCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGATGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcg                                                                                                                                             >  2:5397509/1‑150 (MQ=255)
               gCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgac                                                                                                                                         >  1:2250158/1‑150 (MQ=255)
                 ccGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgacgc                                                                                                                                       <  1:4939884/150‑1 (MQ=255)
                  cGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCcgcgacgcg                                                                                                                                      >  1:4747777/1‑150 (MQ=255)
                     aGACGGCTCTAGCGTTCATGGCCCCTGCCGAGGCGACCCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTGGAGCAGCTCGATCTCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCTCGACGCGAtt                                                                                                                                   <  2:2865082/150‑1 (MQ=255)
                     aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                   <  1:4128650/150‑1 (MQ=255)
                     aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                   <  2:3899511/150‑1 (MQ=255)
                     aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGAtt                                                                                                                                   <  1:1381404/150‑1 (MQ=255)
                     aGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGAta                                                                                                                                   <  2:450634/150‑2 (MQ=37)
                      gACGGCACTTGCGTTCCCGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                                  <  2:2471489/150‑1 (MQ=255)
                        cGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGt                                                                                                                                <  2:3432763/150‑1 (MQ=255)
                         ggCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTc                                                                                                                               >  1:3021648/1‑150 (MQ=255)
                                  gTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCgg                                                                                                                      >  1:998693/1‑150 (MQ=255)
                                       cGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCAATCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGAt                                                                                                                                    <  2:1011685/131‑1 (MQ=25)
                                         gCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGc                                                                                                               <  1:3184325/150‑1 (MQ=255)
                                          ccccTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCg                                                                                                              <  1:4919029/150‑1 (MQ=255)
                                              tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                                          >  1:3180112/1‑150 (MQ=255)
                                              tGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGaa                                                                                                          <  2:5164590/150‑1 (MQ=255)
                                                ccGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTg                                                                                                                                  >  2:1525328/1‑124 (MQ=255)
                                                 cGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCt                                                                                                       <  1:663908/150‑1 (MQ=255)
                                                  gAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATAGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGcgc                                                                                                      <  2:4936770/150‑3 (MQ=255)
                                                       gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCtgtg                                                                                                    <  2:1763016/147‑1 (MQ=255)
                                                       gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAt                                                                                                 >  2:2321482/1‑150 (MQ=255)
                                                       gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAt                                                                                                 >  1:5231274/1‑150 (MQ=255)
                                                       gAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAt                                                                                                 >  2:80458/1‑150 (MQ=255)
                                                         aGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAg                                                                                               <  2:362973/150‑1 (MQ=255)
                                                             ttGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcac                                                                                           <  1:210541/150‑1 (MQ=255)
                                                             ttGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcac                                                                                           >  1:3116135/1‑150 (MQ=255)
                                                              tGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGcaca                                                                                          <  2:4676290/150‑1 (MQ=255)
                                                                   gCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCCTATGCACAgcttc                                                                                     <  2:1631426/150‑5 (MQ=255)
                                                                     aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCAt                                                                                                                             >  1:1190007/1‑108 (MQ=255)
                                                                     aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                   >  1:1770697/1‑150 (MQ=255)
                                                                     aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                   >  2:3227015/1‑150 (MQ=255)
                                                                     aCTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACg                                                                                   >  1:5107440/1‑150 (MQ=255)
                                                                       tCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgc                                                                                 >  2:2646147/1‑150 (MQ=255)
                                                                        cGATCACAAGGACGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgcc                                                                                <  1:4721824/150‑1 (MQ=255)
                                                                          aTCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGgccgc                                                                              >  2:5171195/1‑150 (MQ=255)
                                                                              caAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCg                                                                          >  2:4771246/1‑150 (MQ=255)
                                                                               aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                                        <  2:4720275/150‑1 (MQ=255)
                                                                               aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                                        <  1:4813610/150‑1 (MQ=255)
                                                                               aaGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGc                                                                        >  1:578077/1‑150 (MQ=255)
                                                                                 ggAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCAt                                                                      >  1:1706377/1‑150 (MQ=255)
                                                                                       gggcggCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACggg                                                                <  1:2471676/150‑1 (MQ=255)
                                                                                       gggcggCTGGCTGCTATTGCCGCTTGACCAGCTCGATGGCTGACCTATGGCGCCGCGCGTTTGGTCGAATCCCGCGACGCGATTGGTCATGGACGGGCGCCGGGGGTAAGCTGTGCATCGGCACAGGAAACGGCCGCCCCGCATGACggg                                                                <  1:2292725/150‑1 (MQ=255)
                                                                                        ggcggcTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACgggg                                                               <  2:859957/150‑1 (MQ=255)
                                                                                         gcggcgGGCTGCTATGGCCGCTTGAGCAGCTCGATCGCGGCCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATCGGCACAGGAAACGGCCGCCCCGCATGACGGGgc                                                              <  2:18382/150‑1 (MQ=255)
                                                                                         gcggcTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGgc                                                              <  2:419596/150‑1 (MQ=255)
                                                                                            gctggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGTCATGCGGGgcggc                                                           <  2:4082850/150‑3 (MQ=255)
                                                                                            gctggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGca                                                           <  2:3945875/150‑1 (MQ=255)
                                                                                              tggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTCCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGcaca                                                         >  1:1880336/1‑150 (MQ=255)
                                                                                               ggctgCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGTCATGCGGGgcggccgt                                                        <  1:3698277/150‑6 (MQ=255)
                                                                                                gctgctATTGCCGCTTGAGCAGCTCGATCGCGGACCTAGCGCGCCGCGCGATAGGTCGCACCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGTCATGCGGGgcggccgtt                                                       <  2:218464/150‑7 (MQ=255)
                                                                                                  tgctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCAcaacaa                                                     <  1:2211972/150‑1 (MQ=255)
                                                                                                    ctATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGa                                                   >  2:572626/1‑150 (MQ=255)
                                                                                                     tATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                  <  2:1257103/150‑1 (MQ=255)
                                                                                                     tATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCACATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                  <  1:4764536/150‑1 (MQ=255)
                                                                                                       ttGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACa                                                >  1:2086896/1‑150 (MQ=255)
                                                                                                          ccGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTc                                                                                                                               >  1:1479326/1‑69 (MQ=255)
                                                                                                          ccGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCa                                                                                                  >  2:3311653/1‑98 (MQ=255)
                                                                                                               tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGTTTGGTCATGGACGGGCGCTGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                        <  1:1053348/150‑1 (MQ=255)
                                                                                                               tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                        <  2:1318334/150‑1 (MQ=255)
                                                                                                               tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                        <  2:3166510/150‑1 (MQ=255)
                                                                                                               tGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGAATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTga                                        <  1:1915594/150‑1 (MQ=255)
                                                                                                                   cagcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACa                                    <  1:3935499/150‑1 (MQ=255)
                                                                                                                    agcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACAt                                   >  1:4370716/1‑150 (MQ=255)
                                                                                                                     gcTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATg                                  >  1:33265/1‑150 (MQ=255)
                                                                                                                     acggcaTCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGTTTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATg                                  <  1:259694/145‑1 (MQ=255)
                                                                                                                      cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGc                                 <  1:3756158/150‑1 (MQ=255)
                                                                                                                      cTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGc                                 >  2:3441479/1‑150 (MQ=255)
                                                                                                                       tCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAGAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGcc                                >  1:93484/1‑150 (MQ=255)
                                                                                                                         gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACACGAACAAAACGTGAGAcgcctctg                              >  2:5345832/1‑143 (MQ=255)
                                                                                                                         gATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCAt                              >  1:3062056/1‑150 (MQ=255)
                                                                                                                         gATCGCTGACCTATCGCGCCGCGCGAGAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCAGAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGAcgcctctg                              >  2:4980356/1‑143 (MQ=255)
                                                                                                                          aTCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATg                             >  1:2778542/1‑150 (MQ=255)
                                                                                                                             gCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGGGAGACATGCCATGAAt                          >  1:4459833/1‑150 (MQ=255)
                                                                                                                               tGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGGGAGACATGCCATGAATAc                        >  1:1010289/1‑150 (MQ=255)
                                                                                                                                 aCCTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtg                      >  2:4593767/1‑150 (MQ=255)
                                                                                                                                 aCCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtg                      <  1:3645717/150‑1 (MQ=255)
                                                                                                                                  ccTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt                     >  2:1158576/1‑150 (MQ=255)
                                                                                                                                  ccTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt                     <  1:3472340/150‑1 (MQ=255)
                                                                                                                                  ccTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgt                     <  2:2250244/150‑1 (MQ=255)
                                                                                                                                   cTATCGCGCGGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACtgtg                    >  2:2463260/1‑150 (MQ=255)
                                                                                                                                   cTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACACCAAGAACAAAACGTGAGACATGCCATGAATACtgtg                    >  2:1184860/1‑150 (MQ=255)
                                                                                                                                    tATCGCGACGCGCGATAGGTCGAATCCCGCGACGCGATTGGTCAGGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACCTGCCATGAATACTGTGa                   >  2:4909748/1‑150 (MQ=255)
                                                                                                                                    aaTCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                  >  1:3777937/2‑119 (MQ=25)
                                                                                                                                    aaTCGCGTCGCGGGATGCGACCTATCGCGCGGCGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGaa                                                  <  2:3777936/118‑1 (MQ=25)
                                                                                                                                     aTCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaa                  >  1:5074834/1‑150 (MQ=255)
                                                                                                                                        gcgccgcgcGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGaacaa               >  1:2513328/1‑150 (MQ=255)
                                                                                                                                          gccgcgcGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCCCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACg             >  1:4660646/1‑150 (MQ=255)
                                                                                                                                             gcgcgATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACgtgt          <  1:1258937/150‑1 (MQ=255)
                                                                                                                                                     ggTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTg  >  2:3690503/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CGGCGCGCGCTTCACGCCCGAAGACGGCACTTGCGTTCACGGCCCCTGCCGAGGCGAAGCCTTGGAGGCACTCGATCACAAGGAAGCGGGCGGCTGGCTGCTATTGCCGCTTGAGCAGCTCGATCGCTGACCTATCGCGCCGCGCGATAGGTCGCATCCCGCGACGCGATTGGTCATGGACCGGCGCAGGCGGGAAGCTGTGCATAGGCACAGGAAACGGCCGCCCCGCATGACGGGGCGCACAACAAGAACAAAACGTGAGACATGCCATGAATACTGTGAACAACGTGTTGATCGTG  >  NC_002947/3629907‑3630205

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