Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,002,977 A→T 9.2% D64V (GAT→GTT)  PP_3531 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,002,9770AT 9.2% 362.9 / 22.3 120D64V (GAT→GTT) PP_3531hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (51/58);  new base T (6/5);  total (57/63)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.47e-01

Only 100 of 120 total aligned reads displayed.

TCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGAA  >  NC_002947/4002829‑4003125
                                                                                                                                                    |                                                                                                                                                    
tctGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGAt                                                                                                                                                     <  1:1567658/150‑1 (MQ=255)
    cccTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAAcc                                                                                                                                                 >  2:1016483/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)
         gccgccGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCgg                                                                                                                                            <  1:2208212/150‑1 (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                                                                                                                             <  1:2110122/150‑1 (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)
                                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:2367109/1‑150 (MQ=255)
                                     aCCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGGACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCtgtg                                                                                                                <  1:654209/150‑1 (MQ=255)
                                             ggCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATg                                                                                                        >  2:596001/1‑150 (MQ=255)
                                              gCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGc                                                                                                       >  2:2340915/1‑150 (MQ=255)
                                                aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                             <  1:11600/126‑1 (MQ=255)
                                                aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTgg                                                                                                                             >  2:11600/1‑126 (MQ=255)
                                                aTTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCgg                                                                                                     <  1:1680086/150‑1 (MQ=255)
                                                 ttCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGgc                                                                                                    >  2:1512115/1‑150 (MQ=255)
                                                            gTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGAtgtg                                                                                         >  1:788719/1‑150 (MQ=255)
                                                                gcagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     >  2:2263330/1‑150 (MQ=255)
                                                                 cagcGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                                    <  1:2263330/150‑1 (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:32696/150‑1 (MQ=255)
                                                                              tCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       <  1:596001/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)
                                                                                       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                                              >  2:286735/1‑150 (MQ=255)
                                                                                                       tCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCAt                                              >  1:1503481/1‑150 (MQ=255)
                                                                                                         aGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      <  2:95198/108‑1 (MQ=255)
                                                                                                         aGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAg                                                                                      >  1:95198/1‑108 (MQ=255)
                                                                                                          gCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           <  1:2185812/150‑1 (MQ=255)
                                                                                                          gCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGATGACCAAACCTCGTGGTACCGGACCACCCAGCTGCCATGGc                                           >  2:1230856/1‑150 (MQ=255)
                                                                                                             ccATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCgg                                        <  1:385584/150‑1 (MQ=255)
                                                                                                              cATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGcctcg                                       >  1:2012485/1‑147 (MQ=255)
                                                                                                              cATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGc                                       >  2:1024810/1‑150 (MQ=255)
                                                                                                               aTGAAAGCCAGCGTGCGCACGCTCACTGCCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAggg                                                                                    <  1:2097886/104‑1 (MQ=255)
                                                                                                                tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     >  2:478339/1‑102 (MQ=255)
                                                                                                                tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAgg                                                                                     <  1:478339/102‑1 (MQ=255)
                                                                                                                tGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGa                                     >  1:801315/1‑150 (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)
                                                                                                                    aGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGtca                                                                                                                                        >  2:1738695/1‑45 (MQ=9)
                                                                                                                    aGCCAGCGTGCGCACGCTCACTGACCCGCAGGTTATCCTGCAGGtca                                                                                                                                        <  1:1738695/47‑3 (MQ=9)
                                                                                                                              cgcACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATc                       <  2:2181251/150‑1 (MQ=255)
                                                                                                                              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)
                                                                                                                                     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)
                                                                                                                                       aCTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGaaa              >  2:962224/1‑150 (MQ=255)
                                                                                                                                          cacccgCAGGTTATCCTGCAGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           <  1:2333412/113‑1 (MQ=255)
                                                                                                                                          caCCCGCAGGTTATCCTGCAGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGc                                           >  2:2333412/2‑118 (MQ=255)
                                                                                                                                           acccgCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCt                  <  2:913335/138‑1 (MQ=255)
                                                                                                                                           aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg                    >  1:257193/1‑140 (MQ=255)
                                                                                                                                           aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTg                    <  2:257193/140‑1 (MQ=255)
                                                                                                                                           aCCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGa          <  2:1212440/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                                                                       >  2:2130421/4‑88 (MQ=255)
                                                                                                                                            cccgCAGGTTATCCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTgg                                                                       <  1:2130421/85‑1 (MQ=255)
                                                                                                                                            ccTGCAGGATAACCTGCGGGTGGCCTGCAAACGGCAGTTCAGACGTGTGCTGATGCGGCTGCTGGATGTGGAGACTTACTGCACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGAc         <  1:2340915/150‑1 (MQ=255)
                                                                                                                                              tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATgg                                            >  1:1431947/1‑113 (MQ=255)
                                                                                                                                              tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATgg                                            <  2:1431947/113‑1 (MQ=255)
                                                                                                                                              tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa       <  1:1986909/150‑1 (MQ=255)
                                                                                                                                              tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa       >  1:2091610/1‑150 (MQ=255)
                                                                                                                                              tGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCa       >  1:1276638/1‑150 (MQ=255)
                                                                                                                                                   gATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGaa  >  1:2013436/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCTGCCCTGGCCGCCGATGACTGGCGCCTGGCCTACGACCGCGAAGGCATTCGCGTTTACGTGAGCAGCGCAGGCGCCTCGTCCTACCAGCAGTTTCGTGGTGTCAGCACCATGAAAGCCAGCGTGCGCACGCTCACTGACCTGCAGGATAACCTGCGGGTGGCCTGCAAGTGGCTGTACGCCTGTGACCAGATGCGGCTGCTGGATGTGGAGGGTGACAACACCTGGGTGTACCTGACCACCCAGCTGCCATGGCCGGCGAGGCCGCGGGACATCGTGCTGAAAGTGACCAGCGAA  >  NC_002947/4002829‑4003125

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: