Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,760,260 T→G 7.1% F93V (TTC→GTC)  pvdM → dipeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,760,2600TG 7.1% 357.5 / 12.1 123F93V (TTC→GTC) pvdMdipeptidase
Reads supporting (aligned to +/- strand):  ref base T (62/51);  new base G (4/5);  total (67/56)
Fisher's exact test for biased strand distribution p-value = 7.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 126 total aligned reads displayed.

ACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGCTAC  >  NC_002947/4760114‑4760409
                                                                                                                                                  |                                                                                                                                                     
aCTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCg                                                                                                                                                    >  1:1343304/1‑150 (MQ=255)
  ttCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                                  >  2:1839519/1‑150 (MQ=255)
  ttCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCATTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGAGACGAGGGCCGCGAAGCGGACAAGGATGGCCCGGGCCAGATCGAc                                                                                                                                                  >  1:1121228/1‑150 (MQ=255)
    cGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCt                                                                                                                                                >  1:2180997/1‑150 (MQ=255)
     gACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTg                                                                                                                                               >  1:1442599/1‑150 (MQ=255)
     gACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTg                                                                                                                                               >  1:1765862/1‑150 (MQ=255)
      aCAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTgg                                                                                                                                              <  2:2227708/150‑1 (MQ=255)
       cAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGc                                                                                                                                             >  2:2294544/1‑150 (MQ=255)
           gACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCaaa                                                                                                                                         <  2:1765862/150‑1 (MQ=255)
            aCGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAg                                                                                                                                        <  1:2130441/150‑1 (MQ=255)
            aCGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGAACTGGcccggg                                                                                                                                        <  2:2023343/150‑5 (MQ=255)
                cTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgc                                                                                                                                    >  1:2061025/1‑150 (MQ=255)
                cTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgc                                                                                                                                    <  1:1387556/150‑1 (MQ=255)
                  aCCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgc                                                                                                                                  >  2:1839861/1‑150 (MQ=255)
                   ccccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcg                                                                                                                                 >  1:365641/1‑150 (MQ=255)
                   ccccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGCCACGGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcg                                                                                                                                 >  2:74768/1‑150 (MQ=255)
                    cccATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgcgcgc                                                                                                                                >  2:452760/1‑150 (MQ=255)
                       aTAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCggg                                                                                                                             <  2:632160/150‑1 (MQ=255)
                       aTAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGaa                                                                                                                                                  <  2:1044022/129‑2 (MQ=255)
                       aTAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGGTCGaa                                                                                                                                                  >  1:1044022/1‑128 (MQ=255)
                         agagGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCg                                                                                                                           >  1:864669/1‑150 (MQ=255)
                             gTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTg                                                                                                                       >  1:1283655/1‑150 (MQ=255)
                               cGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTc                                                                                                                     <  1:1528159/150‑1 (MQ=255)
                                     ccAGGCCAATGATCCGAAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGgcgc                                                                                                               >  2:1792656/1‑150 (MQ=255)
                                        ggCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCcgc                                                                                                                                    <  1:1660707/126‑1 (MQ=255)
                                         gCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGcc                                                                                                           >  2:733401/1‑150 (MQ=255)
                                         cccAATGATCTGCAAGCGCGGATCATTTCGTTCGACCGCCATCTCCCCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGcc                                                                                                           <  2:1331784/149‑1 (MQ=255)
                                             ataaTCTACAAGAGCGCATCATGCGGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                       <  2:1343304/147‑1 (MQ=255)
                                             aTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                       <  2:1185327/150‑1 (MQ=255)
                                             aTGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                       >  2:1564024/1‑150 (MQ=255)
                                                  cTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCATAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTc                                                                                                  <  1:2294544/150‑1 (MQ=255)
                                                    gCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCtt                                                                                                <  2:2061025/150‑1 (MQ=255)
                                                    gCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCGGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTAGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCtt                                                                                                >  2:2353616/1‑150 (MQ=255)
                                                     cAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGt                                                                                                                                                       <  1:1753140/94‑1 (MQ=255)
                                                     cAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGt                                                                                                                                                       >  2:1753140/1‑94 (MQ=255)
                                                     cAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTc                                                                                               >  2:378656/1‑150 (MQ=255)
                                                     cAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTc                                                                                               <  2:1442599/150‑1 (MQ=255)
                                                          gCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTg                                                                                          <  1:381544/150‑1 (MQ=255)
                                                              aTCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGccc                                                                                      <  1:1564024/150‑1 (MQ=255)
                                                                 aTTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAg                                                                                   <  1:841676/150‑1 (MQ=255)
                                                                  tttCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGAc                                                                                                                                                  >  2:1560312/1‑86 (MQ=255)
                                                                  tttCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGt                                                                                  >  2:1320057/1‑150 (MQ=255)
                                                                   ttcgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCGAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTc                                                                                 <  2:2302378/150‑1 (MQ=255)
                                                                   ttcgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTc                                                                                 >  2:828112/1‑150 (MQ=255)
                                                                     cgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCt                                                                               <  1:1839861/150‑1 (MQ=255)
                                                                     cgttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCt                                                                               <  1:1839519/150‑1 (MQ=255)
                                                                      gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCgg                                                                                                                              <  2:1878368/102‑1 (MQ=255)
                                                                      gttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGAATGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTg                                                                              >  1:344163/1‑150 (MQ=255)
                                                                       ttcgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTgg                                                                             >  1:664439/1‑150 (MQ=255)
                                                                         cgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTg                                                                                          <  1:1851203/135‑1 (MQ=255)
                                                                         cgACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTg                                                                                          >  2:1851203/1‑135 (MQ=255)
                                                                             aGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAg                                                                                   <  1:1829292/138‑1 (MQ=255)
                                                                             aGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGc                                                                       >  1:669950/1‑150 (MQ=255)
                                                                               ccATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                       <  2:118114/116‑1 (MQ=255)
                                                                               ccATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGa                                                                                                       >  1:118114/1‑116 (MQ=255)
                                                                               ccATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCCCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGa                                                                     <  1:739223/150‑1 (MQ=255)
                                                                                cATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGAc                                                                    >  1:2306631/1‑150 (MQ=255)
                                                                                cATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGAc                                                                    <  2:1522914/150‑1 (MQ=255)
                                                                                 aTCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCt                                                                               >  2:2096998/1‑138 (MQ=255)
                                                                                 aTCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACa                                                                   >  1:1538778/1‑150 (MQ=255)
                                                                                  tctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGc                                                                       >  2:1219546/1‑145 (MQ=255)
                                                                                  tctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACaa                                                                  <  1:1223894/150‑1 (MQ=255)
                                                                                   ctcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAAc                                                                 >  2:1939249/1‑150 (MQ=255)
                                                                                    tcACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACg                                                                <  2:365641/150‑1 (MQ=255)
                                                                                      aCCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                              >  1:914950/1‑150 (MQ=255)
                                                                                      aCCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                              >  1:915230/1‑150 (MQ=255)
                                                                                        ccTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGccc                                                            >  1:2307655/1‑150 (MQ=255)
                                                                                          tGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGc                                                          <  2:914950/150‑1 (MQ=255)
                                                                                               cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGcaccac                                                     <  2:664439/150‑1 (MQ=255)
                                                                                               cTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGcaccac                                                     <  1:733401/150‑1 (MQ=255)
                                                                                                   aCTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCAccacccca                                                 >  2:670585/1‑150 (MQ=255)
                                                                                                         gCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACg                                           >  1:1028065/1‑150 (MQ=255)
                                                                                                          cACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGCCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACgg                                          >  1:1255381/1‑150 (MQ=255)
                                                                                                           aCCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGt                                         >  2:609137/1‑150 (MQ=255)
                                                                                                             cGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGc                                                                       <  1:20539/118‑1 (MQ=255)
                                                                                                               aaGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcg                                     >  2:27108/1‑150 (MQ=255)
                                                                                                                aGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTtcgt                                    <  2:864669/150‑1 (MQ=255)
                                                                                                                    cgcgAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGAt                                <  2:1919621/150‑1 (MQ=255)
                                                                                                                       gAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGcc                             >  1:1943675/1‑150 (MQ=255)
                                                                                                                           cTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGc                         <  2:1112995/150‑1 (MQ=255)
                                                                                                                             gACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGa                       >  2:1662327/1‑150 (MQ=255)
                                                                                                                             gACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGCCGACCACGGCCCGCACCACCCCACCGCCGGTTTCGTCGATGCCGGGCgg                       >  2:2340387/1‑149 (MQ=255)
                                                                                                                              aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGAc                      >  2:1176122/1‑150 (MQ=255)
                                                                                                                              aCAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGCTCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGAc                      >  1:559838/1‑150 (MQ=255)
                                                                                                                                aaGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCg                                              <  1:2325706/124‑1 (MQ=255)
                                                                                                                                aaGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACAt                    >  2:1134707/1‑150 (MQ=255)
                                                                                                                                  ggATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                              <  1:2142290/106‑1 (MQ=255)
                                                                                                                                  ggATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGc                                                              >  2:2142290/1‑106 (MQ=255)
                                                                                                                                  ggATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGc                  >  1:970525/1‑150 (MQ=255)
                                                                                                                                   gATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAATCATGGACCGGCGACAACGACCGGCACCACCCCACCGACGGTTTCGGCGTTGCCGGGCGACATGcc                 >  2:793202/1‑150 (MQ=255)
                                                                                                                                     tGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCa               <  1:828112/150‑1 (MQ=255)
                                                                                                                                        ctttggCCAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGGCCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGa            >  1:1987173/5‑150 (MQ=255)
                                                                                                                                        cccGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGa            >  2:87025/1‑150 (MQ=255)
                                                                                                                                         tttggCCAGGTCGCACTGGCGAAAGCGGGGCGCGGGCGACTGTCCGGCGCTGCCCTGGCCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGaa           <  2:1987173/147‑1 (MQ=255)
                                                                                                                                         ccGGGCCAGTTCGACCTGGGCAGGGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGaa           <  2:559838/150‑1 (MQ=255)
                                                                                                                                         ccGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGaa           >  2:252661/1‑150 (MQ=255)
                                                                                                                                         ccGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGaa           <  1:609137/150‑1 (MQ=255)
                                                                                                                                              ccAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGAcc                                                                                                     >  2:2188838/1‑55 (MQ=255)
                                                                                                                                               cAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTgcgc     <  1:378656/150‑1 (MQ=255)
                                                                                                                                               cAGGTCGAACTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGcc                                                                                                           >  2:1492861/1‑48 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
ACTTCGACAAGGACGGCTACCCCATAGAGGTCGTGCGCCAGGCCAATGATCTGCAAGAGCGGATCATTTCGTTCGACAGCCATCTCACCCTGCCGCTGGACTTCGGCACCGAAGGCCGCGAAGCTGACAAGGATGGCCCGGGCCAGTTCGACCTGGCCAAAGCCGCGCGCGGGCGACTGTCCGGCGCTGCCCTGACCGTCTTCGCTTGGCCCGAGTCCTGGACCGGCGACAACGGCCCGCACCACCCCACCGACGGTTTCGTCGATGCCGGGCGACATGCCCAGGAAGTGCGCTAC  >  NC_002947/4760114‑4760409

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