Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,148,900 G→A G335D (GGC→GAC)  gpr → L‑glyceraldehyde 3‑phosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,148,9000GA83.6% 152.3 / 25.2 67G335D (GGC→GAC) gprL‑glyceraldehyde 3‑phosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (4/7);  new base A (23/33);  total (27/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.41e-01

ACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCG  >  NC_000913/3148659‑3149137
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             
acacCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGa                                                                                                                                                                                                                                       >  2:218316/1‑250 (MQ=255)
  acCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATc                                                                                                                                                                                                                                    >  2:57416/1‑251 (MQ=255)
              tACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGTGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGTAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGc                                                                                                                                                                                                                          <  1:137011/249‑1 (MQ=255)
                 cGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCtt                                                                                                                                                                                                                     <  2:147853/251‑1 (MQ=255)
                      ccAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGAt                                                                                                                                                                                                                >  2:78179/1‑251 (MQ=255)
                      ccAACCTCAACAGCCTGAGGTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCGGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATAGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGAt                                                                                                                                                                                                                <  1:334777/251‑1 (MQ=255)
                       cAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGAt                                                                                                                                                                                                                >  2:116450/1‑250 (MQ=255)
                         aCCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGTACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATaa                                                                                                                                                                                                              >  1:141238/1‑250 (MQ=255)
                           cTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATg                                                                                                                                                                                                           >  1:139798/1‑251 (MQ=255)
                             cAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGAc                                                                                                                                                                                                         <  1:306244/251‑1 (MQ=255)
                                      gcgcTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTtaata                                                                                                                                                                                                <  1:218316/251‑1 (MQ=255)
                                      gcgcTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTtaat                                                                                                                                                                                                 <  1:486621/250‑1 (MQ=255)
                                       cgcTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTtaataa                                                                                                                                                                                               <  1:151362/251‑1 (MQ=255)
                                            aTTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGgc                                                                                                                                                                                          <  1:122668/251‑1 (MQ=255)
                                              tgaatgaaATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGGCCTGTTAATAACGGgcg                                                                                                                                                                                         <  1:388256/250‑1 (MQ=255)
                                               gaatgaaATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGgcgcg                                                                                                                                                                                       <  2:178030/251‑1 (MQ=255)
                                                aatgaaATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGgcgcgc                                                                                                                                                                                      >  1:279079/1‑251 (MQ=255)
                                                        ggCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACcgcg                                                                                                                                                                              <  2:279079/251‑1 (MQ=255)
                                                           acaGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGccc                                                                                                                                                                           <  2:91628/251‑1 (MQ=255)
                                                           acaGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGGTAATAACGGGCGAGCGAACCGCGccc                                                                                                                                                                           >  1:275499/1‑251 (MQ=255)
                                                                agcGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatc                                                                                                                                                                      >  2:390212/1‑251 (MQ=255)
                                                                  cGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatc                                                                                                                                                                      >  1:83999/1‑249 (MQ=255)
                                                                   gTGGACAATCAATGGCGCAAATGGCGTTACGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAatc                                                                                                                                                                      >  1:30371/1‑248 (MQ=255)
                                                                       aCAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGTTTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGCCATTTATCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATTACCTGTTAATACCGGGGGCGCGAACCGCGCCCACatcagatcac                                                                                                                                                               >  2:187830/1‑250 (MQ=255)
                                                                                   gcAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGaa                                                                                                                                                                                   >  2:364096/1‑219 (MQ=255)
                                                                                   gcAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGaa                                                                                                                                                                                   <  1:364096/219‑1 (MQ=255)
                                                                                    cAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGag                                                                                                                                                  <  2:319974/251‑1 (MQ=255)
                                                                                         ggCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTCTGACGagtcagt                                                                                                                                             >  2:71269/1‑251 (MQ=255)
                                                                                           ccttAAGCTGGTTGCTGATAGATGAGGGGTTGACGTCGGTGTTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGGACGTGCAGGCGCTGATTAATCTGTCATTTAGCACCAAGGAGCTGGCGTAGATTGATCAACATATCGCCGTTGACGAGCTGAATCTGTGGCAGGCTGCTTCCGATAAATCACCTGTTAGTATCGTGCGAGTGAACCGGGCCCGAATCAGAGCAGTAATTATGACAAGTCATTCg                                                                                                                                           <  2:278102/249‑1 (MQ=255)
                                                                                            gTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATTACGAGCTGAATCTGTGGc                                                                                                                                                                                                                              >  2:165916/1‑167 (MQ=255)
                                                                                            gTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGc                                                                                                                                                                                                                              <  1:165916/167‑1 (MQ=255)
                                                                                                    ggTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGAt                                                                                                                                   <  1:207122/250‑1 (MQ=255)
                                                                                                         cTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGAGCAGATTAATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGAGCGCGAATCGCGCCCGAATCAGGTCAGTGATTATGACAAGTCAGTCGATCAAGAAAGccc                                                                                                                              >  1:70750/1‑250 (MQ=255)
                                                                                                          tGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGGATCAGATCAGTGATTATGACGAGTCAGTCGGACAAGATAGCCCAt                                                                                                                            >  2:71646/1‑251 (MQ=255)
                                                                                                           gAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATc                                                                                                                           <  1:390212/251‑1 (MQ=255)
                                                                                                              aGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCTATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAAGGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATcaca                                                                                                                        <  1:50645/251‑1 (MQ=255)
                                                                                                                 tgatCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAAAGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACaaa                                                                                                                      >  2:229588/1‑250 (MQ=255)
                                                                                                                   atCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGc                                                                                                                   <  2:145021/251‑1 (MQ=255)
                                                                                                                   atCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGc                                                                                                                   <  1:136454/251‑1 (MQ=255)
                                                                                                                       cgTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCTCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCTGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGaa                                                                                                               <  2:226893/251‑1 (MQ=255)
                                                                                                                                   attgagtgGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCATGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACTAACGTc                                                                                                   <  2:70960/251‑1 (MQ=255)
                                                                                                                                      gattgGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAg                                                                                                 <  1:15318/250‑1 (MQ=255)
                                                                                                                                          ggTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATgg                                                                                             <  1:16325/250‑1 (MQ=255)
                                                                                                                                          ggTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATgg                                                                                             <  1:3550/250‑1 (MQ=255)
                                                                                                                                                 gCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTtaataa                                                                                                                                                                                               >  1:194526/1‑145 (MQ=255)
                                                                                                                                                 gCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTtaataa                                                                                                                                                                                               <  2:194526/145‑1 (MQ=255)
                                                                                                                                                       cgGAGCAACTAGAGGAGCACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGc                                                                                                                                                                                                                          >  2:441238/1‑112 (MQ=255)
                                                                                                                                                       cgGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGc                                                                                                                                                                                                                          <  1:441238/112‑1 (MQ=255)
                                                                                                                                                         gAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACcac                                                                              >  2:142624/1‑250 (MQ=255)
                                                                                                                                                              aCTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGt                                                                        >  1:423127/1‑251 (MQ=255)
                                                                                                                                                                 agagGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGtt                                                                       <  2:139798/249‑1 (MQ=255)
                                                                                                                                                                        aaCGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCGGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCg                                                              >  1:148890/1‑251 (MQ=255)
                                                                                                                                                                           gTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCggg                                                            >  1:307759/1‑250 (MQ=255)
                                                                                                                                                                                   gcTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTAACCGGGACATTTTTc                                                   <  1:71269/251‑1 (MQ=255)
                                                                                                                                                                                        ataatCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATc                                              >  2:416699/1‑251 (MQ=255)
                                                                                                                                                                                          aatCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTTAATCTGTGGCAGGCTTCTTCCGATAAATGACCTGTTAATAACTGGCGCGCGAACCGCGCCCGAATCAGAGCAGTGATTATGACGAGTCAGTCGTTTAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATcc                                             <  1:205184/250‑1 (MQ=255)
                                                                                                                                                                                                 cATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAc                                     <  1:416699/251‑1 (MQ=255)
                                                                                                                                                                                                   tttAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAcgc                                   >  1:424521/1‑251 (MQ=255)
                                                                                                                                                                                                     tAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGAcgcg                                  <  1:71844/250‑1 (MQ=255)
                                                                                                                                                                                                       gCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGGATCTTACGATAAATAACCTGTTAAGAGCGGGCGCGCGAGCCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTAAAGATAGCCCATCTCAAAGGCAGAAAGCACAATCGTCAGATGGATAATGACGTACCGCACCAGTTTGTTATCCGGGGCATTTTTCGCATCCAGAAAGACGCGCAg                               <  2:197879/251‑1 (MQ=255)
                                                                                                                                                                                                            aaGGAGCTGGCGCAGAGTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGt                          <  1:460444/251‑1 (MQ=255)
                                                                                                                                                                                                                    ggCGCAGATTGATCAGCATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGgaa                  <  1:57416/251‑1 (MQ=255)
                                                                                                                                                                                                                         aGATTGATCAGCATATCGGCGATGACGAGCTGAAGCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACAGCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAAt             <  1:481753/251‑1 (MQ=255)
                                                                                                                                                                                                                                  aGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAAt    <  2:83999/251‑1 (MQ=255)
                                                                                                                                                                                                                                    cATATCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCg  <  1:78179/251‑1 (MQ=255)
                                                                                                                                                                                                                                      tatCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCAt                                                                                                                            >  2:258211/1‑127 (MQ=255)
                                                                                                                                                                                                                                      tatCGCCGATGACGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCAt                                                                                                                            <  1:258211/127‑1 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                             
ACACCGAAAATGCTTACCGAAGCCAACCTCAACAGCCTGCGCTTATTGAATGAAATGGCACAGCAGCGTGGACAATCAATGGCGCAAATGGCGTTAAGCTGGTTGCTGAAAGATGATCGCGTGACGTCGGTATTGATTGGTGCCAGCCGCGCGGAGCAACTAGAGGAGAACGTGCAGGCGCTGAATAATCTGACATTTAGCACCAAGGAGCTGGCGCAGATTGATCAGCATATCGCCGATGGCGAGCTGAATCTGTGGCAGGCGTCTTCCGATAAATGACCTGTTAATAACGGGCGCGCGAACCGCGCCCGAATCAGATCAGTGATTATGACGAGTCAGTCGGTCAAGATAGCCCATCACAAATGCAGAAAGCACAAACGTCAGATGGATAATGACGTACCACATCAGTTTGTTATCCGGGACATTTTTCGCATCCATAAAGACGCGCAGTAAGTGAATGGAAGAAATTGCCACAATCG  >  NC_000913/3148659‑3149137

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