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,4980CG97.0% 234.8 / 0.9 69coding (494/891 nt)alsRDNA‑binding transcriptional repressor AlsR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  major base G (38/27);  minor base . (1/1);  total (39/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.15e-01

GCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTTCC  >  NC_000913/4312259‑4312743
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
gCGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTtg                                                                                                                                                                                                                                            <  1:531811/251‑1 (MQ=255)
 cGTTTCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTtgt                                                                                                                                                                                                                                           >  1:200049/1‑251 (MQ=255)
   tttCCGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGc                                                                                                                                                                                                                                         >  1:250278/1‑251 (MQ=255)
      ccGGGGCTGGTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTg                                                                                                                                                                                                                                       <  1:202600/250‑1 (MQ=255)
           gCTGGGGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATTGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATc                                                                                                                                                                                                                                 <  1:43359/251‑1 (MQ=255)
            cTGGTGAGCAAATAATATGATCGGCCCGTTTCGCTATCGGTGAATGGTAGCTATGAGTTATACAAATCATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGAATCAGGATAGGCCTGACAGCGTACGCCAATTGGCAAGAACTTGTGCTGTACATc                                                                                                                                                                                                                                 >  2:406706/1‑250 (MQ=255)
               gTGAGCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGTATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGc                                                                                                                                                                                                                              >  1:210807/1‑250 (MQ=255)
                   gCAAATAATATGATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGcac                                                                                                                                                                                                                            >  1:167336/1‑248 (MQ=255)
                         aatatgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCAtt                                                                                                                                                                                                                  >  2:177508/1‑251 (MQ=255)
                          atatgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCAGAGATAGCAtt                                                                                                                                                                                                                  <  1:82071/251‑1 (MQ=255)
                            atgatCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTg                                                                                                                                                                                                                >  1:369744/4‑251 (MQ=255)
                              gatCGGCCAGTTTCGCTATCGGTGGATGTGAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGAt                                                                                                                                                                                                              <  2:242084/250‑1 (MQ=255)
                                           cGCTATCGGAGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGACTGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATCGCATTTGATCCTCCGGCACCGt                                                                                                                                                                                                 >  2:218378/1‑251 (MQ=255)
                                            gCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCCAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGtc                                                                                                                                                                                                >  2:632513/1‑250 (MQ=255)
                                                   gTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATcc                                                                                                                                                                                         >  1:166028/1‑251 (MQ=255)
                                                             gCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCCCTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTgg                                                                                                                                                                               <  2:20585/251‑1 (MQ=255)
                                                               tATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGcc                                                                                                                                                                             >  1:491851/1‑251 (MQ=255)
                                                                      tataCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGAt                                                                                                                                                                         >  2:428302/1‑248 (MQ=255)
                                                                      tataCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAg                                                                                                                                                                       <  1:615548/250‑1 (MQ=255)
                                                                       ataCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGa                                                                                                                                                                      >  1:5717/1‑250 (MQ=255)
                                                                       ataCAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGa                                                                                                                                                                      >  1:316548/1‑250 (MQ=255)
                                                                          cAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGaaaaa                                                                                                                                                                  <  2:250278/251‑1 (MQ=255)
                                                                          cAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGaaaa                                                                                                                                                                   >  1:408333/1‑250 (MQ=255)
                                                                           aaaTAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCAGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGaaaa                                                                                                                                                                   >  2:67290/1‑249 (MQ=255)
                                                                                   tttGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGggcggc                                                                                                                                                         >  2:405674/1‑251 (MQ=255)
                                                                                              cTTTTTTGCCAGTTCTATGGCCGCTTTTACATCAATGGGTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATACTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtct                                                                                                                                              >  2:152476/1‑251 (MQ=255)
                                                                                               ttttttGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtctc                                                                                                                                             >  1:91751/1‑251 (MQ=255)
                                                                                                tttttGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGAtct                                                                                                                                              <  2:408333/249‑1 (MQ=255)
                                                                                                       cAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACgg                                                                                                                                     >  2:514395/1‑250 (MQ=255)
                                                                                                          ttCTACGTCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACCACATCTCCTTCCTGTAACACCGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCGGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCAGCTGTCTGGCCTGATAGAAAACGCGGGCGGCACGGTGGTTCTCATCAACGtcg                                                                                                                                   >  2:151251/1‑248 (MQ=255)
                                                                                                                         tACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCt                                                                                                                    <  1:387337/250‑1 (MQ=255)
                                                                                                                         tACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCt                                                                                                                    <  1:225766/250‑1 (MQ=255)
                                                                                                                          aCATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCATTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGCCGATCGACTGACCTTc                                                                                                                  >  2:486409/1‑251 (MQ=255)
                                                                                                                              cgcTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAGCTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGGCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAGCGTTGACGATCGACTGACTTTCCAt                                                                                                               <  2:507263/248‑1 (MQ=255)
                                                                                                                                 tGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATg                                                                                                           >  2:88007/1‑251 (MQ=255)
                                                                                                                                  ggTTCTCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAAGAACGAAGCGGACATCAGCATGATGTGAGAATCAGGATAGGCCTGACAGCGTACGACAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTTATCCTCCGGCACCGTACAAATCCCGCTGTCTTGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACTTTGACGATCGACTGACCTTCCAGAAt                                                                                                            <  1:315112/249‑1 (MQ=255)
                                                                                                                                    ttCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTg                                                                                                        >  1:165053/1‑251 (MQ=255)
                                                                                                                                       gCCCGGAATGGGTCACTACCAGCAAAACATCTCCTTCCTGTCACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGGACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGt                                                                                                     >  2:207778/1‑251 (MQ=255)
                                                                                                                                         ccGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTaa                                                                                                   <  2:263229/251‑1 (MQ=255)
                                                                                                                                           ggAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTaa                                                                                                   >  1:596734/1‑249 (MQ=255)
                                                                                                                                              aTGGGTAACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTg                                                                                               >  2:188339/1‑250 (MQ=255)
                                                                                                                                                 ggTCACTACCAGCACAACATCTCCTTCCTGTACCAACGAAGCGGACATCATCAAGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATTGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCCCCGGCCCCGTTCACATCCCGCTGTCTGGCCTTATAGAAAAAGCGGGTGGCACGGTGGATCTCATCAACGTTGACGATCGACTGCCCTACCATAATGGTGCGTAAAGTGATGt                                                                                           <  1:425307/251‑1 (MQ=255)
                                                                                                                                                 ggTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGt                                                                                           <  1:510699/251‑1 (MQ=255)
                                                                                                                                                  gTCACTACTAGTACACCATCTCCTTCGTGTAACAACGGAGCGGACAGCAGCATGATGTGAGCATCAGGATAGGCCTGACAGAGCACGCCAATGCGCAAGAACGTGTGCTGTACATCAGCACAGAGAGCTTTTGATCCTCCGGCACGGTACAAAGCCGGCTGTCTGGCCTGATAGAAAAAGCGGGCTGCACGTTGGATCTCAACAAGGTTGACGAGCGTCTGACTTTCCACAGGGGTGCGTAAAGTGATGtc                                                                                          <  2:16111/251‑2 (MQ=255)
                                                                                                                                                        aCCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAAt                                                                                      <  1:34038/249‑1 (MQ=255)
                                                                                                                                                          cAGCACCACATGCCCTTCCTGTAACAACGTTGCGCACATCATCATGATGTGAGCCTCAGGATCGGCCTGACACCGTCCGCCCATGGGCGAGAACTTGTGCTGTCCAGCAGCACAGATAGCATTTGATCCTCCGGCACCGCACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGCCACGGTGGATCGCCTCCACGTTGGCGATGGACTGACGTTACATAATGGTGCGTCAAGTGATGTTAAATatct                                                                                  >  2:447605/1‑248 (MQ=255)
                                                                                                                                                            gCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTAAAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCtt                                                                                 >  1:481912/1‑250 (MQ=255)
                                                                                                                                                             cacaACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCTGCACCGTACAAACCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCTACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTAt                                                                               <  1:151334/251‑1 (MQ=255)
                                                                                                                                                                 acaTCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGCTAGGCCTGACAGCGTACGCCAATGGACAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGGGCGTAAAGTGATGTTAAATACCTTATTcac                                                                           >  2:219863/1‑251 (MQ=255)
                                                                                                                                                                          ttctGTAACAACGAAGCTGACATCATCATGATGTGAGCGTCCGGATAGGCGTGACAGCGTACGCCAAGGTGCAAGATCTTGTGCTGTCCATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATGTCATCAACGTTGACGATCGACTGACCTGCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATGCGCCACATcc                                                                    <  2:414715/247‑1 (MQ=255)
                                                                                                                                                                              gTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTgcggc                                                              <  1:428302/251‑1 (MQ=255)
                                                                                                                                                                                  caacGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTtc                                                         <  1:177508/251‑1 (MQ=255)
                                                                                                                                                                                               atcatcatGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTGCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCTCTTCATCAAAAGCCaa                                             <  2:491851/251‑1 (MQ=255)
                                                                                                                                                                                                 catcatGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGAAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACt                                           >  2:193592/1‑251 (MQ=255)
                                                                                                                                                                                                    catGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCgg                                        >  2:492338/1‑251 (MQ=255)
                                                                                                                                                                                                        atgTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAgg                                    >  1:290648/1‑251 (MQ=255)
                                                                                                                                                                                                          gtgAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGc                                   <  2:316548/250‑1 (MQ=255)
                                                                                                                                                                                                            gAGCATCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAAt                                <  2:200049/251‑1 (MQ=255)
                                                                                                                                                                                                                aTCAGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCt                            >  2:431579/1‑251 (MQ=255)
                                                                                                                                                                                                                   aGGATAGGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGtt                         >  2:539772/1‑251 (MQ=255)
                                                                                                                                                                                                                    ggATATGCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATTCCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTc                        <  1:406706/251‑1 (MQ=255)
                                                                                                                                                                                                                         ggCCTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGAGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTctgact                   >  2:453709/1‑251 (MQ=255)
                                                                                                                                                                                                                          gCCTGACAGCGTACGCCAATGGGCAAGAACGTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTctgactg                  <  1:426675/251‑1 (MQ=255)
                                                                                                                                                                                                                           ccTGACAGCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTctgactgc                 >  2:524144/1‑250 (MQ=255)
                                                                                                                                                                                                                                  gCGTACGCCAATGGTCAAGAACTTGTGCTGTACATCAGCGCAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAaat          <  1:132120/251‑1 (MQ=255)
                                                                                                                                                                                                                                  gCGTACGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAaat          <  2:69442/251‑1 (MQ=255)
                                                                                                                                                                                                                                       cGCCAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCt     >  1:50286/1‑251 (MQ=255)
                                                                                                                                                                                                                                          cAATGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTTcc  <  1:431579/251‑1 (MQ=255)
                                                                                                                                                                                                                                           aaTGGGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTTcc  >  2:340837/1‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                     
GCGTTTCCGGGGCTGGTGAGCAAATAATATAATCGGCCAGTTTCGCTATCGGTGAATGGTAGCTATGGGTTATACAAATAATCTTTGCCCCGTTCTTTTTTGCCAGTTCTACGGCCGCTTTTACATCACTGGTTCGCCCGGAATGGGTCACTACCAGCACAACATCTCCTTCCTGTAACAACGAAGCGGACATCATCATGATGTGAGCATCAGGATAGGCCTGACAGCGTACGCCAATGCGCAAGAACTTGTGCTGTACATCAGCACAGATAGCATTTGATCCTCCGGCACCGTACAAATCCCGCTGTCTGGCCTGATAGAAAAAGCGGGCGGCACGGTGGATCTCATCAACGTTGACGATCGACTGACCTTCCATAATGGTGCGTAAAGTGATGTTAAATACCTTATTCACCACATCCTGCGGCGCTTCATCAAAAGCCAACTCGGAAGGCAATACCTGTTCTGACTGAGAAAAATAATCTTCC  >  NC_000913/4312259‑4312743

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