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

CCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCG‑TGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  NC_000913/4832‑5306
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
ccAGGCAAATATGGTTTGCTCCTCGGCACGGCGCACCCGGCAAAATTAAAAGAGAGCG‑GGGAAGAGTTTCTCTGTGAAACTGTGGATCTCCCAAAAGATCTGGCATAACGTGCTGATTTG‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATgaag                                                                                                                                                                                                                                   <  2:404046/250‑1 (MQ=255)
                  ttCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCG‑TGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGaaaaacg                                                                                                                                                                                                                 >  1:236404/1‑248 (MQ=255)
                       gggccccgcgctcccgggaaaaTTTAAAG‑GAGGGGTGGAAGCGTTTCCCGGGGAACCGTTGGATCGGCCAAAAGAGCTGGAAGAACGGTCTGATTTA‑CCCTTTCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAggg                                                                                                                                                                                                           <  2:475593/232‑1 (MQ=255)
                           aCCGCGCATCCGGCGAACTTTAAAGAGAGCG‑TGGAAGCGATTCTCGGGGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGaaa                                                                                                                                                                                                        >  2:529777/1‑250 (MQ=255)
                                cctgccGGCGAAATTAAAAGAGGGTG‑GGGAAGCTGTTCCCGGTAAAACGTTGTACCTGCCAAAAGGGCGTGCAGAACTTGGGGATTAT‑CCCTTTTTTTCACAAAATCGTCCCGCCGAGTTTGCTGCGTTGCGTAAATTGATGATGAATTATCATTAAAATCTATTCATTAGCTCAATCAGGTCGGGTTTGCTTTTATGCA‑TCCGGCTTTTTTATGAAGAAATTATGAAGAAAAACGACATGGAAAgaggag                                                                                                                                                                                                  <  2:605638/247‑1 (MQ=255)
                                  aTCCGGCGAAATTTAAAGAGAGCG‑TGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGgaga                                                                                                                                                                                                 >  2:132884/1‑250 (MQ=255)
                                              ttAAAGAGAGCG‑TGGAAGCGATTCTCGGTGAAACTTTGGATCTGCCAAAAGAGCTGGCAGAACTTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATaa                                                                                                                                                                                     <  1:65467/250‑1 (MQ=255)
                                                    agaagg‑TGGAGGCAATTCTCTGTGAACCGTTGGATCTGCAAAAAGAGCTGGCAAAAAGTGGTG‑TTTATCCCGTGTTTTCACAAAAACTGACCGCCGATGTTTCTGCGTTGCGTAAATTGATGATGAATCATTAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCA‑GGATTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCGCAAGAAATGCGGt                                                                                                                                                                              <  2:198663/246‑1 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAt                                                                                                                                                                                   >  2:373966/1‑240 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAt                                                                                                                                                                                   <  1:373966/240‑1 (MQ=255)
                                                           tGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTATCACTAAATGCGGTAACtt                                                                                                                                                                         >  2:402884/1‑250 (MQ=255)
                                                               aGCGATCCTCGGGGAAACGTTGTATCGGACAAAAGAGCTGGCAGAACGTGCTGATTGC‑CCCTTGCTTACACATACTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAACTTCTCAATAAATGCCGTAACTTAgaga                                                                                                                                                                    <  2:217531/251‑1 (MQ=255)
                                                                   aTTCTCGGTGAAACGTTGGATCGGCAAAAAGAGCTGGCAGAACGGGCTGATTTA‑CCCTTGCTTTCACATAACTTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAg                                                                                                                                                                <  1:342785/251‑1 (MQ=255)
                                                                                cGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCAAATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCGGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGCAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAAAAt                                                                                                                                                   >  1:16018/1‑251 (MQ=255)
                                                                                          gCCAAAAGGGCGGGCCGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             <  2:521626/247‑1 (MQ=255)
                                                                                                 gagCTGGCAGAACGTGTTGATTTT‑CCCTTTTTTTCTCAATATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTATATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAAGGCGGTAACTTAGAGAGTAGGAGTGCGGTGAATAACAACCGTCGTTCTCa                                                                                                                                  <  2:315340/251‑1 (MQ=255)
                                                                                                              gTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTcc                                                                                                                     >  2:525338/1‑251 (MQ=255)
                                                                                                                    aTTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGatat                                                                                                               >  2:473191/1‑251 (MQ=255)
                                                                                                                        a‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAc                                                                                                           <  1:409543/251‑1 (MQ=255)
                                                                                                                           ccTTGTTTCCACAAAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCAGCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGTAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGAccc                                                                                                         <  2:623051/251‑1 (MQ=255)
                                                                                                                              tgtttttaaaaagtgtggccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAAGCATCAGTAAAATCTATTCATTATCACAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAGCGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGGGAGTTGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCAt                                                                                                       <  2:236404/232‑1 (MQ=255)
                                                                                                                                        aaaTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATGAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATGACGGGCAATg                                                                                            <  2:595073/250‑1 (MQ=255)
                                                                                                                                          aTCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGAt                                                                                          >  1:537329/1‑251 (MQ=255)
                                                                                                                                              gCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATaaaa                                                                                      <  1:558458/251‑1 (MQ=255)
                                                                                                                                                ccgccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTAATCATTATCTCAATCAGGCCGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAgg                                                                                    <  1:229716/251‑1 (MQ=255)
                                                                                                                                                   ccgATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAACCGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGt                                                                                 <  1:473191/251‑1 (MQ=255)
                                                                                                                                                     gATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTa                                                                                >  1:581382/1‑250 (MQ=255)
                                                                                                                                                      aTTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATaacaac                                                                                                                                             >  1:284484/1‑188 (MQ=255)
                                                                                                                                                      aTTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAAGaacaac                                                                                                                                             <  2:284484/188‑1 (MQ=255)
                                                                                                                                                        tttGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑t                                                                           >  1:487551/1‑251 (MQ=255)
                                                                                                                                                          tgctgcGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑tat                                                                         >  1:329579/1‑249 (MQ=255)
                                                                                                                                                                   gcgtAAATTGAGGAGGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAATAACGACAGGGAAAAAGGAGAAAGTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTGCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGAt                                                                <  1:36974/251‑1 (MQ=255)
                                                                                                                                                                     gtAAATTGATGATGAATAATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGc                                                              >  2:87016/1‑251 (MQ=255)
                                                                                                                                                                       aaaTTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCa                                                             >  1:346485/1‑250 (MQ=255)
                                                                                                                                                                            gatgatgaATCATCAGTAAAAACTATTCATTATCTCAATCAGGCCGGGTTTGCTTGTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCAACGCGTAATCTCCGGATAACGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATcga                                                        <  2:525317/250‑3 (MQ=255)
                                                                                                                                                                               gatgaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCg                                                     <  1:400548/250‑1 (MQ=255)
                                                                                                                                                                                  gaATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTAc                                                  <  1:164638/250‑1 (MQ=255)
                                                                                                                                                                                      catcaGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCa                                             <  2:166611/251‑1 (MQ=255)
                                                                                                                                                                                         caGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACtt                                          <  2:2479/251‑1 (MQ=255)
                                                                                                                                                                                               actCTAGTCATTATCTCACTCAGGCCGGGTTTACTTATATGCAGCC‑GGCTTTTTTTTCAAGAAATTATGGAGAAAAACGACAGGGAAATAGGTGAAATTCTCAATAAATGCGGTAACTAAGAGATTAGGATTGCGGAGAATAACAACCGCCGTACTCATCGCGTAATCTCCGGATATAAACCCATAACGGGCAATGCTAACAGGA‑TCACCATATGAAAAAGCTGAAATCTATCGTACTCGAACTTTCCCTg                                    >  2:109328/3‑251 (MQ=255)
                                                                                                                                                                                                aTCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTgg                                   >  2:178695/1‑251 (MQ=255)
                                                                                                                                                                                                      tCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTgg                             >  1:411235/1‑251 (MQ=255)
                                                                                                                                                                                                        aTTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTGTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTc                           <  2:212253/251‑1 (MQ=255)
                                                                                                                                                                                                                   cAGGCCGGGTTTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGgc                >  2:358045/1‑251 (MQ=255)
                                                                                                                                                                                                                       ccGGGTGTGCTTTTATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGATAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGAGTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAATATGCAATCTATCGTACTCGCACTTTCCCTCGTTCTGGCCGCTCCCATGgcagca            <  2:207439/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ttttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:448274/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tttATGCAGCC‑GGCTTTTTTATGAAGAAATTATGGAGAAAAACGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGTCGTTCTCATCGCGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TATGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCgg  >  1:329905/1‑250 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                
CCAGGCGAATATGGCTTGTTCCTCGGCACCGCGCATCCGGCGAAATTTAAAGAGAGCG‑TGGAAGCGATTCTCGGTGAAACGTTGGATCTGCCAAAAGAGCTGGCAGAACGTGCTGATTTA‑CCCTTGCTTTCACATAATCTGCCCGCCGATTTTGCTGCGTTGCGTAAATTGATGATGAATCATCAGTAAAATCTATTCATTATCTCAATCAGGCCGGGTTTGCTTTTATGCAGCCCGGCTTTTTTATGAAGAAATTATGGAGAAAAATGACAGGGAAAAAGGAGAAATTCTCAATAAATGCGGTAACTTAGAGATTAGGATTGCGGAGAATAACAACCGCCGTTCTCATCGAGTAATCTCCGGATATCGACCCATAACGGGCAATGATAAAAGGAGTAACC‑TGTGAAAAAGATGCAATCTATCGTACTCGCACTTTCCCTGGTTCTGGTCGCTCCCATGGCAGCACAGGCTGCGG  >  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: ‑