Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 830,388 C→T 7.2% A86T (GCG→ACG)  ybiH ← DUF1956 domain‑containing tetR family putative transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913830,3880CT 7.2% 275.6 / 8.5 83A86T (GCG→ACG) ybiHDUF1956 domain‑containing tetR family putative transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (48/29);  new base T (3/3);  total (51/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CAGGTGGTAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGGGCGGCTATC  >  NC_000913/830239‑830536
                                                                                                                                                     |                                                                                                                                                    
cAGGTGGTAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTgcgc                                                                                                                                                      >  2:1591261/1‑150 (MQ=255)
cAGGTGGTAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTgcgc                                                                                                                                                      <  1:1573973/150‑1 (MQ=255)
    tggtAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAAc                                                                                                                                                  <  2:1435138/150‑1 (MQ=255)
    tggtAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAAc                                                                                                                                                  >  2:1094624/1‑150 (MQ=255)
      gtAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACaa                                                                                                                                                >  2:862998/1‑150 (MQ=255)
           cTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGtt                                                                                                                                           <  1:956723/150‑1 (MQ=255)
            tGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTc                                                                                                                                          <  2:805582/150‑1 (MQ=255)
                 tGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCtt                                                                                                                                     >  1:1210748/1‑150 (MQ=255)
                 tGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCtt                                                                                                                                     >  1:1210589/1‑150 (MQ=255)
                 tGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCtt                                                                                                                                     >  1:1156885/1‑150 (MQ=255)
                 tGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCtt                                                                                                                                     <  2:660040/150‑1 (MQ=255)
                   ggAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTcc                                                                                                                                   >  2:32593/1‑150 (MQ=255)
                           ctgctCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATg                                                                                                                           <  2:1283469/150‑1 (MQ=255)
                           ctgctCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATg                                                                                                                           >  1:2994/1‑150 (MQ=255)
                             gctCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcg                                                                                                                         >  2:577749/1‑150 (MQ=255)
                             gctCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcg                                                                                                                         <  1:1650661/150‑1 (MQ=255)
                               tCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcgga                                                                                                                       <  1:507075/150‑1 (MQ=255)
                               tCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGcgga                                                                                                                       <  2:20310/150‑1 (MQ=255)
                                    ggAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGaa                                                                                                                  >  2:408222/1‑150 (MQ=255)
                                            aCTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGc                                                                                                          <  1:748588/150‑1 (MQ=255)
                                              ttGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCa                                                                                                        >  1:1594709/1‑150 (MQ=255)
                                              ttGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCa                                                                                                        >  1:620909/1‑150 (MQ=255)
                                                  tGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATaa                                                                                                    >  2:592715/1‑150 (MQ=255)
                                                      gTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATc                                                                                                >  2:305583/1‑150 (MQ=255)
                                                       ttGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCg                                                                                               >  2:731697/1‑150 (MQ=255)
                                                        tGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCgg                                                                                              <  1:1672218/150‑1 (MQ=255)
                                                         gACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGc                                                                                             <  1:773301/150‑1 (MQ=255)
                                                          aCGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCa                                                                                            >  2:1457763/1‑150 (MQ=255)
                                                            ggTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAAt                                                                                          >  1:1154762/1‑150 (MQ=255)
                                                                          tGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAgg                                                                            <  2:1415681/150‑1 (MQ=255)
                                                                             gcagcTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGa                                                                         >  1:1662635/1‑150 (MQ=255)
                                                                             gcagcTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCACGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGa                                                                         <  2:68885/150‑1 (MQ=255)
                                                                               agcTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAgg                                                                       >  1:504508/1‑150 (MQ=255)
                                                                                gcTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGt                                                                      >  2:815769/1‑150 (MQ=255)
                                                                                 cTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTa                                                                     <  1:975193/150‑1 (MQ=255)
                                                                                  ttAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTAc                                                                    <  1:1506873/150‑1 (MQ=255)
                                                                                       cATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtc                                                               >  1:1304214/1‑150 (MQ=255)
                                                                                       cATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtc                                                               >  1:1304223/1‑150 (MQ=255)
                                                                                       cATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtc                                                               >  1:1304220/1‑150 (MQ=255)
                                                                                        aTGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtct                                                              >  1:1641742/1‑150 (MQ=255)
                                                                                        aTGTTCCAGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtct                                                              >  1:1546242/1‑150 (MQ=255)
                                                                                         tGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtctt                                                             <  2:65429/150‑1 (MQ=255)
                                                                                          gTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtcttc                                                            >  2:1581629/1‑150 (MQ=255)
                                                                                          gTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAAtcttc                                                            >  2:755483/1‑150 (MQ=255)
                                                                                            tCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATcttctt                                                          <  2:141685/150‑1 (MQ=255)
                                                                                             ccTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCttt                                                         <  2:1130926/150‑1 (MQ=255)
                                                                                             ccTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCttt                                                         >  2:428329/1‑150 (MQ=255)
                                                                                              cTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTc                                                        >  1:313973/1‑150 (MQ=255)
                                                                                              cTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTc                                                        >  1:313978/1‑150 (MQ=255)
                                                                                                 cAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGaa                                                     <  1:471372/150‑1 (MQ=255)
                                                                                                  aGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAAc                                                    <  2:636352/150‑1 (MQ=255)
                                                                                                    gcgcGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCg                                                  >  2:523830/1‑150 (MQ=255)
                                                                                                    gcgcGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCg                                                  >  2:411562/1‑150 (MQ=255)
                                                                                                     cgcgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                                                 <  2:1156885/150‑1 (MQ=255)
                                                                                                     cgcgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                                                 <  2:1210748/150‑1 (MQ=255)
                                                                                                     cgcgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGTGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                                                 <  2:1210589/150‑1 (MQ=255)
                                                                                                     cgcgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGa                                                 <  1:1652205/150‑1 (MQ=255)
                                                                                                       cgAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAag                                               >  2:799013/1‑150 (MQ=255)
                                                                                                        gAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGTAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAagt                                              >  2:1293872/1‑150 (MQ=255)
                                                                                                        gAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAagt                                              >  2:83570/1‑150 (MQ=255)
                                                                                                            gATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTCCAAATCTTCTTTCGAACCGAAgtagta                                          >  1:158998/1‑150 (MQ=255)
                                                                                                            gATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAgtagta                                          >  1:1339036/1‑150 (MQ=255)
                                                                                                            gATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAgtagta                                          >  1:926909/1‑150 (MQ=255)
                                                                                                             aTCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGtagtag                                         >  1:677665/1‑150 (MQ=255)
                                                                                                             aTCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGtagtag                                         >  1:605931/1‑150 (MQ=255)
                                                                                                                aGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTg                                      >  1:1021061/1‑150 (MQ=255)
                                                                                                                 gTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCACTCCACTGGGCGCAGGCGAGGTACACATCTTCTTTCGAACCGAAGGAGTAGGTGa                                     >  1:54307/1‑150 (MQ=255)
                                                                                                                 gTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGa                                     <  2:1203403/150‑1 (MQ=255)
                                                                                                                  ttCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGAt                                    >  2:825349/1‑150 (MQ=255)
                                                                                                                    cACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATgg                                  >  1:1284501/1‑150 (MQ=255)
                                                                                                                     aCGGAGGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTGTCGAACCGACGGAGTAGGTGATGGc                                 >  2:1724431/1‑150 (MQ=255)
                                                                                                                           ggCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGat                           >  1:909101/1‑150 (MQ=255)
                                                                                                                            gCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGata                          >  1:616200/1‑150 (MQ=255)
                                                                                                                             cAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGatat                         <  1:1489957/150‑1 (MQ=255)
                                                                                                                             cAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGatat                         >  2:1637754/1‑150 (MQ=255)
                                                                                                                                 ccGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTg                     <  1:930718/150‑1 (MQ=255)
                                                                                                                                    aTCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGTAGTAGTAGGTGATgg                                  >  1:1062040/1‑134 (MQ=255)
                                                                                                                                    aTCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATgg                                  <  2:1062040/134‑1 (MQ=255)
                                                                                                                                       aGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGc               <  2:1572662/150‑1 (MQ=255)
                                                                                                                                          cTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTg            >  1:249665/1‑150 (MQ=255)
                                                                                                                                           tgtgGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGTCGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTgg           <  2:1638518/150‑1 (MQ=255)
                                                                                                                                               gTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGggcgg       <  1:1185425/150‑1 (MQ=255)
                                                                                                                                               gTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGggcgg       <  2:1521088/150‑1 (MQ=255)
                                                                                                                                               gTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGggcga       >  2:1176955/1‑149 (MQ=255)
                                                                                                                                                ttGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGggcggc      >  2:541023/1‑150 (MQ=255)
                                                                                                                                                    gcgAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGGGCGGCTAtc  >  2:1709253/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAGGTGGTAGGCTGCCGTGGGAGAGAGCTGCTCACGGGAGATAAACTTGCTGAGGTTGACGGTGTCATCCTGGGTGAGCAGCTTAATCATGTTCCTGCAGGCGCGAAGGATCAGTTCACGGATGGCAGCCCGATCAGGCTGTGGTTGCGCGAACAAGCGTTCGGCTTCCTCGGCATGCGGACGGAACTGCTCGCCAATAAAATCGGCAATCCACTGGGCGCAGGCGAGGTACAAATCTTCTTTCGAACCGAAGTAGTAGGTGATGGCAGCGATATTCTGCCCGGCCTGGGCGGCTATC  >  NC_000913/830239‑830536

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