Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,312,498 C→G R165P (CGC→CCC)  alsR ← DNA‑binding transcriptional repressor AlsR

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,312,4980CG98.2% 201.2 / NA 55R165P (CGC→CCC) alsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base G (28/26);  total (28/27)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTT  >  NC_000913/4312260‑4312741
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
cGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTtgt                                                                                                                                                                                                                                         <  2:331655/251‑1 (MQ=255)
        gggCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTAc                                                                                                                                                                                                                                  >  1:243190/1‑250 (MQ=255)
              gTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGca                                                                                                                                                                                                                           <  2:3716/251‑1 (MQ=255)
                gAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGCAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGcac                                                                                                                                                                                                                          >  1:356089/1‑250 (MQ=255)
                        aatatgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCa                                                                                                                                                                                                                  <  1:250697/250‑1 (MQ=255)
                         atatgatCGGCCAGTTTGGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTGTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCGAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAtt                                                                                                                                                                                                                <  2:232325/251‑1 (MQ=255)
                           atgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAttt                                                                                                                                                                                                               <  1:418720/247‑1 (MQ=255)
                                    cAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTAAGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCATCCgg                                                                                                                                                                                                     >  2:434536/1‑251 (MQ=255)
                                      gTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCa                                                                                                                                                                                                   <  1:6311/251‑1 (MQ=255)
                                              aTCGGTGCCTGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCCGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACa                                                                                                                                                                                            >  1:306917/1‑250 (MQ=255)
                                                 ggTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATc                                                                                                                                                                                        >  1:399819/1‑251 (MQ=255)
                                                   tgactGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATccc                                                                                                                                                                                      >  1:58840/5‑251 (MQ=255)
                                                   tGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTCCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATccc                                                                                                                                                                                      >  2:345508/1‑251 (MQ=255)
                                                        ggtggcaaTGGGTTATACAAATAATCTGTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTTCCAGCACAACATCTCCTGCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCATGAACGTGTGCTGTACATCAGCACAGATAGCATTTGATCCTTCGTCACCGTGCAAATCCCGCTg                                                                                                                                                                                  <  2:388248/243‑1 (MQ=255)
                                                                 gggTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAttt                                                                                                                                                                                                               >  2:181459/1‑212 (MQ=255)
                                                                 gggTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAttt                                                                                                                                                                                                               <  1:181459/212‑1 (MQ=255)
                                                                      ataCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGAAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGa                                                                                                                                                                    >  2:426349/1‑250 (MQ=255)
                                                                           aataatCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGTTAGaaaaacc                                                                                                                                                              >  2:83940/1‑249 (MQ=255)
                                                                            ataatCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAACTCCCGCTGTCTGGCCTGATAGAAAAAgcg                                                                                                                                                             >  2:115719/1‑251 (MQ=255)
                                                                                     gCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCAc                                                                                                                                                     <  1:300128/250‑1 (MQ=255)
                                                                                         cGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGGTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTgg                                                                                                                                                <  2:429476/251‑1 (MQ=255)
                                                                                               tttttGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtct                                                                                                                                            <  1:104992/249‑1 (MQ=255)
                                                                                               tttttGCCAGTTCTACGGCCGCTTTTACAGCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATTGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtct                                                                                                                                            <  1:105452/249‑1 (MQ=255)
                                                                                                  ttGCCAGTACTACGGACGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCCCAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCtcatc                                                                                                                                        >  2:317716/1‑250 (MQ=255)
                                                                                                       aGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGtt                                                                                                                                  >  1:17289/1‑251 (MQ=255)
                                                                                                              cGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATc                                                                                                                           >  2:263531/1‑251 (MQ=255)
                                                                                                                  cGCTTTTACATCACTGGTTAGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGCAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGAACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGACCGGCGGATCTCATCAACGTTGACGATCGACt                                                                                                                       >  1:405735/1‑251 (MQ=255)
                                                                                                                   gCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGACAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTg                                                                                                                      >  1:87685/1‑251 (MQ=255)
                                                                                                                                       cccGGAATGGGTCACTACCAGCACAACATCTCCTTCCTATAACAACGAAGCGGACATCATCATGATGTGAGCATCATGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTa                                                                                                  <  2:17385/251‑1 (MQ=255)
                                                                                                                                          ggAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGAACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTaa                                                                                                 >  1:164227/1‑249 (MQ=255)
                                                                                                                                           gAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGt                                                                                              <  1:115719/251‑1 (MQ=255)
                                                                                                                                                 gTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGCTAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCACATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGtt                                                                                        >  2:230640/1‑251 (MQ=255)
                                                                                                                                                          aGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCtt                                                                               >  1:335086/1‑251 (MQ=255)
                                                                                                                                                             acaacaTCTCCATCCGGTAACAACGAAGCGTACAACAGCAGGGTGTGATCTTCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACGTGTGCTGTACATCAACACAGATAGCATTTGATCCACCGGCGCCGTACAAAACCCGCTGTCGGGCCTGATAGAAAGAGCGGGCGGGACGATGGATGTCATAAACGTGGACTATCGACTGACCTTCCATGATGGTGCGTAAATGGAGGTTAAATACCTTact                                                                            <  2:390129/251‑3 (MQ=255)
                                                                                                                                                               aacaTCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTTGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTca                                                                          <  2:200584/251‑1 (MQ=255)
                                                                                                                                                                acaTCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTcac                                                                         >  1:96178/1‑251 (MQ=255)
                                                                                                                                                                         tccgtgaacaacGAAGCGGACAACATCAGGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAGTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTTACGATCGACTGTCCTTCCATAATGGTGCGTAAAGTGATGTTAAATTCCTTATGCATCACATCCt                                                                 <  2:306917/244‑1 (MQ=255)
                                                                                                                                                                                     gAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTtcatc                                                    >  1:258230/1‑251 (MQ=255)
                                                                                                                                                                                          ggACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCaaa                                                 >  1:174610/1‑249 (MQ=255)
                                                                                                                                                                                             catcatcatGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCa                                            <  1:345508/251‑1 (MQ=255)
                                                                                                                                                                                               acatcatGATGTGAGCATCAGGATAGGCCTGACAGCGTTCGCAAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCa                                            <  2:449921/248‑1 (MQ=255)
                                                                                                                                                                                                catcatGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACt                                         >  1:449654/1‑251 (MQ=255)
                                                                                                                                                                                                      gatgTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGaa                                    <  1:410138/250‑1 (MQ=255)
                                                                                                                                                                                                       atgTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAgg                                  >  1:14624/1‑251 (MQ=255)
                                                                                                                                                                                                           gcgCATCAGGAGAGGCCTGACAGCGTACGCGAATGGGCAAGATCTTCTGCTGTACATCAGCACAGATAGCATTGAATCCTCCTGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCTCGGTGGATCTCATCAACGTTTACGATCGACTGACCTTCCATAATGGTGCGTAAAGATATGTTAAATACCTTATTCACCACGTCCTGCGGCGGTGCAGCAAAATCCAGCTCTGAAGGCAAt                              <  2:7788/249‑1 (MQ=255)
                                                                                                                                                                                                           gAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGCCGATAGACTGACCTTCCATAATGGTGCGTAAAGGGATGTTAAATACCTTATTCACCACATCATGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAAt                              >  2:417143/1‑251 (MQ=255)
                                                                                                                                                                                                           gAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAAt                              <  2:335086/251‑1 (MQ=255)
                                                                                                                                                                                                            aGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATa                             <  1:11181/251‑1 (MQ=255)
                                                                                                                                                                                                                 cAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGt                        >  1:298553/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTc                      <  1:220473/251‑1 (MQ=255)
                                                                                                                                                                                                                    gATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTct                     <  1:263531/251‑1 (MQ=255)
                                                                                                                                                                                                                       aGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAACCCCGCTGTCGGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTctgac                  >  1:10037/1‑251 (MQ=255)
                                                                                                                                                                                                                          ccTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTctgactg                >  1:270732/1‑250 (MQ=255)
                                                                                                                                                                                                                                     aCGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAaataat     <  1:194758/250‑1 (MQ=255)
                                                                                                                                                                                                                                       gCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCtt  <  2:399819/251‑1 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                   
CGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTT  >  NC_000913/4312260‑4312741

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