Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 815,673 G→T 6.7% intergenic (+126/+66) uvrB → / ← ybhK excision nuclease subunit B/putative transferase YbhK

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913815,6730GT 6.7% 298.0 / 8.1 89intergenic (+126/+66)uvrB/ybhKexcision nuclease subunit B/putative transferase YbhK
Reads supporting (aligned to +/- strand):  ref base G (42/41);  new base T (2/4);  total (44/45)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.92e-01

GTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAG  >  NC_000913/815529‑815822
                                                                                                                                                |                                                                                                                                                     
gTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATcc                                                                                                                                                  <  2:286789/150‑1 (MQ=25)
  ttATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCgg                                                                                                                                                <  1:409053/150‑1 (MQ=255)
   tATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                                                                                                                               <  2:1355566/150‑1 (MQ=255)
    aTCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                              <  2:476511/150‑1 (MQ=255)
    aTCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCa                                                                                                                                              <  2:1430985/150‑1 (MQ=255)
      cgcgGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATa                                                                                                                                            <  1:770285/150‑1 (MQ=255)
             tCGTAACAGGATAGCGAAGACGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGGTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAcg                                                                                                                                     <  1:102238/150‑1 (MQ=255)
              cGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAAcgc                                                                                                                                    >  2:97650/1‑150 (MQ=255)
                    aGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAACCGCGCaaat                                                                                                                              <  1:116021/150‑4 (MQ=18)
                     ggATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCCTACAAACCGCGCaaatt                                                                                                                             <  2:786203/150‑5 (MQ=18)
                     ggATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACttt                                                                                                                             >  2:454025/1‑150 (MQ=255)
                      gATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTg                                                                                                                            <  2:428564/150‑1 (MQ=255)
                        tAGCGAAGAAGACTGATGACAAACGGAATACAGCCTGATGCTCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTc                                                                                                                          <  1:989526/150‑1 (MQ=255)
                        tAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTc                                                                                                                          <  1:378348/150‑1 (MQ=255)
                         aGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCt                                                                                                                         <  2:1314346/150‑1 (MQ=255)
                            gaagaaGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCa                                                                                                                      >  1:1259179/1‑150 (MQ=255)
                            gaagaaGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCa                                                                                                                      >  1:266126/1‑150 (MQ=255)
                             aagaagACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCaa                                                                                                                     <  1:1228516/150‑1 (MQ=255)
                               gaagaCTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATa                                                                                                                   <  1:1385530/150‑1 (MQ=255)
                                    cTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAg                                                                                                              <  2:591130/150‑1 (MQ=255)
                                     tgatgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGc                                                                                                             >  1:580683/1‑150 (MQ=255)
                                      gatgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGcc                                                                                                            <  1:740846/150‑1 (MQ=255)
                                       atgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGcc                                                                                                            <  2:17586/149‑1 (MQ=255)
                                       atgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGcc                                                                                                            >  1:24894/1‑149 (MQ=255)
                                       atgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGcc                                                                                                            <  2:24894/149‑1 (MQ=255)
                                       atgaCAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGcc                                                                                                            >  1:17586/1‑149 (MQ=255)
                                           cAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAt                                                                                                       >  1:35611/1‑150 (MQ=255)
                                             aaCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAttt                                                                                                     >  2:1565933/1‑150 (MQ=255)
                                              aCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtttt                                                                                                    >  1:1931998/1‑150 (MQ=255)
                                              aCGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAAtttt                                                                                                    <  1:1895467/150‑1 (MQ=255)
                                                 gAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCt                                                                                                 >  1:954872/1‑150 (MQ=255)
                                                   aaaCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTc                                                                                               >  1:1848395/1‑150 (MQ=255)
                                                    aaCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCg                                                                                              <  2:1332718/150‑1 (MQ=255)
                                                     aCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCgg                                                                                             <  1:642591/150‑1 (MQ=255)
                                                     aCAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCgg                                                                                             <  2:35611/150‑1 (MQ=255)
                                                        tccTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGAtt                                                                                          >  2:493857/2‑150 (MQ=255)
                                                            gATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCg                                                                                      <  2:14694/150‑1 (MQ=255)
                                                               gcgcTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTa                                                                                   <  2:1851864/150‑1 (MQ=255)
                                                                    aCGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTa                                                                              >  2:521716/1‑150 (MQ=255)
                                                                      gCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTaaa                                                                            >  1:1344637/1‑150 (MQ=255)
                                                                      gCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTaaa                                                                            >  1:1351447/1‑150 (MQ=255)
                                                                         tATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGcc                                                                         <  2:1997740/150‑1 (MQ=255)
                                                                         tATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGcc                                                                         <  1:1323758/150‑1 (MQ=255)
                                                                          aTCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCt                                                                        >  2:802886/1‑150 (MQ=255)
                                                                          aTCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCt                                                                        <  1:454025/150‑1 (MQ=255)
                                                                           tCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTg                                                                       <  2:1504182/150‑1 (MQ=255)
                                                                            cAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGt                                                                      >  2:69500/1‑150 (MQ=255)
                                                                                ccTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACg                                                                  >  2:1588876/1‑150 (MQ=255)
                                                                                  tACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGcc                                                                <  2:266126/150‑1 (MQ=255)
                                                                                     aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCttt                                                             <  1:192449/150‑1 (MQ=255)
                                                                                     aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCttt                                                             >  1:1984194/1‑150 (MQ=255)
                                                                                     aTTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCttt                                                             <  1:199396/150‑1 (MQ=255)
                                                                                      ttttCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtttt                                                            <  2:1778179/150‑1 (MQ=255)
                                                                                      ttttCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCtttt                                                            >  1:344172/1‑150 (MQ=255)
                                                                                       tttCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTc                                                           >  1:800770/1‑150 (MQ=255)
                                                                                          ctcCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAg                                                        >  1:1686186/1‑150 (MQ=255)
                                                                                             cGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAgcgc                                                     >  2:848533/1‑150 (MQ=255)
                                                                                                aaTATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGtt                                                  <  1:97650/150‑1 (MQ=255)
                                                                                                 atatatTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTa                                                 <  2:800770/150‑1 (MQ=255)
                                                                                                     atTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                             >  2:356322/1‑150 (MQ=255)
                                                                                                     atTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                             >  2:1743718/1‑150 (MQ=255)
                                                                                                     atTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                             <  1:56741/150‑1 (MQ=255)
                                                                                                     atTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTa                                             >  2:32903/1‑150 (MQ=255)
                                                                                                      ttGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCGTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaa                                            >  1:1066668/1‑150 (MQ=255)
                                                                                                       acaaagtGCGCGTTGTTTATGCCGGATGCGGCGTGAACGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaac                                           <  1:657428/144‑1 (MQ=25)
                                                                                                        gAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTATCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaaca                                          >  2:1301940/1‑150 (MQ=255)
                                                                                                         aaTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAAcc                                                                                >  1:711982/1‑111 (MQ=255)
                                                                                                         aaTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAAcc                                                                                <  2:711982/111‑1 (MQ=255)
                                                                                                         aaTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaacaa                                         >  1:1021854/1‑150 (MQ=255)
                                                                                                         aaTTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTaacaa                                         >  1:1307421/1‑150 (MQ=255)
                                                                                                               cgcgGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCg                                   >  1:506679/1‑150 (MQ=255)
                                                                                                                  ggTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTc                                >  1:640961/1‑150 (MQ=255)
                                                                                                                    tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGc                                    >  2:801439/1‑144 (MQ=255)
                                                                                                                    tttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGc                                    <  1:801439/144‑1 (MQ=255)
                                                                                                                     ttGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATg                             <  1:1588876/150‑1 (MQ=255)
                                                                                                                      tGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGa                            >  1:1585065/1‑150 (MQ=255)
                                                                                                                         aGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCAGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGa                         >  2:1531163/1‑150 (MQ=255)
                                                                                                                           gCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATa                       >  1:794889/1‑150 (MQ=255)
                                                                                                                            ccGGATGCGGCGTGAACGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATAc                      >  1:837205/1‑150 (MQ=255)
                                                                                                                               gATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGa                   <  2:1585065/150‑1 (MQ=255)
                                                                                                                               gATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTCTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGa                   <  1:1301940/150‑1 (MQ=255)
                                                                                                                                aTAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGaa                  <  1:1692636/150‑1 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAatat               >  1:685801/1‑150 (MQ=255)
                                                                                                                                   aGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAatat               >  2:769404/1‑150 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATc              <  1:493857/150‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATc              <  1:521716/150‑1 (MQ=255)
                                                                                                                                    ggCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATc              >  2:1303804/1‑150 (MQ=255)
                                                                                                                                         tCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTgg         >  1:1845121/1‑150 (MQ=255)
                                                                                                                                         gaaCGCCTTATCCGGCATAAACAACGCGCACTTTGTCTTCAATTGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTgg         >  1:1884753/3‑150 (MQ=255)
                                                                                                                                                gcATCCGTCATAAACAAAGCTCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTATAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGTCGGTCATGACGATACGTAATATCGCTGGGCTCCAg  <  2:437519/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                     
GTTTATCGCGGCATCGTAACAGGATAGCGAAGAAGACTGATGACAAACGGAAAACAGCCTGATGCGCTACGCTTATCAGGCCTACATTTTCTCCGCAATATATTGAATTTGCGCGGTTTGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCATAAACAACGCGCACTTTGTCTTCAATAGGAAGCCGGAATTTTCCTTCGGATTTCCGTTAACCTAAAGCCTGTAACGCCTTTTCCAGCGCGTTATGTAACAACTGGCGGTCATGACGATACGGAATATCGCTGGCCTCCAG  >  NC_000913/815529‑815822

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: