Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,653,262 Δ1 bp 7.5% coding (654/1020 nt) rspB ← putative Zn‑dependent NAD(P)‑binding oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,653,2620A. 7.5% 160.7 / 6.0 53coding (654/1020 nt)rspBputative Zn‑dependent NAD(P)‑binding oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (23/25);  new base . (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

GCAACATGCTGGAAATCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTCCACGCATTTTTTGCCGGAACCACGGCATA  >  NC_000913/1652966‑1653558
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
gCAACATGCTGGAAATTAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCg                                                                                                                                                                                                                                                                                                      <  1:236429/301‑1 (MQ=255)
 tctcttccgatctaTCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCgg                                                                                                                                                                                                                                                                                                     <  2:72590/288‑1 (MQ=255)
         tGGAAATCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGc                                                                                                                                                                                                                                                                                             <  2:238038/301‑1 (MQ=255)
            aaaTCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGt                                                                                                                                                                                                                                                                                          <  1:152299/301‑1 (MQ=255)
              aTCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCggagatcggaaaagg                                                                                                                                                                                                                                                                                        >  1:72500/1‑288 (MQ=255)
               tCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGttat                                                                                                                                                                                                                                                                                       >  2:579313/1‑301 (MQ=255)
               aaaaaaGTAGGGGAAATTAATGTTGTTGGTTTAAGTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTGAGTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAAGTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGttat                                                                                                                                                                                                                                                                                       <  1:325061/299‑1 (MQ=255)
                 aaaCGTTTGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGttatta                                                                                                                                                                                                                                                                                     <  2:257974/301‑1 (MQ=255)
                       ccgatctattaattTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGcc                                                                                                                                                                                                                                                                               <  2:232311/294‑1 (MQ=255)
                              attaattTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGccagaacgg                                                                                                                                                                                                                                                                        >  1:232068/1‑294 (MQ=255)
                                           ttttAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCttt                                                                                                                                                                                                                                                           <  2:202680/300‑1 (MQ=255)
                                            tttAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTc                                                                                                                                                                                                                                                          <  1:578597/301‑1 (MQ=255)
                                                      ccTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCa                                                                                                                                                                                                                                                >  2:600934/1‑301 (MQ=255)
                                                                         tAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATc                                                                                                                                                                                                                             >  2:315018/1‑301 (MQ=255)
                                                                            cGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTTTCTGGCTGTTATTAATCTCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCATTCGTTCATCAATTCGATCGGc                                                                                                                                                                                                                          >  2:427667/1‑301 (MQ=255)
                                                                             gATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGca                                                                                                                                                                                                                         <  1:314614/301‑1 (MQ=255)
                                                                                cGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAAtacaaacgccacca                                                                                                                                                                                                                      >  1:87358/1‑288 (MQ=255)
                                                                                            tGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTa                                                                                                                                                                                                          >  2:556112/1‑301 (MQ=255)
                                                                                               ttttAAGCGTGAAGAGAAAATATAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCACTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACg                                                                                                                                                                                                       <  2:489464/300‑1 (MQ=255)
                                                                                                      cGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGa                                                                                                                                                                                                >  1:281260/1‑301 (MQ=255)
                                                                                                        tGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCCCTCTTTCGCTTTTTCCAGTTGTTCAACAATGCGAACGGACACAATCACATTTTTAACGTTATAGAcc                                                                                                                                                                                              >  1:233141/1‑300 (MQ=255)
                                                                                                              aaaaaTAGTGATTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCtttt                                                                                                                                                                                        <  1:555398/300‑1 (MQ=255)
                                                                                                                 aaTAGTGAGTTCTTTTCAGTTAATTCCTTGCGTAATCACTTCAGACTGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAAt                                                                                                                                                                                     <  1:427103/301‑1 (MQ=255)
                                                                                                                     gagagTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAGCGGTGTCTGGCTGTTATTAATCGCCAAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCtt                                                                                                                                                                                >  1:180605/1‑300 (MQ=255)
                                                                                                                              tttCCGGTTATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACGATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTGCAGGATAGAAGGATGACAAGCCGCATCGATAAGTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCg                                                                                                                                                                        <  1:143006/301‑1 (MQ=255)
                                                                                                                              tttCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCg                                                                                                                                                                        <  2:281563/301‑1 (MQ=255)
                                                                                                                                    gTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCAGCCCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATAGGCAAAAATAACAATTTTAACGTTATAGACGACTTTTAATAACTAAACGAACGTCAAGcc                                                                                                                                                                 >  2:221291/1‑301 (MQ=255)
                                                                                                                                                                       ggAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCCCATTTTTAACGTTATAGACTCCTTTTAATACCTGAAAGATCGTCAGGCAGATTGGACCGGCACCATAAACCAGAACggaaaca                                                                                                                               >  1:67466/1‑296 (MQ=255)
                                                                                                                                                                                          atacGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATg                                                                                                           <  1:600173/301‑1 (MQ=255)
                                                                                                                                                                                                  cGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGATCGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGtt                                                                                                    >  1:11476/1‑301 (MQ=255)
                                                                                                                                                                                                        agaAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTa                                                                                              >  2:280270/1‑301 (MQ=255)
                                                                                                                                                                                                                gCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAAt                                                                                      <  2:84362/301‑1 (MQ=255)
                                                                                                                                                                                                                    tACGGCCTCTTTCAGGATAGAAGGATGTCAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCACCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTa                                                                                  <  2:340246/301‑1 (MQ=255)
                                                                                                                                                                                                                         ccTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTACGGTTACGTTCGCCGCAATGGTAAAAgg                                                                             >  1:509022/1‑301 (MQ=255)
                                                                                                                                                                                                                               tCAGGATAGATGGTTGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGAt                                                                       <  2:24569/301‑1 (MQ=255)
                                                                                                                                                                                                                                    ttAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTa                                                                  <  2:509694/300‑1 (MQ=255)
                                                                                                                                                                                                                                                       aatttaTAATTAATGTCGGCTTGTTGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCAc                                               <  1:546775/296‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTAATTTATTTTGTCGTGATGTCTTTTTCATTGACAATCTTGCCAGCGTTGTTGGGCTGTTATTAATCGCCCAGTCATCCCCGCTCTCTTTCGCTTTTAACATTCGTTCATCAATGCTATCGGCAACAATAACAATTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTc                                         <  1:221048/301‑1 (MQ=255)
                                                                                                                                                                                                                                                            aTAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCAcagata                                          >  2:547448/1‑296 (MQ=255)
                                                                                                                                                                                                                                                               attaatGTCGGTGTTGTGCCTGTTTCAGTGAAAATCTCGCCGATCTGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAg                                       <  1:46111/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                ttaatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCACCAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCAACATAAACAAGAACGGTATCATTTTCAGTCGGTTGACCAGTTCCGGGTAAGTTAACCGCAATGGAAAAAGGATAAATCACTAACGAAAAATGATCGGCACCTACTTCAgg                                      >  1:337845/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 taatGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTAAGGa                                     >  1:293520/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                  aatGTCGGCTTGATTCCTTTTTCAGTGAAAATCTCGCCAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAAt                                   <  2:180802/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                       tgtcTTTATGTCTTTTTCATTGAAAATTTCGCCAAGCGTTTTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTc                               <  2:338308/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                       cGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGAAATTTTc                               >  2:269780/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                               tGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTAAGTCGGTTGACCATGTCCGGTTCAGTAAGCCGCAATGGTAAAAGGTTCGAACATTACCGCATAATTATCGGACACTGCTTAAGGAATTTTCCACGCAtt                       >  1:588443/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                  cTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACAAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCACATATTGAACGGCCACTGCTTCCAGAAATTTCCACGCAttttt                    >  1:114926/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                   tttttCAGTGAAATTTTTGCCATGCGTTGTCTGGCTGTTATTATTCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTCCACGCAtttttt                   <  1:279967/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                       tctgTGAAAATTTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTTAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGAGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGGCCGGTTACGTTAGCCGCAATGGGAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTCCACGCATTTTTTGCCg               <  2:96193/298‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCCCTCTCTTTCGCTTTTTCCAGTCGTTCCTCCCTGGGATCGGCCCCAATCACATTTTTAACTGTATAGACGCCTTTTCATCCCTGAACGCTCGTCAGGCCTCTTTTCCCTTCACCCTCCCCCCTCCCGGTCTCGTTTTCAGTCGGTTGCCCATCTCCGGTTCCTTTCTCAGCAAGTTTCCCCCCTGCGACCATTCCAGCATATTGATCAGACACTGCTTCAGGAATTTTCCACGCATTTTTTGCCGaaa            >  1:96069/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                ttctcGCCAAGCGGTGTTTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGATCGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTCCACGCATTTTTTGCCGGAACCACgg      <  2:11476/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                    cGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTCATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGCACGGCCACTGCTTCAGGAAATTTCCACGCCTTTTTTTCCGGAACCACCGcctt  >  1:458384/1‑298 (MQ=255)
                                                                                                                                                                                                                                                                                                    cGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTAAGGAATTTTCCACGCATTTTTTGCCGGAACCACGGCata  >  1:523583/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                                                                        
GCAACATGCTGGAAATCAAACGTATGGGTAATTAATTTTTCTGGTTTAATTAACCCTTTACTTAACCAGTCGATAACGATCGGGAATTTATTTGCATTTAAGCGTGAAGAGAAAATAGAGAGTTCTTTTCCGGTAATTCCTTGCTGAATCACTTCAGACGGTTCACTGGAGAACCCCATCAATACAATACGTGCCGCTGGAGAAGCCAGCGTTACGGCCTCTTTCAGGATAGAAGGATGACAAGCCGCATCGATAATTAATGTCGGCTTGATGCCTTTTTCAGTGAAAATCTCGCCAAGCGGTGTCTGGCTGTTATTAATCGCCCAGTCAGCCCCGCTCTCTTTCGCTTTTTCCAGTCGTTCATCAATGCGATCGGCAACAATCACATTTTTAACGTTATAGACGCCTTTTAATACCTGAACGATCGTCAGGCCGATTGGACCGGCACCATAAACCAGAACGGTATCATTTTCAGTCGGTTGACCATGTCCGGTTACGTTAGCCGCAATGGTAAAAGGTTCGATCATTACCGCATATTGATCGGCCACTGCTTCAGGAATTTTCCACGCATTTTTTGCCGGAACCACGGCATA  >  NC_000913/1652966‑1653558

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