Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,043,722 T→A 5.8% intergenic (+18/+16) gmk → / ← PP_5297 deoxyguanylate kinase/guanylate kinase/amino acid transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,043,7220TA 5.8% 320.6 / 8.4 103intergenic (+18/+16)gmk/PP_5297deoxyguanylate kinase/guanylate kinase/amino acid transporter
Reads supporting (aligned to +/- strand):  ref base T (43/54);  new base A (3/3);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 105 total aligned reads displayed.

CTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043575‑6043870
                                                                                                                                                    |                                                                                                                                                    
cTACGTCATCATCAATGATGATTTCGACGCGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTg                                                                                                                                                    >  1:1814424/1‑150 (MQ=255)
cTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTg                                                                                                                                                    >  1:536536/1‑150 (MQ=255)
 tACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGCTCGGCTGAAGATCGATGGGGCTGCATTGc                                                                                                                                                   <  1:1594430/150‑1 (MQ=255)
  aCGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCa                                                                                                                                                  <  1:1060832/150‑1 (MQ=255)
         catcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGACGGGGCTGCAATGCAGCCCCat                                                                                                                                           <  1:747556/150‑2 (MQ=255)
          atcaATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAAt                                                                                                                                          >  1:2207662/1‑150 (MQ=255)
            caATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATgg                                                                                                                                        >  2:547863/1‑150 (MQ=255)
            caATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTACGCTCGGATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATgg                                                                                                                                        >  2:2077984/1‑150 (MQ=255)
             aaTGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGc                                                                                                                                                   <  2:2247125/138‑1 (MQ=255)
             aaTGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGc                                                                                                                                                   >  1:2247125/1‑138 (MQ=255)
             aaTGATGATTTCGACACGGCGCTGGCAAACTTAAAAGCAGTGTTTTGCTGGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGt                                                                                                                                       <  2:1814424/150‑1 (MQ=255)
                gatgatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    <  1:1522168/150‑1 (MQ=255)
                gatgatTTCGACACGGCGCTGGCAGACTTCAAAGCAGTGTACCGCACAAATCAGCTGCTGCGGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    <  2:511199/150‑1 (MQ=255)
                 atgatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCGCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGCTTAg                                                                                                                                   >  2:41581/1‑150 (MQ=255)
                   gatTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTCCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTc                                                                                                                                 <  2:1344964/150‑1 (MQ=255)
                    atTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCGGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTcc                                                                                                                                >  1:709990/1‑150 (MQ=255)
                       tCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCGGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGc                                                                                                                             <  2:709990/150‑1 (MQ=255)
                           cacGGCGCTAGCAGACTTGAAAGCAGTGTTCCGCTCGACTCGGCTGCTGCTGAAGAAGCAGACCACCCGCGACGCGGCATTGCTCAAGCAAATGGTCTTCGGAAGATCGATTGGGCTTCATTGCAGCCCCAATGGTTTAGTCCTGCAGCg                                                                                                                         <  2:2017689/150‑1 (MQ=255)
                                cgcTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    <  2:1478899/150‑1 (MQ=255)
                                 gcTGGCAGACTTGAAAGCAGTGTTCCGTTCGAATCGGGTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGTTCAAGCAAATGATAGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAgg                                                                                                                   <  1:170391/150‑1 (MQ=255)
                                   tGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTa                                                                                                                 <  2:2076888/150‑1 (MQ=255)
                                     gCAGATTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGc                                                                                                               <  1:1647140/150‑1 (MQ=255)
                                     gCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCGTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGc                                                                                                               <  2:1770009/150‑1 (MQ=255)
                                     gCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGc                                                                                                               <  2:206033/150‑1 (MQ=255)
                                       aGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  2:1514858/1‑125 (MQ=255)
                                       aGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  1:1514858/125‑1 (MQ=255)
                                         aCTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGGTCGATGGTGCTGCATTGCGGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCctgc                                                                                                           <  2:266114/150‑1 (MQ=255)
                                          cTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCctgct                                                                                                          <  1:2207648/150‑1 (MQ=255)
                                             gAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    <  2:2220313/137‑1 (MQ=255)
                                             gAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    >  1:2220313/1‑137 (MQ=255)
                                                aGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCa                                                                                                    <  2:473676/150‑1 (MQ=255)
                                                  cAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  2:2099408/150‑1 (MQ=255)
                                                  cAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  2:2182878/150‑1 (MQ=255)
                                                  cAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  >  1:2182878/1‑150 (MQ=255)
                                                     tgtTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  2:1822838/1‑111 (MQ=255)
                                                          cGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGcccc                                                                                          >  1:2001045/1‑150 (MQ=255)
                                                           gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      >  2:668931/1‑105 (MQ=255)
                                                           gCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGtt                                                                                                                                      <  1:668931/105‑1 (MQ=255)
                                                            cTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        <  1:2077984/150‑1 (MQ=255)
                                                            cTCGAATCGGCGGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    <  1:631555/122‑1 (MQ=255)
                                                            cTCGAATCGGCGGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAg                                                                                                                    >  2:631555/1‑122 (MQ=255)
                                                               gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAggg                                                                                                 >  2:1876404/1‑138 (MQ=255)
                                                               gAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAggg                                                                                                 <  1:1876404/138‑1 (MQ=255)
                                                                 aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGAc                                                                                                       >  2:2326665/1‑130 (MQ=255)
                                                                 aTCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGAc                                                                                                       <  1:2326665/130‑1 (MQ=255)
                                                                  tCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGt                                                                                  >  2:538135/1‑150 (MQ=255)
                                                                     gctgctgctGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACg                                                                               >  1:160901/1‑150 (MQ=255)
                                                                          tgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTc                                                                                                                                 <  1:1395345/95‑1 (MQ=255)
                                                                          tgctgAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTc                                                                                                                                 >  2:1395345/1‑95 (MQ=255)
                                                                              gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGgc                                                                          >  1:1906927/1‑146 (MQ=255)
                                                                              gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGgc                                                                          <  2:1906927/146‑1 (MQ=255)
                                                                              gaagaaGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGcgcac                                                                      <  2:2001045/150‑1 (MQ=255)
                                                                                  aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAt                                                                                                                                                      >  2:241224/1‑66 (MQ=255)
                                                                                  aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAt                                                                                                                                                      <  1:241224/66‑1 (MQ=255)
                                                                                  aagCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCAt                                                                  <  2:536536/150‑1 (MQ=255)
                                                                                       gcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        <  2:1663822/123‑1 (MQ=255)
                                                                                       gcagcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGa                                                                                        >  1:1663822/1‑123 (MQ=255)
                                                                                         agcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGAc                                                           >  2:1498072/1‑150 (MQ=255)
                                                                                          gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCg                                                                                                                         <  1:695931/87‑1 (MQ=255)
                                                                                          gcagcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCg                                                                                                                         >  2:695931/1‑87 (MQ=255)
                                                                                              gcGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGca                                                       <  1:1501182/150‑1 (MQ=255)
                                                                                               cGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGcat                                                      <  1:2256289/150‑1 (MQ=255)
                                                                                                   aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  >  2:116143/1‑102 (MQ=255)
                                                                                                   aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAgg                                                                                                  <  1:116143/102‑1 (MQ=255)
                                                                                                   aTGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCatcatc                                                  >  1:2102366/1‑150 (MQ=255)
                                                                                                    tGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCAtcatca                                                 >  2:1443200/1‑150 (MQ=255)
                                                                                                      cGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAgg                                               >  2:1522538/1‑150 (MQ=255)
                                                                                                      cGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAgg                                               <  1:1434712/150‑1 (MQ=255)
                                                                                                       ggCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGc                                              <  1:547863/150‑1 (MQ=255)
                                                                                                        gCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTTCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTACCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGcc                                             >  1:2045758/1‑150 (MQ=255)
                                                                                                         cATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCg                                                                            >  2:2299940/1‑118 (MQ=255)
                                                                                                         cATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCg                                                                            <  1:2299940/118‑1 (MQ=255)
                                                                                                              cTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                           <  1:1348850/146‑1 (MQ=255)
                                                                                                              cTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGa                                           >  2:1348850/1‑146 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    <  1:623373/56‑1 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    >  2:623776/1‑56 (MQ=255)
                                                                                                               tCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    <  1:623776/56‑1 (MQ=255)
                                                                                                               tCAAGCAAATGAGCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTa                                                                                                                                    >  2:623373/1‑56 (MQ=255)
                                                                                                                 aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAg                                    >  1:1541697/1‑150 (MQ=255)
                                                                                                                 aaGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAg                                    >  1:1850296/1‑150 (MQ=255)
                                                                                                                      aaTGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCg                               >  2:1481362/1‑150 (MQ=255)
                                                                                                                        agagtggCTGAAGACCGATGTGTCTGCATTGCACACCCAACGGTTTAGTCCTGCAGCGGAAAGGTAGGCTGCTGTTGCAGGTTAGCCCCGAGGAAGTACGCCGGCGAACGCACGCGCGACAGCATCATCAGGCCGCTGCCAAGCACCgag                             <  2:849499/145‑1 (MQ=255)
                                                                                                                         gATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGAGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCgaga                            <  2:2207662/150‑1 (MQ=255)
                                                                                                                          aTCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           <  1:1481362/150‑1 (MQ=255)
                                                                                                                                tGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgc                     >  2:1269649/1‑150 (MQ=255)
                                                                                                                                    aaccattGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           >  2:1061266/7‑140 (MQ=255)
                                                                                                                                    aaccattGGGGCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           <  1:1061266/134‑1 (MQ=255)
                                                                                                                                     aTCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCaa                >  1:488552/1‑150 (MQ=255)
                                                                                                                                          tGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                          <  1:847174/103‑1 (MQ=255)
                                                                                                                                          tGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACa                                                          >  2:847174/1‑103 (MQ=255)
                                                                                                                                          tGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           >  2:2161773/1‑134 (MQ=255)
                                                                                                                                          tGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAgat                           <  1:2161773/134‑1 (MQ=255)
                                                                                                                                           ggggCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgc                     <  2:2276179/139‑1 (MQ=255)
                                                                                                                                           ggggCTGCAATGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATcgc                     >  1:2276179/1‑139 (MQ=255)
                                                                                                                                            gggCTGCATTGCAGCCCCAATGGTTTCGCCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGaa         <  2:2045758/150‑1 (MQ=255)
                                                                                                                                            gggCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGaa         <  1:1269649/150‑1 (MQ=255)
                                                                                                                                                tGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACAcc     >  1:1604751/1‑150 (MQ=255)
                                                                                                                                                  cATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       <  1:785086/66‑1 (MQ=255)
                                                                                                                                                  cATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAg                                                                                       >  2:785086/1‑66 (MQ=255)
                                                                                                                                                   aTTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACAccagc  <  1:1498072/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CTACGTCATCATCAATGATGATTTCGACACGGCGCTGGCAGACTTGAAAGCAGTGTTCCGCTCGAATCGGCTGCTGCTGAAGAAGCAGCAGCAGCGTCATGCGGCATTGCTCAAGCAAATGATCGGCTGAAGATCGATGGGGCTGCATTGCAGCCCCAATGGTTTAGTCCTGCAGCGGCAAGGTAGCCTGCTGACGCAGGGTAGCCCCGAGGAAGTACGCCGGCGCACGCATGCGCGACAGCATCATCAGGCCGATGCCAAGCACCGAGATGATCGCGGCAATCACGAACACCAGC  >  NC_002947/6043575‑6043870

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