Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,121,762 G→A V282V (GTC→GTT)  pyrC ← dihydroorotase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,121,7620GA98.0% 142.9 / NA 50V282V (GTC→GTT) pyrCdihydroorotase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (33/15);  total (34/16)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGC  >  NC_000913/1121519‑1122002
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ataCAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGGTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑g                                                                                                                                                                                                                                           <  1:643055/251‑1 (MQ=255)
    aGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAAGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGt                                                                                                                                                                                                                                        >  1:94449/1‑250 (MQ=255)
        aaaaTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTg                                                                                                                                                                                                                                   <  1:622525/251‑1 (MQ=255)
            tAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCACCCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAAACCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGATCTGCAAAGCATTCATTTCTTCAAAAAAGGTA‑GCGTAACTGccac                                                                                                                                                                                                                               >  2:531360/1‑249 (MQ=255)
               ccTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCAAgcg                                                                                                                                                                                                                            >  1:584921/1‑251 (MQ=255)
                       caacaaGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTggg                                                                                                                                                                                                                    >  1:379574/1‑251 (MQ=255)
                              aaGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTg                                                                                                                                                                                                              >  1:589501/1‑250 (MQ=255)
                                      aTTTTTATTGTTTAACGGACCATCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGACAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGGA‑GCGTAACAGCCCAGCGCGGTTGGGGCGTTGAagcagcac                                                                                                                                                                                                     >  2:357537/1‑249 (MQ=255)
                                      aTTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGccc                                                                                                                                                                                                     >  1:50139/1‑251 (MQ=255)
                                         tttatagcttAACGGACCAGCGTACCGTTTCACCGCCGAGGACTGGAACCAGCGTGTCATCAGTCAGTGCGATGCTTTCCGCAACCTGTTGCTCTTCAACTACCCGTTCGATGAATGTGTCGTTGACCGGAAACCCATCGAAATGCGGGCCGTTTACCGAAAAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑ACGTAACTGCCCAGCGCGGTTGGGGCGTTGAagaagaacgac                                                                                                                                                                                                  >  2:567311/9‑242 (MQ=255)
                                          ttATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCAATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCcgagc                                                                                                                                                                                                 >  2:663797/1‑248 (MQ=255)
                                              tGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGcc                                                                                                                                                                                             >  1:183928/1‑251 (MQ=255)
                                              tGTTTAACGCACCAGCGGACCGTTGCCCCGGCGAGGAAGGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGc                                                                                                                                                                                              <  2:589501/250‑1 (MQ=255)
                                                tttAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGc                                                                                                                                                                                           >  1:82013/1‑251 (MQ=255)
                                                       gACCAGCGTAACATTTCCCCGGCGAGGCCTGCAACCAGCGGGTCATCGGCCCGTGAGATGCATTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGGCGTTGACCGGCAACCCATAGAACCGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGGACAGACGCAACTGcc                                                                                                                                                                                    >  2:351056/1‑250 (MQ=255)
                                                                   gTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATg                                                                                                                                                                        >  1:418712/1‑251 (MQ=255)
                                                                      tCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATgtcg                                                                                                                                                                     >  1:97731/1‑248 (MQ=255)
                                                                          cGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAA‑TACCGTAAGCGTAACTGCCAAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCAGCAACAGCTCTCTTTGCGGTGACGTGc                                                                                                                                                                  >  2:130753/1‑250 (MQ=255)
                                                                            gCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCAATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAtg                                                                                                                                                               <  1:663797/251‑1 (MQ=255)
                                                                                ggAATGGCACCAGCGTTTCATCAGTCAGTTCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgc                                                                                                                                                           <  2:359931/251‑1 (MQ=255)
                                                                                ggAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgc                                                                                                                                                           <  1:87674/251‑1 (MQ=255)
                                                                                 gAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgcg                                                                                                                                                          >  1:551934/1‑251 (MQ=255)
                                                                                   aTGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgcgc                                                                                                                                                         >  1:414271/1‑250 (MQ=255)
                                                                                                      agtcagtGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGCTGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAAt                                                                                                                                      >  1:582613/1‑250 (MQ=255)
                                                                                                                                    ctctTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCAc                                                                                                       >  1:595265/1‑251 (MQ=255)
                                                                                                                                       ttCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                                                                                     >  1:272353/1‑250 (MQ=255)
                                                                                                                                          ccttACCAGTTCGATGAAGGTGCCGTTGCCCGGCAACTCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGGGCTGCAAAGCATTCTTTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGGGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATg                                                                                                 <  2:379574/248‑1 (MQ=255)
                                                                                                                                          aCGTACAAGTTCGATGAATGTGTCGTTGACAGGCAACACATAGAACTGCGGGCCGTTTACAGAACAGAATACATAAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGGA‑GCGTAACTACCCAGCGCGGTTGGGGCGTTGAAGAAGCCCGCGCAGCCGCACCTGCTCTCTTTTCGATGCCGTGCATGTGGCGCCGAAACCGTACCGAGGAATACTTGATTAATACCGCGAGCTACCAGTTCACTCAAtt                                                                                                 >  2:11778/1‑250 (MQ=255)
                                                                                                                                                 agttcaatgatgcggtCGTTGACCGGCATCCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGAGCTGCAAAGCATTAATGTCTTCAAAAACGGTA‑GCGTAACTGCCCTGCACGGTTGGGGCTTTGAAGCAGACCGTGCAGCCGCAACGGCTCTATTGGCGATAACGTGCATGTGGCGCAGAATCCGTACCGAGGAATTCTCGATTAAAACCGCTGGCGACCAGTTCACGCACTGCCTGTTg                                                                                          <  1:387738/237‑1 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAt                                                                                                                                                                <  1:537096/168‑1 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAt                                                                                                                                                                >  2:537096/1‑168 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGa                                                                                                                                >  1:420415/1‑200 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGGAACTGCCCAGCGCGGTTTGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTTCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGa                                                                                                                                <  2:420415/200‑1 (MQ=255)
                                                                                                                                                                 cGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTTCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCTCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTg                                                                          <  2:677677/251‑1 (MQ=255)
                                                                                                                                                                       ccGGCACCCCATAGAACTGCAGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAACGCGGTTGGGGCGTTGAAGAAACACGCGCAGCCGCAACTGATCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGCGGAATACTCGATTAACACCGCTGGCGACCAATTCCCGCAATGCCTGTTGGTGACTATTACGTTTGCGGATg                                                                    >  2:632361/1‑251 (MQ=255)
                                                                                                                                                                          gCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCCTTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGAAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATAATCGATTAAAACCGCGGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTggggtggggt                                                                 >  2:496301/1‑242 (MQ=255)
                                                                                                                                                                                     aaCTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGAAGACAATAc                                                        >  1:673265/1‑249 (MQ=255)
                                                                                                                                                                                              ccGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGGTTGAGGATGGGTAGACAATACAGGTGCGGacg                                             <  1:351056/251‑1 (MQ=255)
                                                                                                                                                                                                  ttACAGAACAGAATGCTTCTAAGGGCTGCAAAGCATTCATTTCTACAAAACCGGTA‑GCGATACTGCCGAGCGCGGTGGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCCCAGAATCCGTACCGAGGAATTCTCGATTAAATCCGCTGGCGACCAGTTCACGCAATGCCGGTTGGTGAATATTACGTTTGAGGATGCGTTGACAATACATGTGCGGacgctc                                          <  2:8717/250‑3 (MQ=255)
                                                                                                                                                                                                        aacagaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGCTTAAAACCGCTGGCGCCCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGACTcc                                    >  1:336217/1‑250 (MQ=255)
                                                                                                                                                                                                           agaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTa                                                                                                                             <  2:504389/158‑1 (MQ=255)
                                                                                                                                                                                                           agaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTa                                                                                                                             >  1:504389/1‑158 (MQ=255)
                                                                                                                                                                                                           agaaTGATTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTGTGCGATGACGTGCATGTGGCGCAGAAGCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGCCCAGTTCACGCGATGCCTGTTGGTGCATATTAAGTTTGAGGATGGGTAGACAATACAGGTGCGGTAGCAAGCCTCCAAcc                                >  2:516868/1‑251 (MQ=255)
                                                                                                                                                                                                              aTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAACAGCCCGCGCAGCCGCACCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGc                             >  2:154050/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttCAAAATGCTGCAAAGCATTCATTTCTTCAACAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGCATACTCGATTAAAACCGCTGGCGACCAGATCACGAAATGCCTGTTGTGGAATATTACTGTTGAGGATGGGTAGACAATACAGGTGCGAGCGAACGCCTCCAACCAGAatat                         >  2:514921/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGAGAATATTACGTTTGCGGATGGGTAGACAATACAGGTGCGGAAGCACACCTCCAACCAGCATATGGtt                     >  2:11688/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGTGGCGTTGAAGCAGCCCGCGCAGCAGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACAGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTg                    >  2:287749/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaGCATTCATTTCTTCAAAAACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCATGTTGGGGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGAAGCACGCCTCCAACCAGCATATGGTTGCGGTTAAAc           >  2:458038/1‑251 (MQ=255)
                                                                                                                                                                                                                                        attgttttataGACGGTA‑GAGTAACTGTCAAGTGTGGCTGGGTCGTTGAAGAAGCACGTGAAGCCGCAACTGCTCTCGTTGCGATGAAGTGAATGTGGCGCAGAATTCGTACCGAGGAATACTCGTTTAAATCCGCTGGCGACCAGTTCACGTAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGATAATACAGGTGCGGACGAATGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtg   <  1:150368/241‑1 (MQ=255)
                                                                                                                                                                                                                                          ttcttcAAAACCGGTC‑GTGTACCTGCCGAGCGCTGTTGGGGCGTTGAAGCAGCCCGAGCATCCGCAACTGCTCTCTTTTCGATGACGTGCATGGTGCGCATAATCCGTTCCGAGGAATCCCCGATTAAAACCGCTGGCGACCAGTTCTCGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCTTATGGTTGCGTTTAAACATCAGAtgc  <  2:582613/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTA‑GCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGC  >  NC_000913/1121519‑1122002

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