Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,971 T→C 6.7% L62P (CTG→CCG)  PP_3531 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,9710TC 6.7% 347.0 / 8.3 115L62P (CTG→CCG) PP_3531hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (49/57);  new base C (5/3);  total (54/61)
Fisher's exact test for biased strand distribution p-value = 4.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 121 total aligned reads displayed.

CACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCA  >  NC_002947/4002825‑4003120
                                                                                                                                                  |                                                                                                                                                     
cACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCa                                                                                                                                                    <  1:1384654/150‑1 (MQ=255)
cACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCa                                                                                                                                                    <  1:925084/150‑1 (MQ=255)
 aCCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAg                                                                                                                                                   >  1:911998/1‑150 (MQ=255)
    tctGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGAt                                                                                                                                                <  1:1567658/150‑1 (MQ=255)
      tGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATaa                                                                                                                                              >  1:706717/1‑150 (MQ=255)
        cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAcc                                                                                                                                            >  1:2181251/1‑150 (MQ=255)
        cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAcc                                                                                                                                            >  2:1016483/1‑150 (MQ=255)
          cTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCGGCAGGATAACCTg                                                                                                                                          <  1:2060324/150‑1 (MQ=255)
           tGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGc                                                                                                                                         <  1:1615412/150‑1 (MQ=255)
            ggCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCg                                                                                                                                        <  2:1259058/150‑1 (MQ=255)
             gccgccGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCgg                                                                                                                                       <  1:2208212/150‑1 (MQ=255)
               cgccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGt                                                                                                                                     >  1:2141600/1‑150 (MQ=255)
                gccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTg                                                                                                                                    <  2:1205730/150‑1 (MQ=255)
                 ccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTgg                                                                                                                                   <  2:706717/150‑1 (MQ=255)
                 ccgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGtca                                                                                                                                   <  1:1008797/150‑3 (MQ=255)
                  cgATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGc                                                                                                                                  <  2:911998/150‑1 (MQ=255)
                   gATGACTGGCGCCTGGCCTACGACCGAGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGcc                                                                                                                                 >  1:571386/1‑150 (MQ=255)
                       aCTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCa                                                                                                                             <  1:316617/150‑1 (MQ=255)
                         tGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAg                                                                                                                           >  2:2110122/1‑150 (MQ=255)
                            cgcCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                        <  2:782174/150‑1 (MQ=255)
                             gcCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGc                                                                                                                       >  2:1845371/1‑150 (MQ=255)
                              ccTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCt                                                                                                                      >  2:1986909/1‑150 (MQ=255)
                              ccTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCt                                                                                                                      <  2:939915/150‑1 (MQ=255)
                                gggCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGt                                                                                                                    <  2:1321271/149‑1 (MQ=255)
                                    acgtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCCCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTCCGTGTGGCCTGCAAGTGGCTGTACGc                                                                                                                >  2:338814/8‑150 (MQ=255)
                                     tACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGcc                                                                                                               >  1:2243709/1‑150 (MQ=255)
                                        gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgt                                                                                                            <  2:2243709/150‑1 (MQ=255)
                                        gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgt                                                                                                            >  2:2367109/1‑150 (MQ=255)
                                              gAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCa                                                                                                                                                    >  1:1328573/1‑104 (MQ=255)
                                              gAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCa                                                                                                                                                    <  2:1328573/104‑1 (MQ=255)
                                                 ggCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATg                                                                                                   >  2:596001/1‑150 (MQ=255)
                                                  gCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGc                                                                                                  >  2:2340915/1‑150 (MQ=255)
                                                    aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                        >  2:11600/1‑126 (MQ=255)
                                                    aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                        <  1:11600/126‑1 (MQ=255)
                                                    aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCgg                                                                                                <  1:1680086/150‑1 (MQ=255)
                                                     ttCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGgc                                                                                               >  2:1512115/1‑150 (MQ=255)
                                                                gTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGAtgtg                                                                                    >  1:788719/1‑150 (MQ=255)
                                                                   aggagcGCAGGGGCCGCGTCCTACCAGCAGTTGCGGGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGGCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                 <  2:104371/147‑1 (MQ=255)
                                                                   agcagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                 <  1:2264259/150‑1 (MQ=255)
                                                                    gcagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                >  2:2263330/1‑150 (MQ=255)
                                                                       gcgcAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTg                                                                             <  1:1657354/150‑1 (MQ=255)
                                                                             gcgcCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTTACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGacaaca                                                                       >  2:330006/1‑150 (MQ=255)
                                                                              cgcCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAacac                                                                      <  2:788719/150‑1 (MQ=255)
                                                                                ccTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACAccg                                                                    <  2:1557430/150‑2 (MQ=255)
                                                                                 cTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTg                                                                   >  2:1932368/1‑150 (MQ=255)
                                                                                  tCGTCCTACCAGCAGTTTGGGGGTTTCACCACCATGAAAGCCAGCGTGGGCACGCTCACTGACCTGCAGGATAACCAGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                  <  1:338814/150‑1 (MQ=255)
                                                                                  tCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                  <  1:596001/150‑1 (MQ=255)
                                                                                  tCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                  <  1:32696/150‑1 (MQ=255)
                                                                                     tCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCaa                                                                                                                            >  1:1788504/1‑89 (MQ=255)
                                                                                     tCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCaa                                                                                                                            <  2:1788504/89‑1 (MQ=255)
                                                                                        tACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTAc                                                            <  1:827562/150‑1 (MQ=255)
                                                                                           cagcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGACTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTg                                                         <  1:847730/150‑1 (MQ=255)
                                                                                            agcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTg                                                                                                                               >  2:2363216/1‑79 (MQ=255)
                                                                                            agcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTg                                                                                                                               <  1:2363216/79‑1 (MQ=255)
                                                                                            agcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCGCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGa                                                        >  1:255389/1‑150 (MQ=255)
                                                                                                           tCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCAt                                         >  1:1503481/1‑150 (MQ=255)
                                                                                                           tCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCAt                                         >  2:286735/1‑150 (MQ=255)
                                                                                                             aGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                 <  2:95198/108‑1 (MQ=255)
                                                                                                              gCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                      <  1:2185812/150‑1 (MQ=255)
                                                                                                              gCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGATGACCAAACCTCGTGGTACCGGACCACCCAGCTGCCATGGc                                      >  2:1230856/1‑150 (MQ=255)
                                                                                                                 ccATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCgg                                   <  1:385584/150‑1 (MQ=255)
                                                                                                                  cATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGc                                  >  2:1024810/1‑150 (MQ=255)
                                                                                                                   aTGAAAGCCAGCGTGCGCACGCTCACTGCCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                               <  1:2097886/104‑1 (MQ=255)
                                                                                                                   aTGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                               >  2:2097886/1‑104 (MQ=255)
                                                                                                                    tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                >  2:478339/1‑102 (MQ=255)
                                                                                                                    tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                <  1:478339/102‑1 (MQ=255)
                                                                                                                    tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                <  1:1928100/150‑1 (MQ=255)
                                                                                                                     gAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                >  2:336191/1‑149 (MQ=255)
                                                                                                                     gAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                <  1:336191/149‑1 (MQ=255)
                                                                                                                       aaGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGc                             <  1:1230856/150‑1 (MQ=255)
                                                                                                                                  cgcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATc                  <  2:2181251/150‑1 (MQ=255)
                                                                                                                                  cgcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAACTGGCTGTACGCCTGTCACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGTCCGCGGGACATc                  >  1:753167/1‑150 (MQ=255)
                                                                                                                                   gcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGa                                                                            >  2:745539/1‑91 (MQ=255)
                                                                                                                                   gcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGa                                                                            <  1:745539/91‑1 (MQ=255)
                                                                                                                                     aCGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg               <  1:1016483/150‑1 (MQ=255)
                                                                                                                                     aCGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg               <  2:2141600/150‑1 (MQ=255)
                                                                                                                                       gCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTa                                                                                                                   <  2:582508/48‑1 (MQ=255)
                                                                                                                                       gCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTa                                                                                                                   >  1:582508/1‑48 (MQ=255)
                                                                                                                                       gCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCGGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCt             >  1:1602024/1‑150 (MQ=255)
                                                                                                                                         tCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGa           >  1:1212440/1‑150 (MQ=255)
                                                                                                                                           ccttaCCCGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCGGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa         <  1:1845371/146‑1 (MQ=255)
                                                                                                                                           aCTGACCTGCAGGATAACCTGCGGGTGGCCTGTAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCgg                        <  2:464961/135‑1 (MQ=255)
                                                                                                                                           aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGacaa                                                                         >  1:744032/1‑86 (MQ=255)
                                                                                                                                           aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGacaa                                                                         <  2:744032/86‑1 (MQ=255)
                                                                                                                                           aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa         >  2:962224/1‑150 (MQ=255)
                                                                                                                                           aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa         <  1:2031390/150‑1 (MQ=255)
                                                                                                                                              caCCCGCAGGTTATCCTGCAGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                      >  2:2333412/2‑118 (MQ=255)
                                                                                                                                               aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg               <  2:257193/140‑1 (MQ=255)
                                                                                                                                               aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg               >  1:257193/1‑140 (MQ=255)
                                                                                                                                               aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGa     <  1:1792954/150‑1 (MQ=255)
                                                                                                                                               aCCCGCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCt             >  1:913335/1‑142 (MQ=255)
                                                                                                                                               aCCCGCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGa     >  2:541148/1‑150 (MQ=255)
                                                                                                                                               aCCCGCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGa     >  2:1383547/1‑150 (MQ=255)
                                                                                                                                                ccTGCAGGATAACCTGCGGGTGGCCTGCAAACGGCAGTTCAGACGTGTGCTGATGCGGCTGCTGGATGTGGAGACTTACTGCACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGAc    <  1:2340915/150‑1 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATgg                                       <  2:1431947/113‑1 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATgg                                       >  1:1431947/1‑113 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGGGAGGCCGGGGGACATCGTGCTGCAAGTGAcaa  >  2:1599964/1‑148 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa  >  1:1276638/1‑150 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa  <  1:1986909/150‑1 (MQ=255)
                                                                                                                                                  tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa  >  1:2091610/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CACCTCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCA  >  NC_002947/4002825‑4003120

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