Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,663 G→T G129G (GGG→GGT 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,6630GT92.2% 160.7 / 1.3 66G129G (GGG→GGTispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (29/30);  total (32/34)
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.58e-01

ATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAG  >  NC_000913/26416‑26892
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCAGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCGAGTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACa                                                                                                                                                                                                                                    <  1:143586/251‑1 (MQ=255)
      ggCGTGAGCGTGGGGCTTTCTTTATTGAGCAGATTAGCGAAGTCCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGc                                                                                                                                                                                                                              <  2:186182/250‑1 (MQ=255)
               gTGGGGCTATATTTATTGAGCATATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCGGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagt                                                                                                                                                                                                                     <  1:179530/251‑1 (MQ=255)
                         cTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGa                                                                                                                                                                                                             >  1:429125/1‑249 (MQ=255)
                             aTTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAgggg                                                                                                                                                                                                        <  1:56036/250‑1 (MQ=255)
                               tGAGCAGATGAGGGAGGTACCGGTCGGCGGGAGCCTTAGTGTCTCCGCACACGGTGTTTCTCAGGCGGGACGTAACGAAGCAAAAGGCCGCGATTTGTCGGTGTTTGATGCCACCTGTCCGCTGGGGACCAAAGTGCATATGGAAGTGGCCCGCGCCAGCGGGCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCGTCCGGAAGTGGGAGGTGCGATGGGCCAGTACAGTAACCCGGAAggggg                                                                                                                                                                                                       <  2:423598/249‑1 (MQ=255)
                                         agcgaatgaCCGGACGGGTCGATCCGGATTTGCGCCGCACACGGGTTTTCTCTGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAATTGCTTATGGAAGTCGCCCGCTCCAGCCGTCGTGGCGAAGAGTCGATTCTCAGCGGCCTCGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGTGGAAGGTATCTGGt                                                                                                                                                                                           <  2:207404/243‑1 (MQ=255)
                                         aGCGAAGTACCGGACGGCCCGATCCTGATTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGt                                                                                                                                                                                           <  1:403460/250‑1 (MQ=255)
                                          gCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTc                                                                                                                                                                                          <  1:192114/251‑1 (MQ=255)
                                              attACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTGCTCAGGCGGTACGTAACAAAGCAAAAAGCCGCGATTTGACGTTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACATTAACCCGGAAGGGAGAATTTATCTGGTCGAAt                                                                                                                                                                                      <  2:211996/249‑1 (MQ=255)
                                              aGTACCGGACGGCGCGATCCTGATTTTCTCCACACACGGTGTTTCTCAGGCGGAACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCACGCGCCAGCAGTCGTGGCGACGAGTCGGTTCTAATCGGCCACGCCGGGCATACGGAAGAGGAAGGTACGATGGGCCAGGACAGTAACCCGGAAGGGGGCATGTATCTGGTCGAAt                                                                                                                                                                                      >  1:400064/1‑251 (MQ=255)
                                                      aCGGCGCGATCCTGATTTCCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                              <  2:290766/251‑1 (MQ=255)
                                                           gcgatctggATTTTCGACGCACACGGTGTTTCGAAGGCGGTCCGTAACTAAGCAAAAAGCCGCGAGTTTACGGTGTTTGATGCCACCTGTCCGCGGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGCACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                                                                                                                                         <  2:238443/243‑1 (MQ=255)
                                                                          cccGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTTCGCTGGTTACCAAAGGGCATATGGAAGTCGCCCGCGCCAGCCGGCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGt                                                                                                                                                          <  2:247639/250‑1 (MQ=255)
                                                                           ccGCACACGGTGTTTCTCAGGCGGTACTTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGTCGGACGATGTGTGGAAACTGACGGTc                                                                                                                                                         <  2:334994/251‑1 (MQ=255)
                                                                               acacGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCCTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGGGGAAACTGACGGTCaaa                                                                                                                                                      >  2:452763/1‑250 (MQ=255)
                                                                                  cGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACg                                                                                                                                                  <  1:346886/251‑1 (MQ=255)
                                                                                    ggttttCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGaa                                                                                                                                                <  2:349340/248‑1 (MQ=255)
                                                                                          ctcAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGGCGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTAGGTTGGGCCAGTACAGTTACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGGAACTGACGGTCCAAAACGAAgaga                                                                                                                                            >  2:279977/1‑249 (MQ=255)
                                                                                              ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCAGAGCATCCGGAAGTGGAAGGTACGATTGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctcc                                                                                                                                      >  1:156564/1‑250 (MQ=255)
                                                                                                 ggTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGGGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCt                                                                                                                                   >  2:233290/1‑251 (MQ=255)
                                                                                                   ggcGTAACGAAGCAAAAAGCCGGGATTGGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGCAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                 <  2:449416/249‑1 (MQ=255)
                                                                                                    aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                 <  1:418132/250‑1 (MQ=255)
                                                                                                     cgaaaCGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGTATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATGATCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTTGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTAt                                                                                                                               <  2:156564/248‑1 (MQ=255)
                                                                                                     cGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTAt                                                                                                                               <  2:152340/251‑1 (MQ=255)
                                                                                                        aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACt                                                                                                                                                                <  2:372476/215‑1 (MQ=255)
                                                                                                        aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACt                                                                                                                                                                >  1:372476/1‑215 (MQ=255)
                                                                                                              gCAAAAATCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAACGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGGACGATGGGACAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCACTCCTTTATGACCCAGa                                                                                                                       >  2:422743/1‑250 (MQ=255)
                                                                                                                  aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGCAGCTCTCCTTTATGACCCAGAccccg                                                                                                                  >  2:138097/1‑248 (MQ=255)
                                                                                                                     gccgcgATTTGACGGGGTTTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTg                                                                                                               <  2:340042/249‑1 (MQ=255)
                                                                                                                         cgATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCgg                                                                                                           >  1:331912/1‑251 (MQ=255)
                                                                                                                           aTTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTg                                                                                                         >  1:294900/1‑251 (MQ=255)
                                                                                                                                   gtgtTTGATGCAACCTGTCCGCTGGTGACCCAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGGTGTGTGGAAACTTACTGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGagtc                                                                                                 >  2:2683/1‑248 (MQ=255)
                                                                                                                                           tGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGAtg                                                                                         >  2:236040/1‑251 (MQ=255)
                                                                                                                                               cccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGAt                                                                                     <  2:221344/250‑1 (MQ=255)
                                                                                                                                                ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGACGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGAAACTGACGGGCAAAAACGAAGCGAAGCTCTCCTTTATGACCCAGACCACGCAGTCGGCGGATGACACGACTGATGTGATc                                                                                    >  2:222308/1‑251 (MQ=255)
                                                                                                                                                ccTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATc                                                                                    >  1:395063/1‑251 (MQ=255)
                                                                                                                                                   gTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAgaga                                                                                                                                            >  1:457606/1‑192 (MQ=255)
                                                                                                                                                          ggTGACCAAATTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCAAAGAGTCGATTCGCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCt                                                                            <  1:232211/249‑1 (MQ=255)
                                                                                                                                                            tGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGGATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAACAACGAAGAGAAGCTCTCCTTTATGGCCCAGACCAAGCTGTCGGTGGATGCCAAGTCTGATGTGATCGACTCGCTGCGt                                                                        >  2:119095/1‑251 (MQ=255)
                                                                                                                                                                  aaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACCGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGCCGGTGAGGGGAAACTGCCGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGCCCACGCTGCCGGGGGATGACCAGTCTTATGTGATCgcgtcgacgcgtaaaccc                                                                  >  2:275136/1‑234 (MQ=255)
                                                                                                                                                                  aaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGGACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAACACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCGGCTGTGTTCGACGCGCTGCGTAAACGc                                                                  >  2:73447/1‑251 (MQ=255)
                                                                                                                                                                      gCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTcc                                                              <  1:236040/251‑1 (MQ=255)
                                                                                                                                                                       cATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTAGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGAGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTccc                                                             >  1:184557/1‑251 (MQ=255)
                                                                                                                                                                         tatGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCg                                                            >  1:261290/1‑250 (MQ=255)
                                                                                                                                                                         tatGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCAGGGCATCCGGGAGTGGCAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACCCGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCg                                                            >  1:198622/1‑250 (MQ=255)
                                                                                                                                                                              aaGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCCCGCTGTCGGTGGATGACACGTCTGATGTGATCGACTCGCTGCGTAAACGCTTCCCGAAAAtt                                                      <  1:138097/251‑1 (MQ=255)
                                                                                                                                                                                   ccccGCGCCAGTCGTCGTGGGTAAGAGTAGATTCTCATCGGCCACGCCGGGCATCCGGAAGATGAAGGTACGATGGGCCTGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTTCTTTATGACCCCGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCgt                                                 <  2:33816/250‑2 (MQ=255)
                                                                                                                                                                                    cccGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACCCGTCTGATGTGGTCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCggg                                                >  2:483825/1‑250 (MQ=255)
                                                                                                                                                                                          ccagccgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGCATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCgcgc                                          >  2:170254/9‑251 (MQ=255)
                                                                                                                                                                                             gccgtcGTGGCGAAGAGGCGATTCGCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCaaa                                       <  2:11148/246‑1 (MQ=255)
                                                                                                                                                                                               cgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGACGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGGGGATGACACGTCTGAAGTGATCGACGCACTGCGTAAACGATTCCCGAAAATTGTCGGTCCGCGCacagg                                     >  1:325957/4‑247 (MQ=255)
                                                                                                                                                                                                 gcGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGGCGAAGCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTTTCTGTGGATGACTCGTCTGATGTGATCGACGCGCTGCGGAAACGCTTCCCGTAAATTGTCGGTCCGCGCAAAGATg                                   <  2:235848/250‑1 (MQ=255)
                                                                                                                                                                                                   gTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCCAAAATGAc                                 >  2:408050/1‑251 (MQ=255)
                                                                                                                                                                                                   gTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGAc                                 >  1:16220/1‑251 (MQ=255)
                                                                                                                                                                                                      gCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCTCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCGTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATc                              <  2:313940/251‑1 (MQ=255)
                                                                                                                                                                                                       cGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGACAATTGTCGGTCCGCGCAAAGATGACATCt                             >  1:36467/1‑251 (MQ=255)
                                                                                                                                                                                                        gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTg                            >  1:41687/1‑251 (MQ=255)
                                                                                                                                                                                                        gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTg                            <  1:452763/251‑1 (MQ=255)
                                                                                                                                                                                                        gaagagTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGACCGATCGGCCAGTACAGAAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGCGGAAACTGACGGTCAAAAACGAAGAGACGCTCTCCTATATGACCCAGACCACGCTGTCGGTGGATTACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCAGCGCAAAGATGACAtgtg                            >  1:198095/1‑248 (MQ=255)
                                                                                                                                                                                                                 aTTCTCATCGGCCACGCCGTGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCAc                   <  1:272539/251‑1 (MQ=255)
                                                                                                                                                                                                                      cATCGGCCACGCCGGGCCTCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTa              <  2:261290/251‑1 (MQ=255)
                                                                                                                                                                                                                       aTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGCCCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACAGCCGCCACGCCACGACTa              >  2:10230/1‑250 (MQ=255)
                                                                                                                                                                                                                         cgggcACGCCGGGCTGCCGGAAGTGGAAGGTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGGCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGTGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAAcc           <  2:138882/247‑1 (MQ=255)
                                                                                                                                                                                                                               cGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAgg     <  2:452882/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGGCATCCGGAAGTGGAAGGTACGGTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGCTGTGTGGAAACTGACGGGCAAAAACGAAGAGAAGCTCTCCTTTATGACCCCGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTAAGCCACGACTAACCGTccgggtg  >  2:247744/1‑245 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                     
ATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAG  >  NC_000913/26416‑26892

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