Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,867,403 T→G 7.1% intergenic (‑56/+76) gllP ← / ← glllA porin‑like protein/gallate dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,867,4030TG 7.1% 659.9 / 7.1 220intergenic (‑56/+76)gllP/glllAporin‑like protein/gallate dioxygenase
Reads supporting (aligned to +/- strand):  ref base T (123/76);  new base G (2/19);  total (125/95)
Fisher's exact test for biased strand distribution p-value = 3.27e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 226 total aligned reads displayed.

CATCGTCGACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCT  >  NC_002947/2867254‑2867552
                                                                                                                                                      |                                                                                                                                                     
cATCGTCGACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCt                                                                                                                                                       >  2:5052662/1‑150 (MQ=255)
 aTCGTCGACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTGTTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTg                                                                                                                                                      <  1:2029903/150‑1 (MQ=255)
  tcgtcgACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTgg                                                                                                                                                     <  1:3480762/150‑1 (MQ=255)
    gtcgACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGCCGATGCGCCGCTggcg                                                                                                                                                   <  1:6453938/150‑3 (MQ=255)
     tcgACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCa                                                                                                                                                  >  2:4807695/1‑150 (MQ=255)
      cgacgaAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGGCGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAg                                                                                                                                                 >  1:7419126/1‑150 (MQ=255)
      cgacgaAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGCCGATGCGCCGCTggcggg                                                                                                                                                 <  2:6177927/150‑5 (MQ=255)
       gacgaAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAgg                                                                                                                                                <  2:8927606/150‑1 (MQ=255)
       gacgaAGCCTTCGGCCAGGGCGTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGGGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTGATTCTCCCGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAgg                                                                                                                                                <  2:6607166/150‑1 (MQ=255)
        acgaAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGt                                                                                                                                               >  2:7200579/1‑150 (MQ=255)
           aaGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCAGCGGCGCGGGGCGGGTTGa                                                                                                                                            >  1:8043184/1‑150 (MQ=255)
            cgCCTTCGGCCCGGGCTTGCAGTGGCAGGCGGGCTCCCAGTGCGAATGTGCGGACGCGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGaa                                                                                                                                           <  2:5864479/149‑1 (MQ=255)
              ccTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATc                                                                                                                                         <  1:8961709/150‑1 (MQ=255)
               cTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCa                                                                                                                                        >  2:5073537/1‑150 (MQ=255)
               cTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCa                                                                                                                                        >  1:4115279/1‑150 (MQ=255)
               cTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCAGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCa                                                                                                                                        >  2:4400616/1‑150 (MQ=255)
                ttCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAg                                                                                                                                       >  2:4669938/1‑150 (MQ=255)
                ttCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAg                                                                                                                                       >  1:150412/1‑150 (MQ=255)
                 tCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAgg                                                                                                                                      <  2:6347171/150‑1 (MQ=255)
                 tCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAgg                                                                                                                                      <  1:8612820/150‑1 (MQ=255)
                   ggCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAgggg                                                                                                                                    <  2:5688700/150‑1 (MQ=255)
                   ggCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCGGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAgggg                                                                                                                                    <  2:1708432/150‑1 (MQ=255)
                      cAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCg                                                                                                                                 >  1:6882788/1‑150 (MQ=255)
                      cAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCg                                                                                                                                 >  1:1206312/1‑150 (MQ=255)
                      cAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCg                                                                                                                                 >  1:1206934/1‑150 (MQ=255)
                        gggCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGt                                                                                                                               <  2:884625/150‑1 (MQ=255)
                           cTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTggg                                                                                                                            >  1:4955221/1‑150 (MQ=255)
                           cTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTggg                                                                                                                            <  2:3632583/150‑1 (MQ=255)
                           cTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTggg                                                                                                                            <  1:6027672/150‑1 (MQ=255)
                           cTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTggg                                                                                                                            <  2:5797032/150‑1 (MQ=255)
                            ttGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGTGGTgggg                                                                                                                           >  1:1902115/1‑150 (MQ=255)
                            ttGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTgggg                                                                                                                           >  1:5246080/1‑150 (MQ=255)
                            ttGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTgggg                                                                                                                           >  2:1980677/1‑150 (MQ=255)
                            ttGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTgggg                                                                                                                           >  1:3967701/1‑150 (MQ=255)
                            ttGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTgggg                                                                                                                           >  1:3967364/1‑150 (MQ=255)
                              gCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTTAATCAGGGGGCGGTGGGGGt                                                                                                                         <  2:2017037/150‑1 (MQ=255)
                               cAGTGGCAGGCTGGCTGCCAGTGCGCATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTg                                                                                                                        >  2:6294482/1‑150 (MQ=255)
                                aGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTgg                                                                                                                       >  2:3224893/1‑150 (MQ=255)
                                aGTGGCAGGCAGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGGGGGGGTgg                                                                                                                       >  2:8545366/1‑150 (MQ=255)
                                 gTGGCTGGCTGGCTGCCAGTGCGATTGTGGGGACGTGCGGGTACTGGGTGCGGCATTTCTTTATTTTTATTCTCCTGGGTTGGGCAATGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGTTTGAATCAGGGGGCGGTGGGGGTggc                                                                                                                      <  1:2556749/150‑1 (MQ=255)
                                  tGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCCGGTATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTggcg                                                                                                                     <  2:7992322/150‑1 (MQ=255)
                                  tGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTggcg                                                                                                                     >  1:2621495/1‑150 (MQ=255)
                                    gCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTggcggc                                                                                                                   >  2:3686462/1‑150 (MQ=255)
                                      tggctggctGCCAGTGCGAATGTGGGGACGAGCGGGTACTGGGTGCGGCATTTAATTATTTTTATTCTCCAGGGATGGGAAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGgcgc                                                                                                                 <  2:3087976/149‑1 (MQ=255)
                                         ctggctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCggg                                                                                                              >  2:2471397/1‑150 (MQ=255)
                                            gctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcg                                                                                                           >  1:3442742/1‑150 (MQ=255)
                                            gctgCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGAcg                                                                                                           >  2:1334906/1‑150 (MQ=255)
                                                   gggctattgtgGGCACGCGCGGGTACTGGGTGGGGCATTTCTTTATTTTTTTCCCCCTGGGTTGGGCAACGCGCGGCCCGCCCGGGGCGCATGGGCGCTGGGCTGGTTGAAACTGGGGGGGGGGGGGGGGGCGGCGCGGGACGCGCCGTc                                                                                                    <  2:4918489/143‑1 (MQ=255)
                                                    tGCGATTGTGGGGACGAGCGGGTACAGGGTGCGGCATTTAATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGGTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCGCca                                                                                                   <  1:7944286/150‑1 (MQ=255)
                                                    tGCGAATGTGGGGACGAGCGGTTAAAGGTTGCGGCATTTCATTATTTTTATTCTCCAGGTATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCGCca                                                                                                   <  2:1554528/150‑1 (MQ=255)
                                                    tGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTca                                                                                                   >  2:5488307/1‑150 (MQ=255)
                                                    tGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTca                                                                                                   >  1:832172/1‑150 (MQ=255)
                                                     gCGAATGTGGGGACGCGCGGGTACCGGGTGCGGAATTTCATTATTTTTATTCTCCCGGTATGGGCAACGCGCTGCCCGCCTGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGACGGCGCGTCCCGCGCCGCcac                                                                                                  <  2:4551666/150‑1 (MQ=255)
                                                       gAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCacac                                                                                                >  1:2904972/1‑150 (MQ=255)
                                                       gAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCacac                                                                                                >  1:1836408/1‑150 (MQ=255)
                                                          tgtgGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGAt                                                                                             >  1:1922513/1‑150 (MQ=255)
                                                             ggggCCGCGGGGGTACTGGGTGCGGAATTCCATTATTTTTTTCCCCCCGGTTGGGGCAACGGGCGGCCCGCCCGGGGCGCAGCGGCGCTGGGCGGGTTGAATCCGGGGGCGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCaa                                                                                          <  2:2053255/150‑1 (MQ=255)
                                                             ggggACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCaa                                                                                          >  1:4561644/1‑150 (MQ=255)
                                                             ggggACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCaa                                                                                          >  2:4084386/1‑150 (MQ=255)
                                                              gggACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAg                                                                                         >  1:744716/1‑150 (MQ=255)
                                                              gggACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGGTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAg                                                                                         >  2:648789/1‑150 (MQ=255)
                                                               ggccgcgCGGGTACCGGGTGCGGCTTTACATTATTTTTTTCCCCCCGGCATGGGCAACGCGGGGCCCGCCCGGGGCGCAGCGGCGCTGGGCCGGTTGAACCTGGGGGGGGTGGGGGGGGCGGCGCGGGACGCGCCGTCACACGATCAAGc                                                                                        <  2:7499236/144‑1 (MQ=255)
                                                                gACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGTATCAGGGGGCGGTGGGGGTGGCGGGGCGGGACGCGCCGTCACTCGATCAAGCg                                                                                       >  2:2250445/1‑150 (MQ=255)
                                                                  cgtgggggAACGGGGGGCGGCTTTTCATTTTTTTTATCCCCCGGGGCGGGCCACCGCGGTGCCCGCCCGGGGCCCATCGGCGGTGGGCTGTTTTATCGGGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCggg                                                                                     <  2:1459039/145‑1 (MQ=255)
                                                                   gAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGt                                                                                    >  1:5195505/1‑150 (MQ=255)
                                                                    accGGGTACCGGGCGCGGCTGTCCATTATGTTCCTGCCCCGGGGCTGGGCACCGCGCGGCCCGCCTGCGGCGCAGGGGCGGTGGGCCGGTTGAACCCGGGGGCGGTGGGGGTGGCGGCGGGGGACGCGCCGTCACACGATCAAGCGGGtt                                                                                   <  2:7130433/148‑1 (MQ=255)
                                                                      ggggAACCGGGGGCGGCATTTCATTATTTTTATGCCCCAGGGTTGGGAAACGCGCGGCCCGCCCGGGGCGCAGCGGCGCGGGGCTGTTTGAACCTGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGc                                                                                 <  2:5439997/149‑1 (MQ=255)
                                                                        gttACTGGGTGGGCCTTTTCTTTTTTTTTTTTCCCCCGGGTTGGGCAACGCGCGGCCCGCCTGCGGCGCATGGGCGGTGGGATGTTTTAATCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCAt                                                                               <  2:1922376/148‑1 (MQ=255)
                                                                           aCAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCgg                                                                            >  2:8318418/1‑150 (MQ=255)
                                                                               ggTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAg                                                                        >  1:5681027/1‑150 (MQ=255)
                                                                                   cGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGCATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGGCAGTTgg                                                                    >  2:2469458/1‑150 (MQ=255)
                                                                                        tttCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCt                                                               >  1:4466270/1‑150 (MQ=255)
                                                                                        ttcatttattTTTTTACCCCTGGGTTGGGAAATGCGCTCCCCGCCCGCGGGGCATGGGGGGGGGGCAGTTTTAACCGGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCt                                                               <  1:6788302/145‑1 (MQ=255)
                                                                                          actttttttttATGCGCCCGGGCTGGGCAACGGGCGGCCCGCCCGGGGGGCTGGGGGGCGGGGCTGGTTGACTCTGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCGCACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCttt                                                             <  2:6874660/144‑1 (MQ=255)
                                                                                            attattTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgc                                                           >  2:1783724/1‑150 (MQ=255)
                                                                                            attattTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACAATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgc                                                           >  1:1011486/1‑150 (MQ=255)
                                                                                             ttattTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgcc                                                          >  2:8336947/1‑150 (MQ=255)
                                                                                               tattttttTCCCCACGGTATGGCAAAAGCGCTGCCCGCCAGGGGCGAAACGGCGCTGGGAATCTTTAACCAGGGGGCGGGGGGGTTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgcca                                                         <  2:3384527/148‑1 (MQ=255)
                                                                                                tttttATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTgccagc                                                       >  1:243184/1‑150 (MQ=255)
                                                                                                  tttATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAATCGGGTTGCATCGGTCAGTTGGGCGCTTTgccagcca                                                     >  1:6225517/1‑150 (MQ=255)
                                                                                                     aTTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCg                                                  >  1:3851493/1‑150 (MQ=255)
                                                                                                         tcCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTa                                              >  2:1846589/1‑150 (MQ=255)
                                                                                                         tcCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTa                                              >  1:2778373/1‑150 (MQ=255)
                                                                                                          ccAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAg                                             >  1:9123743/1‑150 (MQ=255)
                                                                                                            ccgggTTGGGAAACGCGCTGCCCGCCCGGGGCGCAGCGGCGCGGGGCAGTTTTAATCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGc                                            <  2:2323718/148‑1 (MQ=255)
                                                                                                               ggATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGaa                                         >  2:6000515/1‑150 (MQ=255)
                                                                                                                     gCAACGCGCTGCCCGCCCGCGGGGCACCGGCGGGGGGCAGGTTCAACCCGGGGGCGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgct                                   <  2:8197200/150‑1 (MQ=255)
                                                                                                                     gCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgct                                   >  1:962535/1‑150 (MQ=255)
                                                                                                                      cAACGCGCTGCCCGCCCGCGGCGCAGCGGCGCGGGGCTGGTTGAAGCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgctg                                  <  2:1857275/150‑1 (MQ=255)
                                                                                                                      cAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGGCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACctgctg                                  >  2:789438/1‑150 (MQ=255)
                                                                                                                       aaaGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCCGCCACCGAGTAGCGAACCAGCTgg                                 >  1:2278162/1‑150 (MQ=255)
                                                                                                                             cggcccgcccGCGGGGCATCGGCGGTGGGCTGGTTGAAGCTGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACgg                           <  2:7668515/148‑1 (MQ=255)
                                                                                                                                cccGCCTGGGGCGCAGCGGGGGGGGGCTGTTTTATCCTGGGGGGGGTGGGGGTGGGGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTc                        <  1:4696539/150‑1 (MQ=255)
                                                                                                                                cccGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGGTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTc                        >  2:1082298/1‑150 (MQ=255)
                                                                                                                                  cGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCtt                      >  2:5746970/1‑150 (MQ=255)
                                                                                                                                   gcccggggcgCTTGGGCGGGGGGCTGGTTGACTCTGGGGGGGGGGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTc                     <  2:315635/144‑1 (MQ=255)
                                                                                                                                   gCCCGCGGCGCATCGGCGCTGGGCTGTTTGAAGCTGGGGGGGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTc                     <  2:5157108/150‑1 (MQ=255)
                                                                                                                                   gCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTc                     >  1:1048655/1‑150 (MQ=255)
                                                                                                                                            gcATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTcc            >  2:3790489/1‑150 (MQ=255)
                                                                                                                                                    ggtGGGCTGTTTTATTCTGGGGGGGGTGGGGTTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGct    <  2:8827374/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
CATCGTCGACGAAGCCTTCGGCCAGGGCTTGCAGTGGCAGGCTGGCTGCCAGTGCGAATGTGGGGACGAGCGGGTACAGGGTGCGGCATTTCATTATTTTTATTCTCCAGGGATGGGCAAAGCGCTGCCCGCCAGCGGCGCATCGGCGCTGGGCAGGTTGAATCAGGGGGCGGTGGGGGTGGCGGCGCGGGACGCGCCGTCACACGATCAAGCGGGTTGCATCGGTCAGTTGGGCGCTTTGCCAGCCACCGAGTAGCGAACCTGCTGGTTACGGGTCTTCATGAAGTCCTCGATGCTCT  >  NC_002947/2867254‑2867552

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: