Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,781,387 A→C 11.0% intergenic (+77/‑66) 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,3870AC11.0% 426.9 / 6.6 179intergenic (+77/‑66)PP_3338/PP_3339cytochrome bd‑type quinol oxidase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (46/107);  new base C (11/13);  total (59/120)
Fisher's exact test for biased strand distribution p-value = 1.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 194 total aligned reads displayed.

CGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACCGT  >  NC_002947/3781238‑3781536
                                                                                                                                                     |                                                                                                                                                     
 gCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAAGGCGcac                                                                                                                                                      >  2:4314100/1‑150 (MQ=255)
    gcgcCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGcaccac                                                                                                                                                   <  1:4813928/150‑1 (MQ=255)
      gccgccCCAGAGCCAAGGCTTCTTGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACgg                                                                                                                                                 <  1:1856289/150‑1 (MQ=255)
          ccccAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGGCTGCGCCAAAAAGGCGCGCCCCggccca                                                                                                                                             >  2:4016303/1‑146 (MQ=255)
             cAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGaaatagt                                                                                                                                          >  1:7301279/1‑145 (MQ=255)
              agagCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTa                                                                                                                                         <  2:6157943/150‑1 (MQ=255)
                cgCCATGGCTTATTGGTGGTTGGCACACTGTTCATCCGGCCCTGTTGGCTTGGCTGATCCGACTGTTCATGTTGGCCGCCGGCATGCGCTGCATTTGAGCCCTCGCGGTGGTGCGCCTGCTACAAAACGGCGCACCACGGAACTTGTACa                                                                                                                                       <  2:5745576/149‑1 (MQ=255)
                  caaaGGCTTCATGCTGTTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCTACTGTTCATGATGGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGc                                                                                                                                     <  2:4178915/148‑1 (MQ=255)
                    aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGGGGCAAAATGGGGCCCCACGGAACTTGTTCAGccc                                                                                                                                   >  2:1682182/1‑150 (MQ=255)
                    aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATGTACAGccc                                                                                                                                   >  2:1091762/1‑150 (MQ=255)
                    aaGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCAGCGACAAAAATGCGCACCACCGAACATGTAAAAccc                                                                                                                                   >  1:4830261/1‑150 (MQ=255)
                         ttCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGCGCCCCCCCGCACATGTAAAGccccccac                                                                                                                              >  1:5912432/1‑145 (MQ=255)
                               cTGGTTGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCTTGTTGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                                        <  2:5387268/150‑1 (MQ=255)
                               cTGGTCGGCACACTGTTCCTCCTGCCCTGTTGGCTTGGCTTTTCCGAGTTTTCTTGTGGGCCGCCGGCTGGCGCAGCATTGGCGCCCTGGCGTTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGc                                                                                                                        <  1:3560487/150‑1 (MQ=255)
                                 ggTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCgg                                                                                                                      <  2:7403259/150‑1 (MQ=255)
                                      gCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGCCAAAAAGGCGCACCCCGGAGCTTGTACAGCCCGCCACGGATGCGGCGTaa                                                                                                                 >  2:3295428/1‑150 (MQ=255)
                                      gCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAAGGCGCACCACGGAACATGTACAGCCCGCAACCGATGCGGCTTaa                                                                                                                 >  2:2928276/1‑150 (MQ=255)
                                       cacaCTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCTTGTTGGCCGCCGGATTGCGCCGCCATTGCGCCCCCGCGGGGGTGCGCCTGCTACAAAATGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAAt                                                                                                                <  1:2299921/150‑1 (MQ=255)
                                        acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACCGAACATGTACAGCCCGCAACCGATGCGGCTTAATg                                                                                                               >  2:4814878/1‑150 (MQ=255)
                                        acacTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGca                                                                                                                                                       >  2:895880/1‑110 (MQ=255)
                                         cacTGTTCATCCTGCCCTGTTCGATTGGCTGATCCGACTGTTCATGATGGCCGCCGGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgc                                                                                                              <  2:1823422/150‑1 (MQ=255)
                                          acTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAATGCGCACCACCGAACATGTAAAGCCCCCAACCGATTCGGCTTAAAgcg                                                                                                             >  2:153898/1‑150 (MQ=255)
                                           cggTTCATCCGGCCCAGTTCGGTGGGCTGATCCGACTGTTCTTTATGGCCGCGGGCATGCGCTGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                                            <  2:6516643/148‑1 (MQ=255)
                                                cATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCGGCGCCACAAAAGCCCCCCCCCGCACGTGTTAAGCCCGCCACCGATGCGGCGGAAATCGCGccca                                                                                                       >  1:2188379/1‑148 (MQ=255)
                                                 aTCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGCCAAAAAAGCGCCCCCCCGAACTTTTAAAGCCCGCCACCGGTGCGGCTTAATGCGCGcccaa                                                                                                      >  1:5477360/1‑147 (MQ=255)
                                                   ccTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGGACCACGGGACGTTTAAAGCCCCCCACCGAAGCGGCGTAATGCGCGcccaaat                                                                                                    >  2:4642401/1‑145 (MQ=255)
                                                   ccTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCCACCGATGCGGCTTAATGCGCGCCCTAtt                                                                                                    >  2:4732963/1‑150 (MQ=255)
                                                     tGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGGCTGGGTGAAAAATGGGCGCCACGGGAATTTTAAAGGCCGCAAAGGTTTGGGGTTTTTGGGCGGCATAtttg                                                                                                  >  2:971041/1‑148 (MQ=255)
                                                       cccttttCGCTTGGCTTATCCGAGTGTTCATGTTGGCCGCCGGCTTGCGCCGCCTTTGAGCCCTCGCGGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCt                                                                                                <  1:3020306/145‑1 (MQ=255)
                                                       ccctcttCCCTTGGCTTACCCGACTGTTAATTATGCCCCCCGCAATCCCAAGCAATTGACCCCTCGCCGTGTTCCCCCTCCGAAAAAAACCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCt                                                                                                <  2:2781524/145‑1 (MQ=255)
                                                       cccAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACCGAACATTTTAAGCCCCCCACGGATGCGGCGAAATGCGCGCCCAATTCGca                                                                                                >  1:6743694/1‑149 (MQ=255)
                                                         cAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCTGCGGCAAAAAGGCGCACCCCCGAACATGTAAAGCCCCCCACCCCTGCGGCGTAATGCGCGCCAGATTAGcgag                                                                                              >  1:5255696/1‑147 (MQ=255)
                                                            ttCGTTTGGCTTATCCGACTTTTATTCATGGCCGCCGGCATGCCCACCAATTGACCCCTCCCCGTGGGGCGCATCCCAAAAAAAGCCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGt                                                                                           <  1:1216372/150‑1 (MQ=255)
                                                            ttCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGGGCCAGCGGCCAAAAAGCGCACCCCCGCACATGTAAAGCCCCCAAAGGATGCGGGGTTATTGGCGCCAAATTGGcgagggg                                                                                           >  2:3135883/1‑144 (MQ=255)
                                                                 ttGGCTGATCCGACTGTTCATTTTGGCCGCCGGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                                      <  2:1360417/150‑1 (MQ=255)
                                                                 gtGGCTGATCCTACTGTTCATTTTGGCCGCCGGAATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTc                                                                                      <  2:5998767/149‑1 (MQ=255)
                                                                         tCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTc                                                                              <  2:3612093/150‑1 (MQ=255)
                                                                           cGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGCGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcac                                                                            <  2:3224348/150‑1 (MQ=255)
                                                                           cGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCGGCGACAAAAATGCGCACCACCGCACGTGTAAAGCCCGCAAAGGATGCGGCTTAATGCGCGCCAAATTCGCTCGGGTGACTCAAGCTTTaac                                                                            >  2:3364999/1‑150 (MQ=255)
                                                                             agtGTTCATGATGGCCGCCGGCTTGCGCCGCATTGGCGCCCTCGCCTTGGTGCGCCTGCTAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgt                                                                          <  1:550509/148‑1 (MQ=255)
                                                                               tGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGGGCGCCGGCGCCAAAAAGGCGCCCCCCGGCACTTGTACAGCCCCCCCCCGCTGCGGCGTAAAGAGCCCCCCATTAGCGCGGGTGCCACCTGCTTTaacctgc                                                                        >  2:5103865/1‑150 (MQ=255)
                                                                                    ttgttgGCCGCGGGCATGCGCCGCATTTGCGCCCTCGCGGGGGGGCGCCGCCCACAAAAAGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                                   <  1:1334256/146‑1 (MQ=255)
                                                                                    atgatgGCCGCCGGCATGCGCAGCCTTTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAACGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtg                                                                   <  2:5095761/150‑1 (MQ=255)
                                                                                     tgtgggCCGCGGGCATGCGCACCCATTGCCCCCGCGCGGTGGTGCGCCTGCTAAAAAAAGGCCCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCacgtgcacgtgc                                                                  <  1:1589008/146‑1 (MQ=255)
                                                                                      gtgggCCGCCGGCTGGCGCCGCCTTTGCGCCCTCGCGGTGGTGCGCCTGCCAAAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCAcgtgcacgtgca                                                                 <  1:6939115/147‑1 (MQ=255)
                                                                                         ggCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCGACAAAAAGGCGCACCACGGAAAGTTTAAAGCTAGCCAACGAAGCGGCTTAATTCGCGCCCAATTAGCGCGGGTGACTCAAGCTTTaacgtgcacgtgagttt                                                              >  2:3370813/1‑145 (MQ=255)
                                                                                           ccgccgTCTTGCGCAGATTTTGCGCCCGCGCGGTGGTGCGCCTGCCACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCt                                                            <  1:3364999/150‑1 (MQ=255)
                                                                                           ccgccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTcacgtgcacgtgaaagcat                                                            <  2:5114041/150‑8 (MQ=255)
                                                                                            cgccgGCATGCGCAGCAATTGAGCCCTCGCCGTGGTGCGCGTGCGAAAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCtt                                                           <  2:2722357/150‑1 (MQ=255)
                                                                                             gcgggCTTGCGGAGCCTTTGAGCCCTCGCCTTGGTGGGCGTGCGACAAATTGGCGCACCACGGAACTTGTACAGACCGCAACGGATGCGGCTTAATGCGCGCCTTATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:4927051/147‑1 (MQ=255)
                                                                                             gccgccAGGCGCAGCCATGGCCCCCGCGCGGTGGTGCGCCCGCGACAACACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  1:6577512/150‑1 (MQ=255)
                                                                                             gccgGCTTGCGCAGCCATTGTGCCCTCGCCGTGGTGCGCATGCGACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:521521/150‑1 (MQ=255)
                                                                                             gccgGCATGCGCCGCCATTGCGCCCTCGCCGTGGGGCGCGGGTGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:2403507/150‑1 (MQ=255)
                                                                                             gccgGCATGCGCCGCCATTGAGCCCTCGCCGTGGTGCGCCTGCTACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:4665823/150‑1 (MQ=255)
                                                                                             gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:7295579/150‑1 (MQ=255)
                                                                                             gccgGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  2:6147620/150‑1 (MQ=255)
                                                                                               cgGCTTGCGCAGCCATTGAGCCCTCGCCGTGGTGGGCCTGCTACAAAATGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGt                                                        <  2:2071729/150‑1 (MQ=255)
                                                                                                     gcgcAGCCTTTGCGCCCTCGCCGGGGGGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                                  <  2:5702903/150‑1 (MQ=255)
                                                                                                     gcgcAGCCATTGAGCCCTCGCGGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTc                                                  <  2:275979/150‑1 (MQ=255)
                                                                                                       gcAGCCTTTGCGCCCTCGGGGTGGTGCGCCTGCGAAAAAAAGGCGCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgc                                                <  1:1445657/150‑1 (MQ=255)
                                                                                                       gcAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgc                                                <  2:4959311/150‑1 (MQ=255)
                                                                                                          gCCATTGCGCCCTCGCGGGGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCgcctg                                             <  2:1461107/150‑1 (MQ=255)
                                                                                                              ttGCGCCCGCGCCGGGGTGCGCCTGCGACAAAACGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTc                                         <  2:7062805/150‑1 (MQ=255)
                                                                                                                 aGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGcc                                      <  2:4565334/150‑1 (MQ=255)
                                                                                                                     cgcGCCGTGGTGCGCCTGCGACAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                                  <  1:5334001/148‑1 (MQ=255)
                                                                                                                     cgcGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                                  <  2:1938871/148‑1 (MQ=255)
                                                                                                                     cTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                                  <  2:534899/150‑1 (MQ=255)
                                                                                                                     cTCGCCGTGGTGCGCCTACGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGc                                  <  2:4023187/150‑1 (MQ=255)
                                                                                                                         ccGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGtt                              <  2:6021485/150‑1 (MQ=255)
                                                                                                                             gtggcgcATGCAAAAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc                          <  1:314414/147‑1 (MQ=255)
                                                                                                                             ggtgCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc                          <  2:4393134/150‑1 (MQ=255)
                                                                                                                             ggtgCGCCTGCTACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTAGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAAc                          <  2:5081891/150‑1 (MQ=255)
                                                                                                                               tgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          >  2:7239932/1‑116 (MQ=255)
                                                                                                                               tgCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCttt                                                          <  1:7239932/116‑1 (MQ=255)
                                                                                                                                 cgcCGGCGACAAATAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTTCCTCGCCGTGCTGTTCAACctgc                      <  2:5490932/150‑1 (MQ=255)
                                                                                                                                  gcCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACctgct                     <  1:4283517/150‑1 (MQ=255)
                                                                                                                                    cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg                   >  1:3267752/1‑150 (MQ=255)
                                                                                                                                    cTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTgg                   >  2:6082118/1‑150 (MQ=255)
                                                                                                                                       ccaCAAAACGGCCCCCCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGa                <  1:2027838/148‑1 (MQ=255)
                                                                                                                                       cGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGa                >  1:3870368/1‑150 (MQ=255)
                                                                                                                                        gACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGAt               >  1:2518649/1‑150 (MQ=255)
                                                                                                                                        gACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGAt               >  2:5595837/1‑150 (MQ=255)
                                                                                                                                         aCAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATg              <  1:607078/150‑1 (MQ=255)
                                                                                                                                           aagatggCGCACCACGGAACTTGTTCCGCCCGGAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGGCCTTCTTTGTCTGTTCGCCGGGCTCGCCGTGCTGTTCCACCTGCTGGCGCTGcg            <  2:548513/145‑2 (MQ=255)
                                                                                                                                           aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCggg                                                                                            >  2:6674243/1‑70 (MQ=255)
                                                                                                                                           aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc            >  2:5716975/1‑150 (MQ=255)
                                                                                                                                           aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc            >  1:1085768/1‑150 (MQ=255)
                                                                                                                                           aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc            >  2:3611093/1‑150 (MQ=255)
                                                                                                                                           aaaaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGcc            >  1:1085802/1‑150 (MQ=255)
                                                                                                                                             aacggCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGc          <  1:6869649/147‑1 (MQ=255)
                                                                                                                                             aaaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGc          <  1:1196348/150‑1 (MQ=255)
                                                                                                                                              aaGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCt         >  2:2248714/1‑150 (MQ=255)
                                                                                                                                               cggCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt        <  2:5912432/149‑1 (MQ=255)
                                                                                                                                               aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt        <  1:4294034/150‑1 (MQ=255)
                                                                                                                                               aGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCtt        <  1:2382062/150‑1 (MQ=255)
                                                                                                                                                 gcgcACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATgcgc                                                                                                            >  1:6154433/1‑48 (MQ=255)
                                                                                                                                                  cgcACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGAc     <  1:5101708/150‑1 (MQ=255)
                                                                                                                                                   gcACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGAcc    <  2:4853140/150‑1 (MQ=255)
                                                                                                                                                     accacGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACCGt  <  2:6778332/150‑1 (MQ=255)
                                                                                                                                                     accacGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACCGt  <  2:1050373/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCAGCGCCGCCCCAGAGCCAAGGCTTCATGCTGGTAGGCACACTGTTCATCCTGCCCAGTTCGCTTGGCTGATCCGACTGTTCATGATGGCCGCCGGCATGCGCAGCCATTGAGCCCTCGCCGTGGTGCGCCTGCGACAAAAAGGCGCACCACGGAACTTGTACAGCCCGCAACGGATGCGGCTTAATGCGCGCCATATTCGCTCGGGTGACTCATGCTTTCACGTGCACGTGCACATCTTTGTCTGTTCGCCTGCCTCGCCGTGCTGTTCAACCTGCTGGCGATGCCGCTTGACCGT  >  NC_002947/3781238‑3781536

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: