Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,410,666 C→G V129V (GTC→GTG dusB → tRNA‑dihydrouridine synthase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,410,6660CG98.6% 249.3 / NA 73V129V (GTC→GTGdusBtRNA‑dihydrouridine synthase B
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (35/37);  minor base A (1/0);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 4.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATCGCGAATGGCGACATTACTGACCCGCTTAAAGCC  >  NC_000913/3410426‑3410906
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
ggAAAGCGGAAAATCTCGTTTTCCGATGGTGNACATTGATGAACCCGGTATTCGCACCGTGCCAATTGCTGGTAGGGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAAATTATTGAAATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTgg                                                                                                                                                                                                                                        <  2:235497/251‑1 (MQ=255)
  aaaGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTGACCGAGGTCGTGAATGCAGTGGAc                                                                                                                                                                                                                                      <  2:173987/251‑1 (MQ=255)
    aGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGt                                                                                                                                                                                                                                    <  1:363288/251‑1 (MQ=255)
      cGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTc                                                                                                                                                                                                                                  <  1:7446/251‑1 (MQ=255)
              ctcGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGCGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTAcc                                                                                                                                                                                                                          >  1:227261/1‑251 (MQ=255)
                cGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGt                                                                                                                                                                                                                                          <  1:52497/233‑1 (MQ=255)
                cGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGt                                                                                                                                                                                                                                          >  2:52497/1‑233 (MQ=255)
                cGTTTACGGTGTGTGTACATTGATGATCNCGGTATTCGCGACGTGTAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGCGGAAAGCGGAGCCCAGATTATTGATAGCAATATTGGTTGCCCGGCTAAAAAAGGTAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGTTGTCGTTAAATCGATCCTTACCGAGGGCGTGAATGCAATGGACGTTCCTGTTAc                                                                                                                                                                                                                           <  2:269285/248‑1 (MQ=255)
                       ggATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGAtt                                                                                                                                                                                                                 >  1:542836/1‑251 (MQ=255)
                             gTGCACATTGATGAACCCGGTATTAGCACCGTGCAAAGTGCTGGTAGCGATCCGAAAGCAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGCTTATTGATATCAATATGGGTTGCCGGGCTAAAAAAGGGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGACGTCGTTAAATCGAACCTTACCGAGGGCGTGAATGCAGTGGACGTACCGGTAACCCTGAAGATTGGCAc                                                                                                                                                                                                            >  2:117714/1‑251 (MQ=255)
                                cacaTTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCgg                                                                                                                                                                                                         <  2:373025/251‑1 (MQ=255)
                                     tgatgaACCCGGTAATCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGAGTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTGACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTgt                                                                                                                                                                                                     <  2:371373/250‑2 (MQ=255)
                                         gaACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGCGAATCGCAAGCTCGCCGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGGGGTCCTGAATGCAGTGGACGTTCCTGGTTCCCTGAAGATTCGCACCGGCTGGGCAcc                                                                                                                                                                                                >  2:114125/1‑251 (MQ=255)
                                         gaACCCGGTAGTCGCACCGTGCAAATTGCTGNTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTTCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTGGCCCTGAAGATTCGCACCGGCTGGGCAcc                                                                                                                                                                                                <  2:245129/251‑1 (MQ=255)
                                                gTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGAAGCAGCACGTATTAACGTGGAAAGCAGCGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAACAAAATGAATCGCAAGCTCGCAGGTTCAGCCCTCCTGCAGTACCCGGAGGACGTTAAAACGTGCCTTACCGAGGGCGTGAATGCAGAGGAAGTTCCTGTTACACTGCAGATTCGGCCCGGCTGGGCCCCGGaaaa                                                                                                                                                                                          >  2:316308/1‑248 (MQ=255)
                                                  aTTCGCTCCTTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACTTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCg                                                                                                                                                                                       <  2:433784/251‑1 (MQ=255)
                                                    tCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGt                                                                                                                                                                                      >  1:103919/1‑250 (MQ=255)
                                                             gCAAATTGCTGGTAGCGATCCGAAAGAAATGGCCGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGAGTATTGATATCAATATGGGGGGCCCGGCTAACAAAGCGAAGCGCAAGCACGCAGGGGCAGCCCTCTTGCAGTACCCGAATGTCGATAAATCGATCCTTACCGAGGTCGTACATGCAGTGGACGACCCTGTTCCCCTGAAAATTCGCAACGACTGGGCACCggtgaaacgaaaatgcga                                                                                                                                                                              >  2:102549/1‑233 (MQ=255)
                                                              cAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACTCCGTAACTGCGAaga                                                                                                                                                                           <  1:144178/251‑1 (MQ=255)
                                                               aaaTTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAagag                                                                                                                                                                          <  2:207453/251‑1 (MQ=255)
                                                                   tGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGAtt                                                                                                                                                                       <  1:118592/250‑1 (MQ=255)
                                                                    gCTGGGAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCTCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTg                                                                                                                                                                      <  2:212338/250‑1 (MQ=255)
                                                                        gTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGTAGAGATTGCCCa                                                                                                                                                                  >  2:25886/1‑250 (MQ=255)
                                                                         tAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGATGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAACAGATTGCCCAAc                                                                                                                                                                >  1:301304/1‑251 (MQ=255)
                                                                          aGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGTTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCTCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACt                                                                                                                                                               <  1:207438/251‑1 (MQ=255)
                                                                            cGATCCGAAAGAAATGGCAGATGCAGCACGTAGTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTgg                                                                                                                                                             <  1:216282/251‑1 (MQ=255)
                                                                             gATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATCTGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGCGATTGCCCAACTGGc                                                                                                                                                            >  1:182482/1‑251 (MQ=255)
                                                                                  gaaagaaaTGGCAGATGCAGCACGGATTAACGTGGAAAGCGGTGACCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGAAAATGCGCAGAGATTGCCCAACTGGCTgc                                                                                                                                                         >  2:5468/1‑248 (MQ=255)
                                                                                       aaaTGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTTCCCAAATGGCTGAAGACtgg                                                                                                                                                  >  2:363495/1‑250 (MQ=255)
                                                                                        aaTGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACtgtg                                                                                                                                                 <  1:361604/251‑1 (MQ=255)
                                                                                             aaGATGCAGCACGTATGAACGTGGAAAGCGGTGCCCAGATTATTGAGATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACGCCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCAtt                                                                                                                                            <  2:110311/250‑1 (MQ=255)
                                                                                                    agcaCGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGctc                                                                                                                                     <  1:499455/251‑1 (MQ=255)
                                                                                                         gTATTAACGTGGAAAGCGGTGCCCAGATTATTGCTATCAATATGGGTTGCCCGGCTAAAAAAGCGAATCGCAAGCTAGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCCCCGGAACACCGTAACTGCGAAGAGACTGCCCAACTGGCTGAAGCCTGTGGCATTCAGGCTCTGAcc                                                                                                                                >  2:22053/1‑251 (MQ=255)
                                                                                                         gTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGAcc                                                                                                                                <  1:89190/251‑1 (MQ=255)
                                                                                                           aTTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCAt                                                                                                                              >  1:473971/1‑251 (MQ=255)
                                                                                                               aCGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCTGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCTGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCAt                                                                                                                          <  1:65908/251‑1 (MQ=255)
                                                                                                                         cGGTGCCCAGATTCTTGATATCAATATGGGTTGCCCGGCTAAAAAAGCGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTacac                                                                                                                >  2:66055/1‑251 (MQ=255)
                                                                                                                          ggTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAcg                                                                                                               <  2:295320/251‑1 (MQ=255)
                                                                                                                                 aGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGACGAGATTGCCCAACTGGCTGAAGACTGGGGCATTCAGGCTCTGACCATTCATGGCCGTACAAGCGCCtgg                                                                                                        >  2:317007/1‑250 (MQ=255)
                                                                                                                                 aGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCtgt                                                                                                        >  1:536609/1‑251 (MQ=255)
                                                                                                                                  gATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGTCCATTCATGGCCGTACACGCGCCtgtt                                                                                                       <  1:415145/251‑1 (MQ=255)
                                                                                                                                         gATATCAATATGGGTTGCCCTGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCa                                                                                                 <  2:27478/250‑1 (MQ=255)
                                                                                                                                          atatCAATATGGGTTGCCCGGATAAAAAAGTGAAACGCAAGCTCGCAGGTTCAGCCCTCGTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAAt                                                                                               <  2:143671/251‑1 (MQ=255)
                                                                                                                                             tCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATgg                                                                                             >  2:168412/1‑250 (MQ=255)
                                                                                                                                              cAATAGGGGTTGCCAGGCTAAAAAAGGGAATCGCAAGCTCGCAGGTTCAGCCATCTTGCAGTACCCGGATGTCGTTAAATCGATCCGTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGTCTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGga                                                                                            <  2:81111/250‑1 (MQ=255)
                                                                                                                                                    gggTTGCCCGGCTAAAAAAGTGAATCTCAAGCTCGCATGTTCCGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACGGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                                                                                     <  2:336016/251‑1 (MQ=255)
                                                                                                                                                      gTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGCTTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTAAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTg                                                                                     >  1:284238/1‑249 (MQ=255)
                                                                                                                                                         gCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTa                                                                                 <  1:384718/250‑1 (MQ=255)
                                                                                                                                                         gCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTAc                                                                                <  1:12868/251‑1 (MQ=255)
                                                                                                                                                          cccGGCTAACAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGGCGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGAATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTATGTTCAATGGAGAAGCTGAGttcg                                                                               >  2:498689/1‑248 (MQ=255)
                                                                                                                                                          cccGGCTAAAAAAGGGAATCGCAAGCTCGCANGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACg                                                                               <  2:270007/251‑1 (MQ=255)
                                                                                                                                                           ccGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGa                                                                              <  1:114125/251‑1 (MQ=255)
                                                                                                                                                            cGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAACCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCCTTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGAc                                                                             >  2:432566/1‑251 (MQ=255)
                                                                                                                                                                 aaaaaaGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGGCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTATCTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTAt                                                                        <  2:542836/251‑1 (MQ=255)
                                                                                                                                                                   aaaaGTGAATCGCCAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTc                                                                      <  2:311705/251‑1 (MQ=255)
                                                                                                                                                                      aGTGAATCGCAAGCTCGCAGGTTCCGCCCTCTTGCCGTACCCGGATGTCGTTAAAGCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGTCACCGGAAAACCGTAACTGCGAAGAGATTGCCCAACTGGCGGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATTGAGAAGCTGAGTACGACAGTATTCggg                                                                   <  2:540156/251‑1 (MQ=255)
                                                                                                                                                                       gTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCACTGGAGAAGCTGAGTACGACAGTATTCgggg                                                                  >  2:103110/1‑250 (MQ=255)
                                                                                                                                                                          aaTCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTCCACGCGCCTGTTTGTTCAATGGAGCACCTGAGTACGACAGTATTCGGGccgt                                                               >  2:522493/1‑248 (MQ=255)
                                                                                                                                                                              gCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGGGAAGCTGAGTACGCCAGGATTCGGGCAGTTa                                                             >  2:119627/1‑249 (MQ=255)
                                                                                                                                                                               cAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAGGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGCCAGTATTCGGGCAGTTAAGc                                                          >  1:18159/1‑251 (MQ=255)
                                                                                                                                                                                aatctctcAGGTTCCGCCCTCGTGCAGTACCCGGATGCCGTTAAATCGATCCTTACCGGGGTCGTGAATGCAGTGGACGGTCCTGTTACCCTGAAGATTCTCACCGGATTGTCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTTGCTGGAGACTTTGGCATTCAGGCTCTTACCATTCATGGCCGTACACTCGCCTGTTTGTTCAATGTAGAAGATGACTTCGTCAGTATCCGGGCAGTTaatca                                                         <  2:324727/244‑4 (MQ=255)
                                                                                                                                                                                aaGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACAAGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGc                                                          >  1:56202/1‑250 (MQ=255)
                                                                                                                                                                                        aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGCGGTCTGGAAAGCAGTGGACGTTCCTGTTGCCCTGAAGATTCGCACCGGCTGGGCACCGGAACCACGTAACTGCGACGCGATTTCCCAACTGGCTGAAGACTGTGGCATTCGGGGTCTGACCATTCATGGCCGGACGCGCGCCTGTTTGTTCAGTTGCGGAGCGGGATAAGACAGAATTAGGGCAGTTAAGCAGAAAGttt                                                 >  2:393061/1‑251 (MQ=255)
                                                                                                                                                                                        aGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGttt                                                 >  1:244035/1‑251 (MQ=255)
                                                                                                                                                                                             cAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCAtt                                            <  1:328673/251‑1 (MQ=255)
                                                                                                                                                                                              aGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTc                                           >  1:490972/1‑251 (MQ=255)
                                                                                                                                                                                                   tctTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGtt                                      >  1:60805/1‑251 (MQ=255)
                                                                                                                                                                                                         aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCCTTCCGGTTATcgcg                                >  1:284470/1‑251 (MQ=255)
                                                                                                                                                                                                         aGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATcgcg                                >  1:115707/1‑251 (MQ=255)
                                                                                                                                                                                                          gTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTGCATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGCCAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATctcg                                >  2:308242/1‑247 (MQ=255)
                                                                                                                                                                                                                    gTCGTTAAATCGATCCTTACCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATCGCGAATGGCGACAt                     <  1:310693/251‑1 (MQ=255)
                                                                                                                                                                                                                                      aCCGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATCGCGAATGGCGACATTACTGACCCGCTTAAAGc   >  1:210723/1‑251 (MQ=255)
                                                                                                                                                                                                                                        cGAGGTCGTGAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATCGCGAATGGCGACATTACTGACCCGCTTAAAGcc  >  1:518802/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                
GGAAAGCGACAAATCTCGTTTACGGATGGTGCACATTGATGAACCCGGTATTCGCACCGTGCAAATTGCTGGTAGCGATCCGAAAGAAATGGCAGATGCAGCACGTATTAACGTGGAAAGCGGTGCCCAGATTATTGATATCAATATGGGTTGCCCGGCTAAAAAAGTGAATCGCAAGCTCGCAGGTTCAGCCCTCTTGCAGTACCCGGATGTCGTTAAATCGATCCTTACCGAGGTCGTCAATGCAGTGGACGTTCCTGTTACCCTGAAGATTCGCACCGGCTGGGCACCGGAACACCGTAACTGCGAAGAGATTGCCCAACTGGCTGAAGACTGTGGCATTCAGGCTCTGACCATTCATGGCCGTACACGCGCCTGTTTGTTCAATGGAGAAGCTGAGTACGACAGTATTCGGGCAGTTAAGCAGAAAGTTTCCATTCCGGTTATCGCGAATGGCGACATTACTGACCCGCTTAAAGCC  >  NC_000913/3410426‑3410906

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: