Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,624 T→G S116S (TCT→TCG 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,6240TG95.9% 142.8 / ‑3.6 49S116S (TCT→TCGispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base G (25/22);  total (26/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.55e-01

GTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGA  >  NC_000913/26379‑26860
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
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)
                                                                                                                                                         agccgcgATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCt                                                                                <  1:216448/248‑1 (MQ=255)
                                                                                                                                                          gccgcgATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCt                                                                                                   <  1:382968/229‑1 (MQ=255)
                                                                                                                                                                   tGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGCGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGCAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGGATCGCCGGCCGAAGCGTGGAAACTGACGGTCCAATACGGAGAGAAGCTATCCTTTATGCCCCAGCCCACGCTGTCGGTGGAt                                                                      >  2:377044/1‑251 (MQ=255)
                                                                                                                                                                        gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGacac                                                                 >  1:37504/1‑251 (MQ=255)
                                                                                                                                                                                aTGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGAt                                                          <  1:100240/251‑1 (MQ=255)
                                                                                                                                                                                      ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATc                                                    <  1:198411/251‑1 (MQ=255)
                                                                                                                                                                                      ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGAATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGagc                                                    >  1:243327/1‑249 (MQ=255)
                                                                                                                                                                                       cTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTAAGATGGGCCAGTATAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGAAACGTCTGATGTGATCg                                                   >  1:42445/1‑251 (MQ=255)
                                                                                                                                                                                         gTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGAc                                                 <  1:294981/251‑1 (MQ=255)
                                                                                                                                                                                          ggccgttGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCg                                                   <  1:346440/243‑1 (MQ=255)
                                                                                                                                                                                             cGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGAcgcg                                              >  1:381444/1‑251 (MQ=255)
                                                                                                                                                                                                 ggTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGGGGATGACACGTCTGATGTGATCGACGCGCTGCg                                         >  1:156064/1‑251 (MQ=255)
                                                                                                                                                                                                        aaaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACg                                   >  1:321567/1‑251 (MQ=255)
                                                                                                                                                                                                                   ggAAGTCGCCCGCGCCAGACGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTAGCTGGTCGAATCGCCGGACGATGTGAGGAAAATGACGGACAAAAACGAAGAGAAGCTCTCCGTTAGAAACCAGACCACGCGGTCGGGGGATGAAACGTCTGATGTGATAGACGCGCTGCGTAAACGATTGCCGaaa                         >  2:350203/1‑250 (MQ=255)
                                                                                                                                                                                                                         cGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTc                   <  1:8192/250‑1 (MQ=255)
                                                                                                                                                                                                                                         cGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGa  <  1:339434/251‑1 (MQ=255)
                                                                                                                                                                                                                                           tGGCGAAGAGTCGATTCGCATCGGCCCCGCCGGGCATCCGGAAGTGGAAGGGACGATGGGCCAGTTCAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTTACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGATGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATg   <  2:64606/248‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
GTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGA  >  NC_000913/26379‑26860

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