Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,534 T→C S86S (AGT→AGC ispH → 4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091326,5340TC95.2% 103.8 / ‑3.7 42S86S (AGT→AGCispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (23/17);  total (25/17)
Fisher's exact test for biased strand distribution p-value = 5.06e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

GTTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGA  >  NC_000913/26288‑26769
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
gtttgCCAACCCGGGGGGTTTTTTTGCCGGGGTATCCCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACTAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCcgcg                                                                                                                                                                                                                                           <  2:314251/247‑1 (MQ=255)
            gcgTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGgtgt                                                                                                                                                                                                                               >  1:85968/1‑251 (MQ=255)
                                          tATCAGCATTATTGAAAACACACTGGCCATTTACGCCGCACCGATATATGTCCGTAACGAAGTGGTACATAACCGCTATGTGGTCGGTAGCATGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGACGTACCGGACGGCGCGATCCTGATTTTATCCGCACCCGGTGATTCTCAGGCGGTACGTAACGATGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCAGCTGTCCGCAGGTGACCaaa                                                                                                                                                                                                >  2:163117/1‑250 (MQ=255)
                                                          aaCGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGc                                                                                                                                                                               <  1:41269/251‑1 (MQ=255)
                                                              cgcTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTCTGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGAGTTTGATGCCACCTGTCCGCTGGTGACCAAGGTGCATCTGGAAGTCGCCcgg                                                                                                                                                                           >  2:339434/1‑250 (MQ=255)
                                                                     cATTTACGGCGCACCGATATATGCCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAgcc                                                                                                                                                                    >  1:238330/1‑249 (MQ=255)
                                                                            tgcgcATCGATATATGGCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGTGAAGTACCGGACGGCGCGATCCTTAGTTTCTCCGCACACGGTGTTTCTCAGGCGGTTCTTAACGAAGCAAAAAGCCGCGATTTGACGgt                                                                                                                                                                                                                                 <  2:214820/184‑1 (MQ=255)
                                                                            ggCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGgt                                                                                                                                                                                                                                 >  1:214820/1‑185 (MQ=255)
                                                                                 aCCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATATTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCACGATTTGACGGTGTTTGATGCCACCTGCCCGCTGGTGCCCAACGAGCATATGGAAGTCGCCCGCGCCATTCGCCGTGGCga                                                                                                                                                          >  2:103871/1‑249 (MQ=255)
                                                                                    gATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgcagcgt                                                                                                                                                     >  1:249332/1‑245 (MQ=255)
                                                                                           gTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATtct                                                                                                                                              >  1:143707/1‑251 (MQ=255)
                                                                                                          gTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCgg                                                                                                                               >  1:254319/1‑251 (MQ=255)
                                                                                                                aaCCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACTGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAATCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAGAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGAGTCTCATCGGCCACGCCGGGCatcc                                                                                                                         <  1:150993/251‑4 (MQ=255)
                                                                                                                     cTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTCTCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAACGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGGTTCTCATCGGCCACGCCGGGCATCCggc                                                                                                                      >  2:363155/1‑248 (MQ=255)
                                                                                                                        tgtgGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTgg                                                                                                                 >  1:303726/1‑251 (MQ=255)
                                                                                                                              cGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggta                                                                                                           <  1:295588/251‑3 (MQ=255)
                                                                                                                                               gTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagt                                                                                          >  1:225805/1‑251 (MQ=255)
                                                                                                                                                        ccTTTATTGAGCAGATTCTCGAAGTACCGGGCTGCGCGATCATGATTTTCTCCGCACACCGTGGTTCTCAGGCCGTACGTACAGAAGCAAAAAGCCGTGATTGGCCGGTGTTCGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAATGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                 <  1:255845/250‑1 (MQ=255)
                                                                                                                                                           ttATAGAGCAGATTAGCGAAGTACCGGACTGCGCGATCCGGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTGTACGGGTGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCAGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAgg                                                                               <  2:113837/250‑1 (MQ=255)
                                                                                                                                                               tGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGa                                                                           <  1:64902/250‑1 (MQ=255)
                                                                                                                                                                               gTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATc                                                          <  1:116129/251‑1 (MQ=255)
                                                                                                                                                                                   cGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCTTATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCg                                                      >  1:9277/1‑251 (MQ=255)
                                                                                                                                                                                            cgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCCGAACAGTACCCCGGAAAGGGGAATGTATCTGGTCGAATCGCCGGACGAtgtg                                             >  2:332989/1‑251 (MQ=255)
                                                                                                                                                                                             gATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGACCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCCGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGCATCGCCGGACGAtgtgg                                            >  2:198411/1‑250 (MQ=255)
                                                                                                                                                                                                ccTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGa                                         <  1:76728/251‑1 (MQ=255)
                                                                                                                                                                                                  tGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGaa                                        >  1:36155/1‑250 (MQ=255)
                                                                                                                                                                                                   gATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAACGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAAc                                      >  1:134307/1‑251 (MQ=255)
                                                                                                                                                                                                       ttCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGa                                   <  1:345866/250‑1 (MQ=255)
                                                                                                                                                                                                       ttCTCAGCACAAGGTGGTGCTCAGGCGGTACGATACGAGGCAAAAAGCTGCGATTTGACTGTGTTGGATGCCACCGATGCGGTGGTGACCAAAGTGCATATGGAAGTCTCCCGCGTCAGCCGTCGGGGCGAAGAGGCGATTCTCATCGGCCACTCCGGGCCTCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTGGAATCGCCGGCCGATGTGTGGAAACTGa                                   <  2:312589/250‑1 (MQ=255)
                                                                                                                                                                                                         ctcctcACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTTTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGCTGGGCAAGTACAGTAACCCGGAAGGGGGACTGGAGCTGGTCGAATCGCCggtcgtagagtggaacctgaatg                                >  2:314928/1‑230 (MQ=255)
                                                                                                                                                                                                         ctcCGCACTCGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCAGCGATTTGACGGTGTTTGATGCCACCGTGCCGCTGGTGACTAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCAGCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTAGAATCg                                                         <  2:110077/226‑1 (MQ=255)
                                                                                                                                                                                                         ctcCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCg                                                         >  1:110077/1‑226 (MQ=255)
                                                                                                                                                                                                            cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGACGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGACGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGACTGTATCTGGGCGAATCGCCGGACGATGTGTGGAAACTGTCGGTc                              >  2:80072/1‑250 (MQ=255)
                                                                                                                                                                                                              cacaCGTTGTTTCTCAGGCGGTACTTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCaaa                           <  2:85968/251‑1 (MQ=255)
                                                                                                                                                                                                                     tgtTTCTCAGGCGGTACGTAACGAAGCAAAAAACCGCGATTTGACGGAGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAAAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGCTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGa                      >  1:221242/1‑249 (MQ=255)
                                                                                                                                                                                                                       tttCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAag                    >  1:237944/1‑249 (MQ=255)
                                                                                                                                                                                                                           tcAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGc              >  1:110547/1‑251 (MQ=255)
                                                                                                                                                                                                                                  gTACGTAACGACGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCACCCGTCGTGGCGAAGAGCCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAGGGTACGAGGGGCCAGTACAGTAACCCGGAAGGGGGAATGTAACTGGTCGAATAGCCGGACGATGTGTGGAATCTGACGGACAACAAAGAAGGGaaacccccctt       >  2:32282/1‑242 (MQ=255)
                                                                                                                                                                                                                                    aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCAGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt      <  1:117016/250‑1 (MQ=255)
                                                                                                                                                                                                                                     cgcaaCGAAGCAAAAAGCCGCTAGTTGACGGTTGTTGATGCCACCGGTCCTCTGGTGACCAAAGTGAATATGGAAGCCGCCCGCGCCAGCCGGCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGTACGATGTGTGGAAACTGACGGTAAAAAACGAAGAGAATCTCTCCTTTAt    <  2:249332/248‑1 (MQ=255)
                                                                                                                                                                                                                                     cGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTAt    >  1:260422/1‑251 (MQ=255)
                                                                                                                                                                                                                                        aaCGAATCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATTGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTTGAATCTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGa  <  1:101901/250‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                             
GTTGGCCAACCCGCGTGGTTTTTGTGCCGGGGTAGACCGCGCTATCAGCATTGTTGAAAACGCGCTGGCCATTTACGGCGCACCGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGA  >  NC_000913/26288‑26769

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