Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,980 A→T 9.2% N65I (AAC→ATC)  PP_3531 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,9800AT 9.2% 346.9 / 19.6 118N65I (AAC→ATC) PP_3531hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (50/57);  new base T (6/5);  total (56/62)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.99e-01

Only 100 of 119 total aligned reads displayed.

TGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGAA  >  NC_002947/4002831‑4003125
                                                                                                                                                     |                                                                                                                                                 
tGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATaa                                                                                                                                                   >  1:706717/1‑150 (MQ=255)
  cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAcc                                                                                                                                                 >  1:2181251/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)
       gccgccGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCTTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCgg                                                                                                                                            <  1:1249563/150‑1 (MQ=255)
       gccgccGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCgg                                                                                                                                            <  1:288148/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)
                      cgcCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                             <  1:2110122/150‑1 (MQ=255)
                       gcCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGc                                                                                                                            >  2:1845371/1‑150 (MQ=255)
                        ccTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCt                                                                                                                           <  2:939915/150‑1 (MQ=255)
                          gggCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGt                                                                                                                         <  2:1321271/149‑1 (MQ=255)
                              cTACGACAGCGCAGGCATAAGCGTTTACATGAGCAGCGCGGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGACCTGCAAGTGGCTGTACGc                                                                                                                     >  2:1353787/1‑150 (MQ=255)
                              acgtaaacgcgAATGCCTTCGCGTTTACGTGAGCAGCGCAGGCGCCCCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTCCGTGTGGCCTGCAAGTGGCTGTACGc                                                                                                                     >  2:338814/8‑150 (MQ=255)
                               tACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGcc                                                                                                                    >  1:1557430/1‑150 (MQ=255)
                               tACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGcc                                                                                                                    >  1:2243709/1‑150 (MQ=255)
                                  gACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgt                                                                                                                 <  2:2243709/150‑1 (MQ=255)
                                   aCCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGGACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgtg                                                                                                                <  1:654209/150‑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)
                                                       tACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGAt                                                                                            <  2:249095/150‑1 (MQ=255)
                                                             aggagcGCAGGGGCCGCGTCCTACCAGCAGTTGCGGGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGGCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      <  2:104371/147‑1 (MQ=255)
                                                             agcagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      <  1:2264259/150‑1 (MQ=255)
                                                              gcagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     >  2:2263330/1‑150 (MQ=255)
                                                               cagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                                    <  1:2263330/150‑1 (MQ=255)
                                                                       gcgcCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTTACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGacaaca                                                                            >  2:330006/1‑150 (MQ=255)
                                                                        cgcCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAacac                                                                           <  2:788719/150‑1 (MQ=255)
                                                                           cTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTg                                                                        >  2:1932368/1‑150 (MQ=255)
                                                                            tCGTCCTACCAGCAGTTTGGGGGTTTCACCACCATGAAAGCCAGCGTGGGCACGCTCACTGACCTGCAGGATAACCAGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       <  1:338814/150‑1 (MQ=255)
                                                                            tCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       <  1:32696/150‑1 (MQ=255)
                                                                              gTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGgt                                                                     >  2:2099181/1‑150 (MQ=255)
                                                                               tCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCaa                                                                                                                                 <  2:1788504/89‑1 (MQ=255)
                                                                               tCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCaa                                                                                                                                 >  1:1788504/1‑89 (MQ=255)
                                                                                  tACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTAc                                                                 <  1:827562/150‑1 (MQ=255)
                                                                                     cagcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGACTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTg                                                              <  1:847730/150‑1 (MQ=255)
                                                                                      agcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTg                                                                                                                                    >  2:2363216/1‑79 (MQ=255)
                                                                                      agcagTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCGCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGa                                                             >  1:255389/1‑150 (MQ=255)
                                                                                                     tCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCAt                                              >  2:286735/1‑150 (MQ=255)
                                                                                                     tCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCAt                                              >  1:1503481/1‑150 (MQ=255)
                                                                                                       aGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      >  1:95198/1‑108 (MQ=255)
                                                                                                       aGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      <  2:95198/108‑1 (MQ=255)
                                                                                                        gCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           <  1:2185812/150‑1 (MQ=255)
                                                                                                           ccATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCgg                                        <  1:385584/150‑1 (MQ=255)
                                                                                                            cATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGcctcg                                       >  1:2012485/1‑147 (MQ=255)
                                                                                                            cATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGc                                       >  2:1024810/1‑150 (MQ=255)
                                                                                                             aTGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                                    >  2:2097886/1‑104 (MQ=255)
                                                                                                              tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     <  1:478339/102‑1 (MQ=255)
                                                                                                              tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     >  2:478339/1‑102 (MQ=255)
                                                                                                              tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                     >  1:801315/1‑150 (MQ=255)
                                                                                                               gAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                     >  2:336191/1‑149 (MQ=255)
                                                                                                                 aaGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGc                                  <  1:1230856/150‑1 (MQ=255)
                                                                                                                  aGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGtca                                                                                                                                        >  2:1738695/1‑45 (MQ=9)
                                                                                                                  aGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGtca                                                                                                                                        <  1:1738695/47‑3 (MQ=9)
                                                                                                                            cgcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAACTGGCTGTACGCCTGTCACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGTCCGCGGGACATc                       >  1:753167/1‑150 (MQ=255)
                                                                                                                             gcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGa                                                                                 <  1:745539/91‑1 (MQ=255)
                                                                                                                             gcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGa                                                                                 >  2:745539/1‑91 (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                                                                              <  2:744032/86‑1 (MQ=255)
                                                                                                                                     aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGacaa                                                                              >  1:744032/1‑86 (MQ=255)
                                                                                                                                     aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCgg                             >  1:464961/1‑135 (MQ=255)
                                                                                                                                     aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa              >  2:962224/1‑150 (MQ=255)
                                                                                                                                     aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa              <  1:2031390/150‑1 (MQ=255)
                                                                                                                                        cacccgCAGGTTATCCTGCAGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           <  1:2333412/113‑1 (MQ=255)
                                                                                                                                        caCCCGCAGGTTATCCTGCAGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           >  2:2333412/2‑118 (MQ=255)
                                                                                                                                         acccgCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCt                  <  2:913335/138‑1 (MQ=255)
                                                                                                                                         aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg                    <  2:257193/140‑1 (MQ=255)
                                                                                                                                         aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg                    >  1:257193/1‑140 (MQ=255)
                                                                                                                                         aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGa          <  2:1212440/150‑1 (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)
                                                                                                                                          cccgCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       <  1:2130421/85‑1 (MQ=255)
                                                                                                                                          cccgCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       >  2:2130421/4‑88 (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:1986909/150‑1 (MQ=255)
                                                                                                                                            tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa       >  1:2091610/1‑150 (MQ=255)
                                                                                                                                                 gATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGaa  >  1:2013436/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGAA  >  NC_002947/4002831‑4003125

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: ‑