Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 868,055 A→C 6.1% intergenic (+150/+27) PP_0750 → / ← mqo‑I hypothetical protein/malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947868,0550AC 6.1% 966.7 / 29.5 303intergenic (+150/+27)PP_0750/mqo‑Ihypothetical protein/malate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (153/129);  new base C (11/10);  total (164/139)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 303 total aligned reads displayed.

TAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAGA  >  NC_002947/867910‑868203
                                                                                                                                                    |                                                                                                                                                    
tAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCCCACCCCCACGAGTCGCGCCCCACCGGCCGGCCATACCGAAACAGTGGGTGCGGGCAATCCCCCGACGTGGGCGGTACAGGCG‑‑‑GAACAAGCCCCACGCACAAAAAAAACCCCCAgggg                                                                                                                                                  >  2:3406895/1‑150 (MQ=255)
 aGCCTGTGCGGGCCTCTTCGCGGGCGGGCCCGCCCCCACAGGAGTCGCGCCGCACCTGCGGGCTTTACCGACCCTTTGGGAGGGGGCATGCCCGCGAAGTGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGc                                                                                                                                                 <  1:6773443/150‑1 (MQ=255)
     aGTGCCGGCCTCGTCGCGGGCGCGCCCGCTCCCACAGGAGTCCCGCCCCACCGGCCGGCCTTACCGCCCCGGGGGGAGCGGGCAAGCCCGCCCAGAGGGCGGTCCAGGCG‑‑‑GATCAAGCGCCCGGCCCAAAAAACGCCCCGCGGGGCggtg                                                                                                                                             >  2:7306849/1‑147 (MQ=255)
       tGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                                                                           <  2:5839898/150‑1 (MQ=255)
       tGCCGGCCTCTACGCGGGCGCCCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGAATTACCGAACCTGTGGGTGCGTGAATGCCCTCTAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGCGttttt                                                                                                                                           <  1:604125/150‑1 (MQ=255)
       ggCCGGCCCCTTCCGGGGCGCCCCCGCCCCCACCGGAGCCGCGCCGCACCTGGGGGCATTACCGAACGTGTGGGCGCGGGAATGCCGGCGAAGTGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttttt                                                                                                                                           <  1:3354407/149‑1 (MQ=255)
        gCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCTCGGGGCGtttttt                                                                                                                                          <  2:1026133/150‑2 (MQ=255)
          cGGCCTCTTCGCGGGCGCGCCCGCTCCCACTGGAGTCGCGCCGCACCTGCGGGCTTTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCTCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCg                                                                                                                                        <  2:3604462/150‑1 (MQ=255)
           ggCCTGTGCGGGGGCGCGCCCGCTCCCCCAGGTGTCGCGCCGCACCTGGGGGAATTACCGACCCTGTGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGt                                                                                                                                       <  1:198647/150‑1 (MQ=255)
           ggCCTCTTCGCGGGCGCGCCCGCTCCCACCGGAGTCGCGCCGCACCTGCGGGCATTACCGACCCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGt                                                                                                                                       <  2:6541538/150‑1 (MQ=255)
           ggCCTCGTCGCGGGCGCGCCCGCTCCCCCCGGAGTCGCGCCGCACCGGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGt                                                                                                                                       <  1:5265057/150‑1 (MQ=255)
             ccTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGttt                                                                                                                                     <  1:4303215/150‑1 (MQ=255)
                 ttcgcgggcatgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactccggtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGa                                                                                                                                 >  2:3179383‑M2/72‑150 (MQ=255)
                  tCGCGGGCGCGCCCGCCCCCCCCGGCGTCGCGCCGCACCTGCGGGCATTCCCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGAc                                                                                                                                <  1:4114372/150‑1 (MQ=255)
                     cgGGCGCGCCCGCTCCCACAGGAGTAGCCCCCCACCTGCGGGCCTTACCGAACCGGTGGGCGCGGGCAATCCCCCGCAGAGGGCGGTACAGGCA‑‑‑GATCAAGCGCCAGGCACAAAAAAAGCCCCGAGGGGCGGTTTTGCGGTTGAGACCgg                                                                                                                             >  2:3987792/1‑150 (MQ=255)
                     cgGGCGCGCCCCCTCCCACAGGAGTCGCGCCCCACCATCGGGCCATACCGCACCCGTGGGAGCGGGCATGCCCCCGCAGAGGGCGGTACAGGCA‑‑‑GATCAAGCGCCAGGCACAAAAAAACCCCCCCGGGGCGGTTTTGCGGTTCAGAACgg                                                                                                                             >  2:1591530/1‑150 (MQ=255)
                       ggCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATGACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCt                                                                                                                           <  2:2092090/150‑1 (MQ=255)
                        gcgcgcCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGACGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCtt                                                                                                                          >  2:4898669/1‑150 (MQ=255)
                          atgcccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTaa                                                                                                                        >  1:1583096‑M2/63‑150 (MQ=255)
                           cgcCCGCTCCCCCCGGCGTCGCGCCGCACCGGGGGGAATTACCGACCGTGGGGGTGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAg                                                                                                                       <  1:2965446/150‑1 (MQ=255)
                            GCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGc                                                                                                                      <  2:6639832‑M1/150‑1 (MQ=255)
                              ccgctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGt                                                                                                                    >  1:2868843‑M2/59‑150 (MQ=255)
                               CGCTCCCCCCGGAGGCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTa                                                                                                                   >  2:7149209‑M1/1‑150 (MQ=255)
                                gctcccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAg                                                                                                                  >  2:6033927‑M2/57‑150 (MQ=255)
                                GCTCCCCCAGGAGCCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAg                                                                                                                  >  1:1825254‑M1/1‑150 (MQ=255)
                                 ctcccacaggttcggtaatgcccgcaggtgcggcgcaactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGt                                                                                                                 >  2:554840‑M2/56‑150 (MQ=255)
                                  TCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTg                                                                                                                <  1:6915842‑M1/150‑1 (MQ=255)
                                  TCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTg                                                                                                                <  1:6346447‑M1/150‑1 (MQ=255)
                                    ccacaggttcggtaatgcccgtaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGCCGGCACTGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                              >  1:3018563‑M2/53‑150 (MQ=255)
                                    ccacaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                              >  1:2643653‑M2/53‑150 (MQ=255)
                                    CCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                              <  2:1826316‑M1/150‑1 (MQ=255)
                                    CCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                              <  2:47271‑M1/150‑1 (MQ=255)
                                     aACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGGGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                             <  2:6375745‑M1/149‑1 (MQ=255)
                                     CACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAATCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                             <  2:3241507‑M1/150‑1 (MQ=255)
                                     CACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGctg                                                                                                             <  1:6380487‑M1/150‑2 (MQ=255)
                                     CACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCtt                                                                                                             <  1:5429123‑M1/150‑1 (MQ=255)
                                      accggAGTCGCGCCGCACCTGCGGGCATTCCCGACCCGGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                            <  1:2314613/147‑1 (MQ=255)
                                      accggAGGCGCGCCGCACCTGCGGGCATTACCGAACCGGTGGGCGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                            <  1:5750160/147‑1 (MQ=255)
                                      acaggttcggtaatgcccgcaggtgcggcgcgactcctgtgggagcgggcGCGCCCGCGAAGAGGCCGGCACTGGCTACAGATCAGCCACGGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCt                                                                                                              >  1:3143586‑M2/51‑150 (MQ=255)
                                      ACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAATCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTg                                                                                                            <  1:508560‑M1/150‑1 (MQ=255)
                                       CAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCGTGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGttt                                                                                                                                             >  2:7395279‑M1/1‑116 (MQ=255)
                                                    gcACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGgcttct                                                                                              <  2:2272150/150‑1 (MQ=255)
                                                        cggCGGGCATTACCGAACCTGTGGGAGCGGGCCTGCCCGCGACGCGGGCGGTACAGGCT‑‑‑GATCACGCGCCAGGCACACAACACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                                          <  1:4742518/148‑1 (MQ=255)
                                                        cTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGctctct                                                                                          <  1:7355399/150‑1 (MQ=255)
                                                            gggCATTACCGACCCTGGGGGCGCGGGCATGCCCGCGACGAGGGCGGTACCGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCa                                                                                      <  1:1777114/150‑1 (MQ=255)
                                                            gggCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCa                                                                                      >  2:5453404/1‑150 (MQ=255)
                                                                  tACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCa                                                                                >  2:2669854/1‑150 (MQ=255)
                                                                   aCCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAg                                                                               >  2:6289036/1‑150 (MQ=255)
                                                                    ccGAACCTGTGGGAGCGGGCATGCCCGCTAAGAGGGCGGTAAAGGCT‑‑‑GATAAAGCGCCAGGAACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                              <  1:4026930/150‑1 (MQ=255)
                                                                    ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                              >  2:6394973/1‑150 (MQ=255)
                                                                    ccGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAgg                                                                              >  2:6070032/1‑150 (MQ=255)
                                                                            gtgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtc                                                                      <  1:5966360‑M2/138‑1 (MQ=255)
                                                                             tgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc                                                                                                                            <  2:2367600‑M2/84‑1 (MQ=255)
                                                                             tgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGc                                                                                                                            >  1:2367600‑M2/12‑95 (MQ=255)
                                                                             tgggagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGcgtcg                                                                     >  2:726368‑M2/12‑150 (MQ=255)
                                                                                gagcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTa                                                                  >  2:827436‑M2/9‑150 (MQ=255)
                                                                                  gcgggcATGCCCGCGAAGCGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATc                                                                <  1:6217415‑M2/144‑1 (MQ=255)
                                                                                  gcgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATc                                                                >  1:495751‑M2/7‑150 (MQ=255)
                                                                                   cgggcATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCg                                                               <  1:4918984‑M2/145‑1 (MQ=255)
                                                                                           cccGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGAt                                                       >  2:5051164/1‑150 (MQ=255)
                                                                                            ccGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATc                                                      >  2:4434425/1‑150 (MQ=255)
                                                                                            ccGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATc                                                      <  2:2868843/150‑1 (MQ=255)
                                                                                              gcgAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCtt                                                    >  2:7037505/1‑150 (MQ=255)
                                                                                              gcgAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCtt                                                    >  1:3457675/1‑150 (MQ=255)
                                                                                              gcgAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCtt                                                    >  1:3937487/1‑150 (MQ=255)
                                                                                                gacgcgggCGGTACAGGCT‑‑‑GATCAAGCGCCCGGCACACAACACGCCCCGCGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgt                                                  <  1:4485120/145‑1 (MQ=255)
                                                                                                gAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgt                                                  <  1:1090393/150‑1 (MQ=255)
                                                                                                 aaGAGGGCGGTCCCGGCT‑‑‑GAGCCCGCGCCAGGCACAACCCACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                 <  1:321873/150‑1 (MQ=255)
                                                                                                 aaGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                 >  2:4105167/1‑150 (MQ=255)
                                                                                                 aaGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                 <  2:6989320/150‑1 (MQ=255)
                                                                                                 aaGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtgtg                                                 <  1:7202842/150‑1 (MQ=255)
                                                                                                    agGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTa                                              >  1:1148796/1‑150 (MQ=255)
                                                                                                      ggCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACa                                            >  1:6625646/1‑150 (MQ=255)
                                                                                                              aGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTtg                                                   >  1:844045/1‑135 (MQ=255)
                                                                                                                gCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGTGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATc                                  >  1:5883139/1‑150 (MQ=255)
                                                                                                                gCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGTCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATc                                  >  2:3085178/1‑150 (MQ=255)
                                                                                                                 cT‑‑‑GATCAAGCGACAGGCACAAAAAACGCCACGAGGGGCGTGTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCGTGTGGTACAGCGCAGCAGCg                                 >  2:2240438/1‑150 (MQ=255)
                                                                                                                  t‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCgg                                >  2:6857274/1‑150 (MQ=255)
                                                                                                                      gATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGt                               <  1:4826147/150‑1 (MQ=255)
                                                                                                                       aTCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTAt                                                                 >  1:3160141/1‑115 (MQ=255)
                                                                                                                       aTCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGtt                              >  1:1434528/1‑150 (MQ=255)
                                                                                                                       aTCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGtt                              <  1:7430513/150‑1 (MQ=255)
                                                                                                                          aaGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGcc                           >  1:5550055/1‑150 (MQ=255)
                                                                                                                             cgcCAGGCACAAAAAACGCCCCGCGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGTAGCATCGGTTGCCAgg                        <  1:2448896/150‑1 (MQ=255)
                                                                                                                             ccacgggCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAgg                        <  2:380061/145‑1 (MQ=255)
                                                                                                                               ccAGGCACAAAAAACGACCCGAGGGGCGTTTTAGCGTATCCGACAGGCTTAACCGTAGTGATTTGCCGCCGGGCTGGTCTAAGCCAGGTCCCGGGCCGCGTCGTTATCGGCAATGATCGTGTGGGACCGCCCAGCAGCGGTTGCCAGGGt                      >  1:723949/1‑150 (MQ=255)
                                                                                                                                cAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGtc                     >  1:2123064/1‑150 (MQ=255)
                                                                                                                                    cacaAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                       >  2:2793454/1‑128 (MQ=255)
                                                                                                                                    cacaAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCgca                                       <  1:2793454/128‑1 (MQ=255)
                                                                                                                                     acaAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCg                                 <  2:1111639/133‑1 (MQ=255)
                                                                                                                                     acaAAAAACGCCCCGAGGGGCGTGTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCCCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTtctc                >  1:6019134/1‑150 (MQ=255)
                                                                                                                                      caAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcg               >  1:928847/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcgc              >  1:561668/1‑150 (MQ=255)
                                                                                                                                       aaaaaaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTcgc              >  1:3575865/1‑150 (MQ=255)
                                                                                                                                         aaaaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGa            >  1:3711673/1‑150 (MQ=255)
                                                                                                                                         aaaaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGa            <  2:1460314/150‑1 (MQ=255)
                                                                                                                                           aaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGc          >  1:548475/1‑150 (MQ=255)
                                                                                                                                           aaCGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGc          >  2:133465/1‑150 (MQ=255)
                                                                                                                                               ccccGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCggg      <  2:52161/150‑1 (MQ=255)
                                                                                                                                                 ccGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGaa    <  1:3800081/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TAGCCAGTGCCGGCCTCTTCGCGGGCGCGCCCGCTCCCACAGGAGTCGCGCCGCACCTGCGGGCATTACCGAACCTGTGGGAGCGGGCATGCCCGCGAAGAGGGCGGTACAGGCT‑‑‑GATCAAGCGCCAGGCACAAAAAACGCCCCGAGGGGCGTTTTTGCGTTTCAGACCGGCTTAAGCGTAGTGCTTGGCCGCCGGGCTGCTCTCTGCCAGGTCCAGGGCCGCGTCGTTATCGGCACTGATCTTGTGGTACAGCGCAGCATCGGTTGCCAGGGTCTTCTCGCGAGCCGGGAAGA  >  NC_002947/867910‑868203

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: