Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,517,557 T→A 6.0% E54V (GAA→GTA) ‡ fecR ← regulator for fec operon, periplasmic

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,517,5570TA 6.0% 586.2 / 18.0 182E54V (GAA→GTA) ‡fecRregulator for fec operon, periplasmic
Reads supporting (aligned to +/- strand):  ref base T (84/87);  new base A (4/7);  total (88/94)
Fisher's exact test for biased strand distribution p-value = 5.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.75e-01

Only 100 of 191 total aligned reads displayed.

GTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAAC  >  NC_000913/4517410‑4517706
                                                                                                                                                   |                                                                                                                                                     
 tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCaa                                                                                                                                                    <  1:1219696/150‑1 (MQ=255)
 tttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGCGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAACCTGGTTGCGCAGGTTTTCaa                                                                                                                                                    >  1:968742/1‑150 (MQ=255)
  ttCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAAc                                                                                                                                                   <  2:1280125/150‑1 (MQ=255)
     cGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTGATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGGCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctg                                                                                                                                                <  2:1087790/150‑1 (MQ=255)
     cGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctg                                                                                                                                                <  2:2263769/150‑1 (MQ=255)
     cGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCAAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGAGAGGCGGGTATCGTGCAACGCCCGGCTGGCACCGTCGCCAGCCACCCCGCCAAGCTGGTTGCGAAGGGTTTCAACctg                                                                                                                                                >  1:170611/1‑150 (MQ=255)
       aCTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgct                                                                                                                                              <  1:259430/150‑1 (MQ=255)
       aCTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACctgct                                                                                                                                              <  2:1611525/150‑1 (MQ=255)
       aCTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcg                                                                                                                                              <  2:1911912/150‑2 (MQ=255)
        cTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgc                                                                                                                                             <  2:133006/150‑3 (MQ=255)
          gCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGcc                                                                                                                                           >  2:173515/1‑150 (MQ=255)
          gCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCAACCCGCCACGCTGGTTGCGCAGGTTTTCAACCTGCTGcc                                                                                                                                           >  2:902936/1‑150 (MQ=255)
           ccAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgcaac                                                                                                                                          <  2:2419410/150‑6 (MQ=255)
           ccAGAGCTGCAAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTGAAAACctgcgcaac                                                                                                                                          <  1:2437151/150‑6 (MQ=255)
            cAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAg                                                                                                                                         >  1:978933/1‑150 (MQ=255)
             agagCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAgg                                                                                                                                        >  2:266487/1‑150 (MQ=255)
              gagCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGc                                                                                                                                       <  2:323983/150‑1 (MQ=255)
                  tGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                                   >  1:1754606/1‑150 (MQ=255)
                  tGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                                   >  2:844596/1‑150 (MQ=255)
                  tGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCACCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCAc                                                                                                                                   >  2:170092/1‑150 (MQ=255)
                   gCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACt                                                                                                                                  >  2:474345/1‑150 (MQ=255)
                      aGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGt                                                                                                                               >  2:1830843/1‑150 (MQ=255)
                      aGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGt                                                                                                                               <  2:1117381/150‑1 (MQ=255)
                         cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTat                                                                                                                            >  2:1662443/1‑150 (MQ=255)
                         cTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTat                                                                                                                            >  1:1963838/1‑150 (MQ=255)
                            cGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcct                                                                                                                         >  1:1568460/1‑150 (MQ=255)
                            cGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcct                                                                                                                         >  1:950789/1‑150 (MQ=255)
                             gCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctg                                                                                                                        >  2:310637/1‑150 (MQ=255)
                             gCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctg                                                                                                                        >  1:1525947/1‑150 (MQ=255)
                              ccAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctga                                                                                                                       >  2:947623/1‑150 (MQ=255)
                                aGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTatcctgatc                                                                                                                     <  1:1005992/150‑1 (MQ=255)
                                    ccGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGt                                                                                                                 >  1:2196670/1‑150 (MQ=255)
                                    ccGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGt                                                                                                                 >  1:936802/1‑150 (MQ=255)
                                      gAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTc                                                                                                               >  1:2272243/1‑150 (MQ=255)
                                      gAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTc                                                                                                               <  1:2382963/150‑1 (MQ=255)
                                       aGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCa                                                                                                              <  1:2210738/150‑1 (MQ=255)
                                       aGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCa                                                                                                              <  1:878139/150‑1 (MQ=255)
                                        gCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCAt                                                                                                             >  2:47408/1‑150 (MQ=255)
                                         cAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATa                                                                                                            <  2:1640086/150‑1 (MQ=255)
                                          aaCAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATAc                                                                                                           <  1:1873294/150‑1 (MQ=255)
                                              gcagTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACt                                                                                                       >  2:934347/1‑150 (MQ=255)
                                                  tAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTg                                                                                                   <  2:1845634/150‑1 (MQ=255)
                                                      ccTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAg                                                                                               >  2:571358/1‑150 (MQ=255)
                                                       cTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgc                                                                                              >  1:246481/1‑150 (MQ=255)
                                                          tcatcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgcgcg                                                                                           >  2:1777111/1‑150 (MQ=255)
                                                           catcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAgcgcgc                                                                                          >  1:1367278/1‑150 (MQ=255)
                                                            atcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCCGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCt                                                                                         <  2:2108464/150‑1 (MQ=255)
                                                             tcaCGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGGCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCtt                                                                                        <  1:1474255/150‑1 (MQ=255)
                                                                cGTGTCGGCGGGGGTTGGGGGTTTCGTTCAACGGCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCt                                                                                     <  1:2047980/150‑1 (MQ=255)
                                                                      ggcggGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTtgtg                                                                               >  1:1694826/1‑150 (MQ=255)
                                                                       gcggGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTgg                                                                              <  1:295230/150‑1 (MQ=255)
                                                                         gggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGc                                                                            <  2:1872800/150‑1 (MQ=255)
                                                                          ggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCt                                                                           >  2:69329/1‑150 (MQ=255)
                                                                          ggTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCt                                                                           >  2:668090/1‑150 (MQ=255)
                                                                            tGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGa                                                                         >  2:1170397/1‑150 (MQ=255)
                                                                                 cGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACgcgc                                                                    <  2:1787630/150‑1 (MQ=255)
                                                                                 cGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACgcgc                                                                    <  1:1462917/150‑1 (MQ=255)
                                                                                  gggTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCt                                                                   >  1:1009362/1‑150 (MQ=255)
                                                                                  gggTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCt                                                                   >  2:2054017/1‑150 (MQ=255)
                                                                                    gTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcg                                                                 >  2:769817/1‑150 (MQ=255)
                                                                                     tATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcgc                                                                >  1:756267/1‑150 (MQ=255)
                                                                                      aTCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTcgcc                                                               <  1:266487/150‑1 (MQ=255)
                                                                                         gTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCt                                                            >  1:272915/1‑150 (MQ=255)
                                                                                         gTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCt                                                            >  1:809060/1‑150 (MQ=255)
                                                                                            cAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAg                                                         >  1:539909/1‑150 (MQ=255)
                                                                                            cAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAg                                                         <  1:2124005/150‑1 (MQ=255)
                                                                                             aaCGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGc                                                        <  2:1754606/150‑1 (MQ=255)
                                                                                                gCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACg                                                     <  2:1128311/150‑1 (MQ=255)
                                                                                                 cccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACgg                                                    <  1:2454302/150‑1 (MQ=255)
                                                                                                  ccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGc                                                   <  1:844596/150‑1 (MQ=255)
                                                                                                  ccGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGc                                                   <  2:1458402/150‑1 (MQ=255)
                                                                                                       tGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATAcc                                              <  2:238855/150‑1 (MQ=255)
                                                                                                          cAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                           >  2:1509552/1‑150 (MQ=255)
                                                                                                          cAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGt                                           >  1:131409/1‑150 (MQ=255)
                                                                                                            aCGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTgg                                         <  1:851707/150‑1 (MQ=255)
                                                                                                               tCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGaa                                      >  2:2030681/1‑150 (MQ=255)
                                                                                                                 gCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGc                                    <  1:668090/150‑1 (MQ=255)
                                                                                                                      gCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                               >  2:407014/1‑150 (MQ=255)
                                                                                                                      gCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAAc                               <  1:524427/150‑1 (MQ=255)
                                                                                                                       cacaCCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACg                              >  1:1904715/1‑150 (MQ=255)
                                                                                                                            cGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAgcg                         >  2:2455532/1‑150 (MQ=255)
                                                                                                                                aaGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTg                     <  2:246481/150‑1 (MQ=255)
                                                                                                                                  gCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAc                   <  1:583684/150‑1 (MQ=255)
                                                                                                                                       ttGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggcgg              <  1:1647185/150‑1 (MQ=255)
                                                                                                                                       ttGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGAcggcgg              <  2:1188276/150‑1 (MQ=255)
                                                                                                                                        tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg             <  1:351533/150‑1 (MQ=255)
                                                                                                                                        tGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCggg             <  1:1956294/150‑1 (MQ=255)
                                                                                                                                          cgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGaa           >  2:2010850/1‑150 (MQ=255)
                                                                                                                                          cgcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCGGCACGCAGCGCCTGACGGCGGGaa           >  1:1050496/1‑150 (MQ=255)
                                                                                                                                           gcAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAAt          >  1:411218/1‑150 (MQ=255)
                                                                                                                                            cAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATc         >  2:807075/1‑150 (MQ=255)
                                                                                                                                                ttttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGtt     >  1:683332/1‑150 (MQ=255)
                                                                                                                                                ttttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGtt     >  2:537631/1‑150 (MQ=255)
                                                                                                                                                ttttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGtt     <  1:1620712/150‑1 (MQ=255)
                                                                                                                                                 tttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATc         >  1:1585221/1‑145 (MQ=255)
                                                                                                                                                 tttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATc         <  2:1585221/145‑1 (MQ=255)
                                                                                                                                                 tttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTa    >  2:381210/1‑150 (MQ=255)
                                                                                                                                                  ttCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTaa   <  2:1807857/150‑1 (MQ=255)
                                                                                                                                                   tCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAAc  <  2:2312971/150‑1 (MQ=255)
                                                                                                                                                   tCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAAc  <  1:474345/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTTTCCGACTGCCAGAGCTGCCAGCCTCCGCCAGCGCCGAGCAACAGCAGTAATCCTTTCATCACGTGACGGCGGGTGAGGCGGGTATCGTGCAACGCCCGGCTGGCAACGTCGCCAGGCACACCGCCAAGCTGGTTGCGCAGGTTTTCAACCTGCTGCCAGGCCCACTGGTTATCCTGATCCTGTTCATACCACTGTTGCCAGCGCGCTTCCTGTTGTGGGCTGACGCGCTCGCCGCTTAGCACGGCATACCAGTGGGAAGCTGAACGCAGCGCCTGACGGCGGGAATCGGTTAAC  >  NC_000913/4517410‑4517706

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