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→GTTpyrCdihydroorotase
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

ATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGC  >  NC_000913/1121519‑1122002
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ataCAGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGGTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAg                                                                                                                                                                                                                                           <  1:643055/251‑1 (MQ=255)
    aGGTAAAATAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAAAGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGt                                                                                                                                                                                                                                        >  1:94449/1‑250 (MQ=255)
        aaaaTAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTg                                                                                                                                                                                                                                   <  1:622525/251‑1 (MQ=255)
            tAACCTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCACCCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAAACCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGATCTGCAAAGCATTCATTTCTTCAAAAAAGGTAGCGTAACTGccac                                                                                                                                                                                                                               >  2:531360/1‑249 (MQ=255)
               ccTAATGACAACAAGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCAAgcg                                                                                                                                                                                                                            >  1:584921/1‑251 (MQ=255)
                       caacaaGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTggg                                                                                                                                                                                                                    >  1:379574/1‑251 (MQ=255)
                              aaGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTg                                                                                                                                                                                                              >  1:589501/1‑250 (MQ=255)
                                      aTTTTTATTGTTTAACGGACCATCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGACAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGGAGCGTAACAGCCCAGCGCGGTTGGGGCGTTGAagcagcac                                                                                                                                                                                                     >  2:357537/1‑249 (MQ=255)
                                      aTTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGccc                                                                                                                                                                                                     >  1:50139/1‑251 (MQ=255)
                                         tttatagcttAACGGACCAGCGTACCGTTTCACCGCCGAGGACTGGAACCAGCGTGTCATCAGTCAGTGCGATGCTTTCCGCAACCTGTTGCTCTTCAACTACCCGTTCGATGAATGTGTCGTTGACCGGAAACCCATCGAAATGCGGGCCGTTTACCGAAAAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAACGTAACTGCCCAGCGCGGTTGGGGCGTTGAagaagaacgac                                                                                                                                                                                                  >  2:567311/9‑242 (MQ=255)
                                          ttATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCAATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCcgagc                                                                                                                                                                                                 >  2:663797/1‑248 (MQ=255)
                                              tGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGcc                                                                                                                                                                                             >  1:183928/1‑251 (MQ=255)
                                              tGTTTAACGCACCAGCGGACCGTTGCCCCGGCGAGGAAGGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGc                                                                                                                                                                                              <  2:589501/250‑1 (MQ=255)
                                                tttAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGc                                                                                                                                                                                           >  1:82013/1‑251 (MQ=255)
                                                       gACCAGCGTAACATTTCCCCGGCGAGGCCTGCAACCAGCGGGTCATCGGCCCGTGAGATGCATTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGGCGTTGACCGGCAACCCATAGAACCGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGGACAGACGCAACTGcc                                                                                                                                                                                    >  2:351056/1‑250 (MQ=255)
                                                                   gTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATg                                                                                                                                                                        >  1:418712/1‑251 (MQ=255)
                                                                      tCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATgtcg                                                                                                                                                                     >  1:97731/1‑248 (MQ=255)
                                                                          cGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAATACCGTAAGCGTAACTGCCAAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCAGCAACAGCTCTCTTTGCGGTGACGTGc                                                                                                                                                                  >  2:130753/1‑250 (MQ=255)
                                                                            gCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCAATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAtg                                                                                                                                                               <  1:663797/251‑1 (MQ=255)
                                                                                ggAATGGCACCAGCGTTTCATCAGTCAGTTCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgc                                                                                                                                                           <  2:359931/251‑1 (MQ=255)
                                                                                ggAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgc                                                                                                                                                           <  1:87674/251‑1 (MQ=255)
                                                                                 gAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgcg                                                                                                                                                          >  1:551934/1‑251 (MQ=255)
                                                                                   aTGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGgcgc                                                                                                                                                         >  1:414271/1‑250 (MQ=255)
                                                                                                      agtcagtGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGCTGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAAt                                                                                                                                      >  1:582613/1‑250 (MQ=255)
                                                                                                                                    ctctTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCAc                                                                                                       >  1:595265/1‑251 (MQ=255)
                                                                                                                                       ttCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGc                                                                                                     >  1:272353/1‑250 (MQ=255)
                                                                                                                                          ccttACCAGTTCGATGAAGGTGCCGTTGCCCGGCAACTCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGGGCTGCAAAGCATTCTTTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGGGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATg                                                                                                 <  2:379574/248‑1 (MQ=255)
                                                                                                                                          aCGTACAAGTTCGATGAATGTGTCGTTGACAGGCAACACATAGAACTGCGGGCCGTTTACAGAACAGAATACATAAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGGAGCGTAACTACCCAGCGCGGTTGGGGCGTTGAAGAAGCCCGCGCAGCCGCACCTGCTCTCTTTTCGATGCCGTGCATGTGGCGCCGAAACCGTACCGAGGAATACTTGATTAATACCGCGAGCTACCAGTTCACTCAAtt                                                                                                 >  2:11778/1‑250 (MQ=255)
                                                                                                                                                 agttcaatgatgcggtCGTTGACCGGCATCCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGAGCTGCAAAGCATTAATGTCTTCAAAAACGGTAGCGTAACTGCCCTGCACGGTTGGGGCTTTGAAGCAGACCGTGCAGCCGCAACGGCTCTATTGGCGATAACGTGCATGTGGCGCAGAATCCGTACCGAGGAATTCTCGATTAAAACCGCTGGCGACCAGTTCACGCACTGCCTGTTg                                                                                          <  1:387738/237‑1 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAt                                                                                                                                                                <  1:537096/168‑1 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCAt                                                                                                                                                                >  2:537096/1‑168 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGa                                                                                                                                >  1:420415/1‑200 (MQ=255)
                                                                                                                                                              tgtCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGGAACTGCCCAGCGCGGTTTGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTTCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGa                                                                                                                                <  2:420415/200‑1 (MQ=255)
                                                                                                                                                                 cGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTTCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCTCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTg                                                                          <  2:677677/251‑1 (MQ=255)
                                                                                                                                                                       ccGGCACCCCATAGAACTGCAGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAACGCGGTTGGGGCGTTGAAGAAACACGCGCAGCCGCAACTGATCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGCGGAATACTCGATTAACACCGCTGGCGACCAATTCCCGCAATGCCTGTTGGTGACTATTACGTTTGCGGATg                                                                    >  2:632361/1‑251 (MQ=255)
                                                                                                                                                                          gCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCCTTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGAAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATAATCGATTAAAACCGCGGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTggggtggggt                                                                 >  2:496301/1‑242 (MQ=255)
                                                                                                                                                                                     aaCTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGAAGACAATAc                                                        >  1:673265/1‑249 (MQ=255)
                                                                                                                                                                                              ccGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGGTTGAGGATGGGTAGACAATACAGGTGCGGacg                                             <  1:351056/251‑1 (MQ=255)
                                                                                                                                                                                                  ttACAGAACAGAATGCTTCTAAGGGCTGCAAAGCATTCATTTCTACAAAACCGGTAGCGATACTGCCGAGCGCGGTGGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCCCAGAATCCGTACCGAGGAATTCTCGATTAAATCCGCTGGCGACCAGTTCACGCAATGCCGGTTGGTGAATATTACGTTTGAGGATGCGTTGACAATACATGTGCGGacgctc                                          <  2:8717/250‑3 (MQ=255)
                                                                                                                                                                                                        aacagaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGCTTAAAACCGCTGGCGCCCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGACTcc                                    >  1:336217/1‑250 (MQ=255)
                                                                                                                                                                                                           agaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTa                                                                                                                             <  2:504389/158‑1 (MQ=255)
                                                                                                                                                                                                           agaaTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTa                                                                                                                             >  1:504389/1‑158 (MQ=255)
                                                                                                                                                                                                           agaaTGATTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTGTGCGATGACGTGCATGTGGCGCAGAAGCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGCCCAGTTCACGCGATGCCTGTTGGTGCATATTAAGTTTGAGGATGGGTAGACAATACAGGTGCGGTAGCAAGCCTCCAAcc                                >  2:516868/1‑251 (MQ=255)
                                                                                                                                                                                                              aTGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAACAGCCCGCGCAGCCGCACCTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGc                             >  2:154050/1‑251 (MQ=255)
                                                                                                                                                                                                                  ttCAAAATGCTGCAAAGCATTCATTTCTTCAACAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGCATACTCGATTAAAACCGCTGGCGACCAGATCACGAAATGCCTGTTGTGGAATATTACTGTTGAGGATGGGTAGACAATACAGGTGCGAGCGAACGCCTCCAACCAGAatat                         >  2:514921/1‑251 (MQ=255)
                                                                                                                                                                                                                      aaGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGAGAATATTACGTTTGCGGATGGGTAGACAATACAGGTGCGGAAGCACACCTCCAACCAGCATATGGtt                     >  2:11688/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGTGCTGCAAAGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGTGGCGTTGAAGCAGCCCGCGCAGCAGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACAGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTg                    >  2:287749/1‑251 (MQ=255)
                                                                                                                                                                                                                                aaGCATTCATTTCTTCAAAAACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCATGTTGGGGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGAAGCACGCCTCCAACCAGCATATGGTTGCGGTTAAAc           >  2:458038/1‑251 (MQ=255)
                                                                                                                                                                                                                                        attgttttataGACGGTAGAGTAACTGTCAAGTGTGGCTGGGTCGTTGAAGAAGCACGTGAAGCCGCAACTGCTCTCGTTGCGATGAAGTGAATGTGGCGCAGAATTCGTACCGAGGAATACTCGTTTAAATCCGCTGGCGACCAGTTCACGTAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGATAATACAGGTGCGGACGAATGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGAtg   <  1:150368/241‑1 (MQ=255)
                                                                                                                                                                                                                                          ttcttcAAAACCGGTCGTGTACCTGCCGAGCGCTGTTGGGGCGTTGAAGCAGCCCGAGCATCCGCAACTGCTCTCTTTTCGATGACGTGCATGGTGCGCATAATCCGTTCCGAGGAATCCCCGATTAAAACCGCTGGCGACCAGTTCTCGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCTTATGGTTGCGTTTAAACATCAGAtgc  <  2:582613/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                 
ATACAGGTAAAATAACCTAATGACAACAGGAAGCTACGATTTTTATTGTTTAACGGACCAGCGTACCGTTTCCCCGGCGAGGAATGGCACCAGCGTGTCATCAGTCAGTGCGATGCTTTCAGCAACCTGTTGCTCTTCACGTACCAGTTCGATGAATGTGTCGTTGACCGGCAACCCATAGAACTGCGGGCCGTTTACAGAACAGAATGCTTCAAAGTGCTGCAAAGCATTCATTTCTTCAAAGACGGTAGCGTAACTGCCCAGCGCGGTTGGGGCGTTGAAGCAGCCCGCGCAGCCGCAACTGCTCTCTTTGCGATGACGTGCATGTGGCGCAGAATCCGTACCGAGGAATACTCGATTAAAACCGCTGGCGACCAGTTCACGCAATGCCTGTTGGTGAATATTACGTTTGAGGATGGGTAGACAATACAGGTGCGGACGCACGCCTCCAACCAGCATATGGTTGCGGTTAAACATCAGATGC  >  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: