Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,374 G→C 7.8% intergenic (+64/‑79) PP_3338 → / → PP_3339 cytochrome bd‑type quinol oxidase/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,781,3740GC 7.8% 384.5 / 4.4 176intergenic (+64/‑79)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (79/75);  new base C (2/13);  total (81/95)
Fisher's exact test for biased strand distribution p-value = 5.72e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 185 total aligned reads displayed.

GCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGA  >  NC_002947/3781225‑3781522
                                                                                                                                                     |                                                                                                                                                    
gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCg                                                                                                                                                      >  1:4797818/1‑150 (MQ=255)
gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCg                                                                                                                                                      >  1:5623694/1‑150 (MQ=255)
gCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCg                                                                                                                                                      >  2:2429374/1‑150 (MQ=255)
   tCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGaga                                                                                                                                                   >  1:2395618/1‑148 (MQ=255)
       ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaaa                                                                                                                                               <  1:914171/150‑1 (MQ=255)
       ggAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACaaaaa                                                                                                                                               <  2:3644446/150‑1 (MQ=255)
          aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTTGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                            <  1:5505822/150‑1 (MQ=255)
          aGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgc                                                                                                                                            <  1:6778617/150‑1 (MQ=255)
           gCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGgcg                                                                                                                                           <  2:4705786/150‑1 (MQ=255)
              gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAAGgcgccc                                                                                                                                        >  1:2975187/1‑148 (MQ=255)
              gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGGCAAAAAGGCGcac                                                                                                                                        >  2:684482/1‑150 (MQ=255)
                aGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcacca                                                                                                                                      <  1:4152624/150‑1 (MQ=255)
                   gccgccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCAAAAACGCGcccccacg                                                                                                                                   >  1:2575766/1‑146 (MQ=255)
                      gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGGGCCCAAAAAGCGcccccccgacc                                                                                                                                >  1:755140/1‑143 (MQ=255)
                      gccCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCAGCGCCAAAAAGGCGCACCACCGAAc                                                                                                                                >  2:3318446/1‑150 (MQ=255)
                      gccCCAGAGCCAAGGCATCATGCAGGTAGGCACACTTTACATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATAGAGCCCTCGCCGTGGGGCGCCAGCGAAAAAAAGGCGCACCACGGAAc                                                                                                                                >  2:3318514/1‑150 (MQ=255)
                         ccAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGACAAAAAAGCGCACCACCGAACGTg                                                                                                                             >  1:4817363/1‑150 (MQ=255)
                          cAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAcaagt                                                                                                                            >  2:3664763/1‑146 (MQ=255)
                               ccATGGCTTCTTGGTGGTTGGAACCCTTTTCCTCGTGCCCTGTTCGGTTGGCTTTTCCTACTGTTCTTGCGGGCCGCCGGCATGCGCCGCCTTTGCGCCCGCGCGGGGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGc                                                                                                                       <  2:142486/150‑1 (MQ=255)
                               ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCCCGGAACATGTACAGc                                                                                                                       >  1:2805365/1‑150 (MQ=255)
                               ccAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGACAAAAAGGCGCACCACCGAACTTGTACAGc                                                                                                                       >  1:6778474/1‑150 (MQ=255)
                                       tCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGGCAAAAAAGCGCACCACGCAACGAGTACACCCCGCAACg                                                                                                               >  2:1302005/1‑150 (MQ=255)
                                            gtgttAGGCACCCTGTTCATCCGGCCCAGTTCGTTTGGCTGACCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                          <  1:5632333/146‑1 (MQ=255)
                                            cTGGTAGGCACACTTTACATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCCCAGCCATTGAGCCCACGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                          <  2:5651183/150‑1 (MQ=255)
                                             tGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGCACATGTACAGCCCGCCACGGATGCg                                                                                                         >  2:2758293/1‑150 (MQ=255)
                                              ggTAGGCACACTGTTCATCCTGCCCTGTTCGCTTGGCTGAGCCGAGTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                                        <  2:4052992/150‑1 (MQ=255)
                                                   gCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGCGCACCCCGGCCCTTGTAAAGCCCGCCACGGGTGCGGcggga                                                                                                   >  2:4037812/1‑146 (MQ=255)
                                                    cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGGACATTTAAAGCCCGCAAACGGTTCGGCGTAAt                                                                                                  >  2:1189352/1‑150 (MQ=255)
                                                    cacaCTGTTCAGCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCAGGCTGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                  <  2:2242925/150‑1 (MQ=255)
                                                     acacTGTTCTTCCGGCCCTTTTGGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCATTTGAGCCCTCGCCTTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATg                                                                                                 <  2:4755101/150‑1 (MQ=255)
                                                     acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATg                                                                                                 >  2:6768032/1‑150 (MQ=255)
                                                     acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGGGCACCACCGCACATGTAAAGCCCCCAACCGATGCGGCGTAATg                                                                                                 >  2:1679972/1‑150 (MQ=255)
                                                         tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAAGCGCACCACCGAACATGTACAGCCCCCAACCGATGCGGCTTAAAtcgcg                                                                                             >  2:958855/1‑150 (MQ=255)
                                                         tGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAAAATGTACAGCCCGCAACGGATGCGGCTTAATgcgcg                                                                                             >  2:3866792/1‑150 (MQ=255)
                                                          gTTCAGCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgcgc                                                                                            <  2:6731676/150‑1 (MQ=255)
                                                              acccggCCCTGTTCGCTGGGCTGAGCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCAATTGAGCCCGCGCGGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCata                                                                                        <  2:6396377/145‑1 (MQ=255)
                                                              aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGGCAAAAAGGCGCACCACGGGACTTGTACAGCCCGCAACCGATGCGGCTTAATGCGCGCCaat                                                                                        >  2:554356/1‑148 (MQ=255)
                                                                 ctccccAGTTCGCTTGGCTCACCCGACGGTACAGGATGGCCGCCGGCATGCCCAGCAATTGAGCCCTCGCCGGGGGGCGCCTGCCACAAAAAGCCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTc                                                                                     <  1:4569999/147‑1 (MQ=255)
                                                                  tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCCTATTCg                                                                                    >  2:6535870/1‑150 (MQ=255)
                                                                  tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACATGTACAGCCCGCAACGGATGCGGCGTAAAGCGCGCCATAttag                                                                                    >  2:2828871/1‑148 (MQ=255)
                                                                     cccgTTCGCTTGGCTGAGCCTACTTTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                                 <  2:4181602/147‑1 (MQ=255)
                                                                     ccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                                 <  2:6285895/150‑1 (MQ=255)
                                                                     ccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTc                                                                                 <  2:6285884/150‑1 (MQ=255)
                                                                      cAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGGCAAAAAAGCCCCCCACCGCACATGTACAGCCCGCCACCGATGCGGCGGAAAGCGCGCCCAATTAGCTCg                                                                                >  1:3873134/1‑150 (MQ=255)
                                                                       ttttCGCTTGGATTATCCTACTTTTCTTGTTGGCCGCGGGCTTGCGCAGCATTTGAGCCCTCGCGGTGGTGCGCCTGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCgg                                                                               <  1:3246548/148‑1 (MQ=255)
                                                                       aGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCGGCGACAAAAAAGCGCCCCACCGAACTTGTACAGCCCCCAACGGAGGCGGCGTAATGCGCGCCCTATTAGGTCgg                                                                               >  1:2699179/1‑150 (MQ=255)
                                                                        tttCGCTTGGCTGATCCGACTGTCATTGATGGCCGCGGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCggg                                                                              <  2:2938275/149‑1 (MQ=255)
                                                                         ttCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGGCAAAAAGGGGCCCCAACGAACGAGTACAGCCCGCAACCGATGCGGCTTAATTCGCGCCCTATTCGCACGGGt                                                                             >  2:722941/1‑150 (MQ=255)
                                                                           cGCTTGGCTTACCCGACTGTTCAGGCGGGCCGCCGGCAGGCGCAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAACGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGa                                                                           <  2:5770233/150‑1 (MQ=255)
                                                                            gCTGGGCTGAGCCGAGTGTTCATGCGGGCCGCCGGCTTGCGCAGCCATTGCGCCCGCGCGGTGGGGCGCCTGCCACAAACCGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGAc                                                                          <  2:2283442/150‑1 (MQ=255)
                                                                              ttGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAAAGGATGCGGCTTAATGCGCGCCAAATTCGCTCGGGTGACTc                                                                        >  2:526823/1‑150 (MQ=255)
                                                                              tgggCTGTTCCGCCTGTTCCTGAGGGCCGCCGGCTTGCGCAGCCTTTGCGCCCTCGCGGTGGGGCGCCTGCGCCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                        <  2:6428565/148‑1 (MQ=255)
                                                                                gggtgttCCGAGTGTTCCTGCGGGCCGCCGGCATGCGCCGCAATTGCCCCCGGGCGGTGGTGCGCCCGCCAAAAAACGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCAt                                                                      <  1:2540838/144‑1 (MQ=255)
                                                                                   ggATCCGACTGTTCATGCGGGCCGCCGGCATGCGCAGCCATTGAGCCCGCGCCGTGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCt                                                                   <  2:6036892/149‑1 (MQ=255)
                                                                                        cgccTGTTCATGTTGGCCGCCGGCAGGCCCAGCCATTGAGCCCGCGCCGGGGGGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                              <  2:325056/147‑1 (MQ=255)
                                                                                        cGACTGTTCTTGTGGGCCGCGGGCCGGCGCCGCCTTTGCCCCCCCGCGGGGGGGCGCCGGCCACAAACCGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                              <  1:3159216/150‑1 (MQ=255)
                                                                                              ttCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAAAATGAAAAGACCGCAACGGAAGCCCCTTAAAGCGCGCCAAATTAGCTCGGGTGACCAATGCATAcaccagcac                                                        >  2:6263151/1‑150 (MQ=255)
                                                                                               tCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCACCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacg                                                       >  1:1830847/1‑150 (MQ=255)
                                                                                                cATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCCAAAAGGCGCCCCCCCGGCCTAGTAAAGCCCCCAACGGCTGCGGCGTAATGAGCGCCAAAAAAGCTCGGGGGGCGCATGCGTTgacgggcaccg                                                      >  2:2391283/1‑148 (MQ=255)
                                                                                                 atgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCCCCACCGAACTTTTACAGCCCCCAACCGATGCGGCTTAATGCGCGCCAAATTCGCTCGGGTGTCTCATTCTTTaacgtgcacgtg                                                     >  1:4506306/1‑150 (MQ=255)
                                                                                                  tgatgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTTCGCCTTCGGCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCCGCTTTATGCGCGCCATATTCGCTCGGGTGACTCATGCGGTCacgtgcacgtgc                                                    >  2:5701305/1‑150 (MQ=255)
                                                                                                    atgGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACgtgcacgtgcac                                                  <  2:6652914/150‑1 (MQ=255)
                                                                                                        ccgccgGCTTGCGCAGCATTTGAGCCCTGGCGGTGTTGCGCCTGCTAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCt                                              <  1:1531081/150‑1 (MQ=255)
                                                                                                         cgccgGCATGCGCAGCAATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCtt                                             <  2:1384629/150‑1 (MQ=255)
                                                                                                          gcgggCTTGCGCTGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                            <  2:1710708/147‑1 (MQ=255)
                                                                                                          gccgGAATGCGCAGCATTTTCGCCCTCGCGGTGTTGCGCCTGCTACAAAATGGCGAACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                            <  1:5313260/150‑1 (MQ=255)
                                                                                                           ccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTg                                           <  2:5022369/150‑1 (MQ=255)
                                                                                                            gggCTTGCCCAGCATTTGCGCCCTGGCGGTGGTGCGCCTGCCAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                          <  1:3525141/149‑1 (MQ=255)
                                                                                                            cggctggcgcAGCCTTTGACCCCTCGCGGGGGTGCGCCTGCGAAAAACCGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                          <  1:2029030/144‑1 (MQ=255)
                                                                                                            cggctggcgcAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                          <  2:5211215/144‑1 (MQ=255)
                                                                                                            cgGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCTGCGACAAAAAGGGGCACCACCGGACATGTACAGCCCCCCACGGCTGCGGCTTAATGCCCGCCATATTCGCTCGGGTGCATCCGGCGTTCACGTGCACGTGCACATCTTTGt                                          >  2:4394348/1‑150 (MQ=255)
                                                                                                              gcttGCGCTGCATTTGAGCCCTGGCGGTGGGGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCt                                        <  1:2065832/147‑1 (MQ=255)
                                                                                                                 tGCGCCGCCTTTGCCCCGTCGCGGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGtt                                     <  1:6563500/150‑1 (MQ=255)
                                                                                                                  gcgcCGCCATTGAGCCCGCGCCGTGGTGCGCCTGCTAAAAAAAGGCGAACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                    <  2:1840756/150‑1 (MQ=255)
                                                                                                                  gcgcAGCCATTGCGCCCTCGCGGTGGTGCGCCTGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                    <  2:255003/150‑1 (MQ=255)
                                                                                                                  gcgcAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGCACCACCGAACGTGTACAGCCCCCAACCGATTCGGCTTAAAGCGCGCCAAATTCGCTCGGGTGGCTCATGCTTTAACGTGCACGTGCACAAATTTTTCTGTTc                                    >  2:2161838/1‑150 (MQ=255)
                                                                                                                      agcatttGCGCCCTCGCGGTGGTGCGCCTGCCACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcct                                <  1:3664763/145‑1 (MQ=255)
                                                                                                                         ctttGCGCCCTCGCGGTGGTGCGCCTGCCAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctgcc                             <  1:5391068/148‑1 (MQ=255)
                                                                                                                              cgCCCTCGCGGTGGTGCGCCTGCTAAAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                        <  1:1796599/149‑1 (MQ=255)
                                                                                                                              aGCCCTCGCCGGGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                        <  2:4610695/150‑1 (MQ=255)
                                                                                                                                cccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt                      <  2:3931456/150‑1 (MQ=255)
                                                                                                                                cccTCGCCGTGGTGCGCCTGCGAAAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGt                      <  1:5546765/150‑1 (MQ=255)
                                                                                                                                 ccTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTg                     <  2:3065100/150‑1 (MQ=255)
                                                                                                                                  cTCGCGGTGGTGCGCTTGCGACAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                    <  1:6467477/150‑1 (MQ=255)
                                                                                                                                   tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt                   <  2:2631755/150‑1 (MQ=255)
                                                                                                                                   tCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCt                   <  1:6297695/150‑1 (MQ=255)
                                                                                                                                        gtggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCa              <  1:5061231/150‑1 (MQ=255)
                                                                                                                                         tggtgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCaa             <  2:2716931/150‑1 (MQ=255)
                                                                                                                                         tggggCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGGTCGCCTGCCTCGCCGTGCTGTTCaa             <  2:1179505/150‑1 (MQ=255)
                                                                                                                                             gcgcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctg         >  2:833114/1‑150 (MQ=255)
                                                                                                                                                ccTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgctg      >  2:2347376/1‑150 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     >  2:1014248/1‑150 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     >  2:3817131/1‑150 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     <  2:6309929/150‑1 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     >  2:6107627/1‑150 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     >  1:4066798/1‑150 (MQ=255)
                                                                                                                                                 cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg     >  2:5027782/1‑150 (MQ=255)
                                                                                                                                                  tGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGc    >  2:6265515/1‑150 (MQ=255)
                                                                                                                                                  tGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGGTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGc    >  2:864625/1‑150 (MQ=255)
                                                                                                                                                     gACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                                >  2:5138206/1‑87 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCATCTGGGAAGCCGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGA  >  NC_002947/3781225‑3781522

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