Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,951,976 C→T intergenic (+1419/‑225) ilvX → / → ilvM uncharacterized protein IlvX/acetolactate synthase II subunit IlvM

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,951,9760CT98.6% 241.0 / NA 69intergenic (+1419/‑225)ilvX/ilvMuncharacterized protein IlvX/acetolactate synthase II subunit IlvM
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (31/37);  minor base A (1/0);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 4.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

CGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCAACATCAGGTCAATG  >  NC_000913/3951729‑3952222
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
cGGAAAATTTCATCACCTCCAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGc                                                                                                                                                                                                                                                     <  1:512591/251‑1 (MQ=255)
 ggAAAATTTCATCACCTCCAGCGGCTTAGGTACCATTGGTTGTGGTTTACCGGCGTCGGTTGGCGCGCAAGGCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGTACCGTAAAACGCAAGCATTTACC‑GTTTAAAATCTTCTTACTCGATTACCAACGGTTAGGGAGGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCg                                                                                                                                                                                                                                                    <  1:134587/251‑1 (MQ=255)
   aaaaTTTCATCACCTCCAGCGGCTGAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGaa                                                                                                                                                                                                                                                  <  2:351471/251‑1 (MQ=255)
                  ccAGCGGCTTAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGGCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGAt                                                                                                                                                                                                                                   <  2:307672/251‑1 (MQ=255)
                     ggggcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAAc                                                                                                                                                                                                                                <  2:223543/246‑1 (MQ=255)
                       ggcttaagcaCCATGGGTTTTGGTTTACCGGCGGCGGTTGGGGCGGAAGACCCGCGTCCAAACGTTACCGTGGGCGGTATCTCCGGGGCCCGATCTTTCATGTTGGATGTGGACAG‑GCTGGGCACCGGAAAAACCAAGCAGTTACCAGTGGAAAA‑CGGCTTACCCGATAACCCACGGTTAGGGGTGGTGTGACAATGGCCACAACTGTTTTTCCACGGACGAGATAGCGGAACCACCCTTACCGATacacc                                                                                                                                                                                                                              >  2:477413/10‑247 (MQ=255)
                        gcttAGGCACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGGCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGGGCA‑AGAGCTGGGCACCGTAAAAAGCAAGCAGTTACC‑GTTGGAAAACGGCTTACTCGATAACCAACGGTTAGGGGTGGGTTGAAAATGGCAACAAATGGTTTTCCAGGGACGATATAGCGAAAACACCCTTACCGATAAcccc                                                                                                                                                                                                                             >  2:153931/3‑251 (MQ=255)
                                   aTGGGTTTTGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCAt                                                                                                                                                                                                                  <  1:188442/251‑1 (MQ=255)
                                          ttGGTTTACCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGc                                                                                                                                                                                                            <  1:421657/250‑1 (MQ=255)
                                                 aCCGGCGGCGGTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGcc                                                                                                                                                                                                     <  1:476189/250‑1 (MQ=255)
                                                    ggcggcggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCg                                                                                                                                                                                                 <  2:178976/251‑1 (MQ=255)
                                                     gcggcggTTGGCGCGCAAGTCGCGCGACCAAACGATTCCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCTCCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCgg                                                                                                                                                                                                <  1:293882/251‑1 (MQ=255)
                                                          ggTTGGCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATcc                                                                                                                                                                                           >  1:128251/1‑251 (MQ=255)
                                                              ggCGCGCAAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTgg                                                                                                                                                                                       <  1:425910/251‑1 (MQ=255)
                                                                   gcaAGTCGCGCGACCAAACGATACCGTCGTCTGTATCTCCGGAGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGCGATGGTTCGACAATGGCAGCCACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCTTGTTAGCCAGCGCCTTCGGCATCCCAGGCCAac                                                                                                                                                                                  >  2:385785/2‑251 (MQ=255)
                                                                           cgcgACCAAACGATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCAcc                                                                                                                                                                          >  2:43731/1‑251 (MQ=255)
                                                                           cgcgACCAAACGATACCGTCGTCTGTATCGCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGTCCAACACATCAcc                                                                                                                                                                          <  1:263510/251‑1 (MQ=255)
                                                                                 caaaCGATACCGTCGCCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCgtaa                                                                                                                                                                     >  1:458041/3‑247 (MQ=255)
                                                                                      gATACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGAcc                                                                                                                                                                >  2:542309/1‑250 (MQ=255)
                                                                                       aTACCGTCGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCCACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCCCCCGTAACGAcc                                                                                                                                                                >  2:347633/1‑249 (MQ=255)
                                                                                            gtcgTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGa                                                                                                                                                         <  1:361016/248‑1 (MQ=255)
                                                                                               gTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCNAACACATCACCCGTAAAGACCAGGTTGAAGc                                                                                                                                                      >  2:270090/1‑251 (MQ=255)
                                                                                                    tATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCAc                                                                                                                                                 >  1:129408/1‑251 (MQ=255)
                                                                                                              gACGGATCTTCCATGATGAATGTGCA‑AGAGCTGGCCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATAACCCGTACAGACCAGGTTGAAGCGGCACTCGACACCAt                                                                                                                                       >  2:159228/1‑251 (MQ=255)
                                                                                                                 ggatttcgcATGATGGATGTGAA‑AGAGCTGGGCACTGTAAAACGCAAGCAGTGACC‑GGTGAAAACCGTCTGCCTCGATAACCAGAGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACTGATACCCCCGATTTCATCATGTGAGCCAGCGCCGTCGGCCTCCTTGGCCAACACATCCTCCGTAAAGACCAGGTTTCAGCGGCACTCGACACCATGCt                                                                                                                                    <  2:170828/243‑1 (MQ=255)
                                                                                                                          atgatgAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTg                                                                                                                            >  1:175640/1‑250 (MQ=255)
                                                                                                                          atgatgAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTg                                                                                                                            >  1:512417/1‑250 (MQ=255)
                                                                                                                               gaATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGc                                                                                                                      >  1:39964/1‑251 (MQ=255)
                                                                                                                                   gtgCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCAt                                                                                                                   <  1:176203/250‑1 (MQ=255)
                                                                                                                                              tGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCAt                                                                                                        >  1:225458/1‑251 (MQ=255)
                                                                                                                                                 gCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGtc                                                                                                     >  1:154595/1‑251 (MQ=255)
                                                                                                                                                     cGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCaa                                                                                                 <  1:416683/251‑1 (MQ=255)
                                                                                                                                                     cGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCa                                                                                                  >  1:411778/1‑250 (MQ=255)
                                                                                                                                                      gTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCCCAAt                                                                                                >  1:47677/1‑251 (MQ=255)
                                                                                                                                                        aaaaCGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATc                                                                                               <  1:211567/250‑1 (MQ=255)
                                                                                                                                                          aaCGCAAGCAGTTACC‑GTTGAAAATCGTCTTNCTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCTACTGGTTTGCCAGGAACGATATAGCGAAACCACCCTTGCCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTCCATGTCCCAATcgac                                                                                            <  2:234146/251‑1 (MQ=255)
                                                                                                                                                           aCGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATcgacg                                                                                           <  1:136572/251‑1 (MQ=255)
                                                                                                                                                              cAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAAc                                                                                        <  2:90575/251‑1 (MQ=255)
                                                                                                                                                               aaGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACt                                                                                       >  2:133702/1‑251 (MQ=255)
                                                                                                                                                                  cAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTg                                                                                     >  1:13321/1‑250 (MQ=255)
                                                                                                                                                                  cAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGAAAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCCTCCCTGGCCAACACATCACCCGTAAAGACCAGGATGAAGCGGCCCTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCCATCGACGAACTTg                                                                                     >  2:288810/1‑250 (MQ=255)
                                                                                                                                                                  cAGTTACC‑GTTGAAAATAGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATATTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGACACTCGACACCATTCTGACCAGTGATGGGCCATACCTGCTTCAGGTCTCAATCGACGCACTTg                                                                                     >  2:299726/1‑250 (MQ=255)
                                                                                                                                                                         c‑GTTGAACATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGTAACGATATAGCGAAACCCCCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACTCATTACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTc                                                                             <  1:275858/251‑1 (MQ=255)
                                                                                                                                                                         c‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTc                                                                             <  1:351651/251‑1 (MQ=255)
                                                                                                                                                                                 aaTCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGCTTTCCTCATGTTAGCCAGCGCCTTAGGCATCCCTGGCCAACACATCCCCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGCGCACGTCTGGCCGc                                                                      >  2:375745/1‑251 (MQ=255)
                                                                                                                                                                                    cGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTgg                                                                   <  1:288810/251‑1 (MQ=255)
                                                                                                                                                                                     gTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGt                                                                  >  2:166784/1‑251 (MQ=255)
                                                                                                                                                                                     gTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGt                                                                  >  1:427930/1‑251 (MQ=255)
                                                                                                                                                                                         tACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGCTTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTgccg                                                              >  1:293832/1‑251 (MQ=255)
                                                                                                                                                                                         tACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTgcc                                                               <  1:153931/250‑1 (MQ=255)
                                                                                                                                                                                               aTAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTggg                                                        >  1:335705/1‑250 (MQ=255)
                                                                                                                                                                                               aTAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTggg                                                        >  1:142694/1‑250 (MQ=255)
                                                                                                                                                                                                     accGGTTAGGGATGGTTCGACAATGGCAGCAACGGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGt                                                  <  2:405724/249‑1 (MQ=255)
                                                                                                                                                                                                      aCGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGATCGATATAGCGAAACCTCCCTTACCGATAACCCCGTTTTCCTCATGTTAGCCAGCGCCTTCGCCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGCCACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGATCTTTAGAACGTCTTGCCGCTTGTGCCGCCTGGTGCCAGTa                                                 <  1:125082/251‑1 (MQ=255)
                                                                                                                                                                                                         gTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAAtt                                              <  1:170384/251‑1 (MQ=255)
                                                                                                                                                                                                             gggATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGGGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGa                                          <  2:494076/251‑1 (MQ=255)
                                                                                                                                                                                                              ggATGGTTCGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGa                                          >  1:85787/1‑250 (MQ=255)
                                                                                                                                                                                                                     gcGACAATGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCAGACCTGCTTCATGTCTCAAGCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTgg                                  <  2:489178/250‑1 (MQ=255)
                                                                                                                                                                                                                          aaTGGCAGCAACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGaaa                             <  1:542309/251‑1 (MQ=255)
                                                                                                                                                                                                                               cagcaACTTGTTTTCCAGGAACGATATAGCGAAACCCCCCTTACCGATACCCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCCTCCCTGGCCAACACATCACCCGTAAAGACCAGGTTTAAGCGGCACTCGACACCATGCTGAACAGTGATGGTCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGGAATTCAGAAATGTTGGAGAAAGTATc                        <  2:350562/251‑1 (MQ=255)
                                                                                                                                                                                                                                agcaACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCa                       <  1:181873/251‑1 (MQ=255)
                                                                                                                                                                                                                                  caACTGTTTTTCCAGTAACGAGATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGGCCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGTGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTTGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatg                     <  1:72534/251‑1 (MQ=255)
                                                                                                                                                                                                                                  caACTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatg                     <  1:292776/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatgat                   <  1:375745/251‑1 (MQ=255)
                                                                                                                                                                                                                                    aCTGTTTTTCCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACAGCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGGCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCatgat                   <  2:227192/251‑1 (MQ=255)
                                                                                                                                                                                                                                        tttttcCAGGAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCAGCTTCATGTCTCAATCTACGAACTTGAGAACGTCTGGCCGATGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCa                >  2:79916/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 gAACGATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCAACATCAGGTc      <  1:80256/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   aCGATAAAGCGAAACCAACATGACCGGTAAACCCGATTTCCTCATGATAGCCAGCGCCTTCGGCATCCCTGGCCAACACATTAACCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTATTTAGGAACAGTTGGAGAAATTATCAAGATGCAACATCAGGTCa     >  2:79925/1‑250 (MQ=255)
                                                                                                                                                                                                                                                     gATATAGCGAAACCACCCTTACCGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCCTGGCCAACACATCACCCGTAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGTGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCAACATCAGGTCAATg  >  1:507967/1‑251 (MQ=255)
                                                                                                                                                                                                                                                         |                                                                                                                                                                                                                                                      
CGGAAAATTTCATCACCTCCAGCGGTTTAGGTACCATGGGTTTTGGTTTACCGGCGGCGGTTGGCGCACAAGTCGCGCGACCGAACGATACCGTTGTCTGTATCTCCGGTGACGGCTCTTTCATGATGAATGTGCA‑AGAGCTGGGCACCGTAAAACGCAAGCAGTTACC‑GTTGAAAATCGTCTTACTCGATAACCAACGGTTAGGGATGGTTCGACAATGGCAGCAACTGTTTTTTCAGGAACGATACAGCGAAACCACCCTTACTGATAACCCCGATTTCCTCATGTTAGCCAGCGCCTTCGGCATCCATGGCCAACACATCACCCGGAAAGACCAGGTTGAAGCGGCACTCGACACCATGCTGAACAGTGATGGGCCATACCTGCTTCATGTCTCAATCGACGAACTTGAGAACGTCTGGCCGCTGGTGCCGCCTGGCGCCAGTAATTCAGAAATGTTGGAGAAATTATCATGATGCAACATCAGGTCAATG  >  NC_000913/3951729‑3952222

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