Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 5,067 Δ1 bp intergenic (+47/‑167) thrC → / → yaaX threonine synthase/DUF2502 domain‑containing protein YaaX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009135,0670C.97.9% 202.5 / ‑3.5 47intergenic (+47/‑167)thrC/yaaXthreonine synthase/DUF2502 domain‑containing protein YaaX
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (21/25);  total (21/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/4832‑5306
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
ccAGGCAAATATGGTTTGCTCCTCGGCACGGCGCACCCGGCAAAATTAAAAGAGAGCGGGGAAGAGTTTCTCTGTGAAACTGTGGATCTCCCAAAAGATCTGGCATAACGTGCTGATTTGCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATgaag                                                                                                                                                                                                                                   <  2:404046/250‑1 (MQ=255)
                  ttCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGaaaaacg                                                                                                                                                                                                                 >  1:236404/1‑248 (MQ=255)
                       gggccccgcgctcccgggaaaaTTTAAAGGAGGGGTGGAAGCGTTTCCCGGGGAACCGTTGGATCGGCCAAAAGAGCTGGAAGAACGGTCTGATTTACCCTTTCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAggg                                                                                                                                                                                                           <  2:475593/232‑1 (MQ=255)
                           aCCGCGCATCCGGCGAACTTTAAAGAGAGCGTGGAAGCGATTCTCGGGGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGaaa                                                                                                                                                                                                        >  2:529777/1‑250 (MQ=255)
                                cctgccGGCGAAATTAAAAGAGGGTGGGGAAGCTGTTCCCGGTAAAACGTTGTACCTGCCAAAAGGGCGTGCAGAACTTGGGGATTATCCCTTTTTTTCACAAAATCGTCCCGCCGAGTTTGCTGCGTTGCGTAAATTGATGATGAATTATCATTAAAATCTATTCATTAGCTCAATCAGGTCGGGTTTGCTTTTATGCATCCGGCTTTTTTATGAAGAAATTATGAAGAAAAACGACATGGAAAgaggag                                                                                                                                                                                                  <  2:605638/247‑1 (MQ=255)
                                  aTCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGgaga                                                                                                                                                                                                 >  2:132884/1‑250 (MQ=255)
                                              ttAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACTTTGGATCTGCCAAAAGAGCTGGCAGAACTTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATaa                                                                                                                                                                                     <  1:65467/250‑1 (MQ=255)
                                                    agaaggTGGAGGCAATTCTCTGTGAACCGTTGGATCTGCAAAAAGAGCTGGCAAAAAGTGGTGTTTATCCCGTGTTTTCACAAAAACTGACCGCCGATGTTTCTGCGTTGCGTAAATTGATGATGAATCATTAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCAGGATTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCGCAAGAAATGCGGt                                                                                                                                                                              <  2:198663/246‑1 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAt                                                                                                                                                                                   >  2:373966/1‑240 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAt                                                                                                                                                                                   <  1:373966/240‑1 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTATCACTAAATGCGGTAACtt                                                                                                                                                                         >  2:402884/1‑250 (MQ=255)
                                                               aGCGATCCTCGGGGAAACGTTGTATCGGACAAAAGAGCTGGCAGAACGTGCTGATTGCCCCTTGCTTACACATACTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAACTTCTCAATAAATGCCGTAACTTAgaga                                                                                                                                                                    <  2:217531/251‑1 (MQ=255)
                                                                   aTTCTCGGTGAAACGTTGGATCGGCAAAAAGAGCTGGCAGAACGGGCTGATTTACCCTTGCTTTCACATAACTTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAg                                                                                                                                                                <  1:342785/251‑1 (MQ=255)
                                                                                cGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCAAATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCGGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGCAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAAt                                                                                                                                                   >  1:16018/1‑251 (MQ=255)
                                                                                          gCCAAAAGGGCGGGCCGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             <  2:521626/247‑1 (MQ=255)
                                                                                                 gagCTGGCAGAACGTGTTGATTTTCCCTTTTTTTCTCAATATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTATATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAGGCGGTAACTTAGAGAGTAGGAGTGCGGTGAATAACAACCGTCGTTCTCa                                                                                                                                  <  2:315340/251‑1 (MQ=255)
                                                                                                              gTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTcc                                                                                                                     >  2:525338/1‑251 (MQ=255)
                                                                                                                    aTTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGatat                                                                                                               >  2:473191/1‑251 (MQ=255)
                                                                                                                        aCCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAc                                                                                                           <  1:409543/251‑1 (MQ=255)
                                                                                                                           ccTTGTTTCCACAAAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCAGCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGTAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAccc                                                                                                         <  2:623051/251‑1 (MQ=255)
                                                                                                                              tgtttttaaaaagtgtggccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAAGCATCAGTAAAATCTATTCATTATCACAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAGCGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGGGAGTTGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCAt                                                                                                       <  2:236404/232‑1 (MQ=255)
                                                                                                                                        aaaTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATGAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATGACGGGCAATg                                                                                            <  2:595073/250‑1 (MQ=255)
                                                                                                                                          aTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          >  1:537329/1‑251 (MQ=255)
                                                                                                                                              gCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                      <  1:558458/251‑1 (MQ=255)
                                                                                                                                                ccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAATCATTATCTCAATCAGGCCGGGTGTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    <  1:229716/251‑1 (MQ=255)
                                                                                                                                                   ccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAACCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                 <  1:473191/251‑1 (MQ=255)
                                                                                                                                                     gATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                >  1:581382/1‑250 (MQ=255)
                                                                                                                                                      aTTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             >  1:284484/1‑188 (MQ=255)
                                                                                                                                                      aTTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAGaacaac                                                                                                                                             <  2:284484/188‑1 (MQ=255)
                                                                                                                                                        tttGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCt                                                                           >  1:487551/1‑251 (MQ=255)
                                                                                                                                                          tgctgcGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCtat                                                                         >  1:329579/1‑249 (MQ=255)
                                                                                                                                                                   gcgtAAATTGAGGAGGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAATAACGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTGCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGAt                                                                <  1:36974/251‑1 (MQ=255)
                                                                                                                                                                     gtAAATTGATGATGAATAATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGc                                                              >  2:87016/1‑251 (MQ=255)
                                                                                                                                                                       aaaTTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCa                                                             >  1:346485/1‑250 (MQ=255)
                                                                                                                                                                            gatgatgaATCATCAGTAAAAACTATTCATTATCTCAATCAGGCCGGGTTTGCTTGTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCAACGCGTAATCTCCGGATAACGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATcga                                                        <  2:525317/250‑3 (MQ=255)
                                                                                                                                                                               gatgaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCg                                                     <  1:400548/250‑1 (MQ=255)
                                                                                                                                                                                  gaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTAc                                                  <  1:164638/250‑1 (MQ=255)
                                                                                                                                                                                      catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCa                                             <  2:166611/251‑1 (MQ=255)
                                                                                                                                                                                         caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACtt                                          <  2:2479/251‑1 (MQ=255)
                                                                                                                                                                                               actCTAGTCATTATCTCACTCAGGCCGGGTTTACTTATATGCAGCCGGCTTTTTTTTCAAGAAATTATGGAGAAAAACGACAGGGAAATAGGTGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGAGAATAACAACCGCCGTACTCATCGCGTAATCTCCGGATATAAACCCATAACGGGCAATGCTAACAGGATCACCATATGAAAAAGCTGAAATCTATCGTACTCGAACTTTCCCTg                                    >  2:109328/3‑251 (MQ=255)
                                                                                                                                                                                                aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                   >  2:178695/1‑251 (MQ=255)
                                                                                                                                                                                                      tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                             >  1:411235/1‑251 (MQ=255)
                                                                                                                                                                                                        aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCGGCTTGTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                           <  2:212253/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGCCGGGTTTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                >  2:358045/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccGGGTGTGCTTTTATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGATAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAATATGCAATCTATCGTACTCGCACTTTCCCTCGTTCTGGCCGCTCCCATGgcagca            <  2:207439/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttttATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:448274/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tttATGCAGCCGGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:329905/1‑250 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
CCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCGTGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTACCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACCTGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/4832‑5306

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