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,6630GT95.7% 208.6 / NA 71G129G (GGG→GGTispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  major base T (33/34);  minor base A (3/0);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 2.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGTACGCGCCCTG  >  NC_000913/26426‑26906
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
ggAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagt                                                                                                                                                                                                                                        <  1:290244/250‑1 (MQ=255)
  ggCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtac                                                                                                                                                                                                                                      <  1:1879/250‑1 (MQ=255)
       ggggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagta                                                                                                                                                                                                                                  >  1:96062/1‑250 (MQ=255)
        gggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAc                                                                                                                                                                                                                                >  1:58327/1‑251 (MQ=255)
         ggCTATCTTTATTGGGCAGATTAGCGAAGTACCGGCCGGCGCGATCCTGATTTTCCCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGGCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAcc                                                                                                                                                                                                                               <  2:43184/251‑1 (MQ=255)
                        gCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACTGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGt                                                                                                                                                                                                                <  1:51101/251‑1 (MQ=255)
                           gATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATc                                                                                                                                                                                                             <  1:233382/251‑1 (MQ=255)
                            tttAGAGAAGAACCGGACGGCGCGATCCTGCTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAAAGAAGCAAAAAGCCGCGATTTGACGGTGTTTGCTGCCACCTGTCCGCTGGTGACCAACGTGAATCTGGAAGTAGCACGCGCCAGCCGTCGTGGCGGAGAGTCGCTTCTCATAGGCCACGCCGGGCATCCGGAAGTGGAAGGTACAATGGGCCAGCACAGTAACCCGGCAggggggctgattct                                                                                                                                                                                                            >  2:239633/2‑242 (MQ=255)
                                 cGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCg                                                                                                                                                                                                       <  1:204726/251‑1 (MQ=255)
                                     gTACCGGACGGCGCGATCCTTATTTTCTCCGCACACGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATc                                                                                                                                                                                                   <  2:194581/251‑1 (MQ=255)
                                      gaCCGGACGGCGAGATTCTGTTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTTGCGTTGTGGCGGTGTTTGATGCCACCGGTCCGCTGGGGTCCAAAGGGCATATGGAAGCCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTCCAGTAACCCGGAAGGGGGAATGTATCTGGTAGaagcg                                                                                                                                                                                                  <  2:376735/250‑4 (MQ=255)
                                       aCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCg                                                                                                                                                                                                  <  1:17025/250‑1 (MQ=255)
                                         cGGACAGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCg                                                                                                                                                                                               >  1:489675/1‑251 (MQ=255)
                                           gACGGCGCGATCCGTATTTCTTCCGCACCCGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGAGTCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCTCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCAGGCCGGGCATCCGGAAGTGGAAGGTACGATGTGCCAGTACAGTAGCCCGGAAGGGGGAATTTATCTGGTCGAATTGCCGGa                                                                                                                                                                                             <  2:640303/251‑1 (MQ=255)
                                            aCGGCGCGATCCTGATTTTCTCCGCTCACGTTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                                            <  1:524642/251‑1 (MQ=255)
                                            aCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                                                                                            >  1:511134/1‑251 (MQ=255)
                                                cgcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtg                                                                                                                                                                                        <  1:618366/251‑1 (MQ=255)
                                                 gcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                                                                                                                                                       >  1:194208/1‑251 (MQ=255)
                                                   gATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGGATGTATCTGGTCGAATCGCCGGACGAtgtgt                                                                                                                                                                                     >  1:567684/1‑251 (MQ=255)
                                                          atttctgccGCACACGGTGTTTCTCAGGCGGTCCGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAGCCGTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGGATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAAca                                                                                                                                                                              <  2:33752/244‑2 (MQ=255)
                                                           atttatgccgactaggttgttCTTAGGCGGGACGTAAAGAAGCAAAAAGCAGCGATTTGACGGGGTTTTATGCAACCGGCCTGCTGGTGCCCAAAGTGCATATGGAAGGCGCCCGCGCCAGCTTGCAGGGCGAAGAGTCGATTCCCAACGGACACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCTGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTGGAAACTg                                                                                                                                                                             <  2:126036/232‑1 (MQ=255)
                                                            tttCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTg                                                                                                                                                                             <  1:16671/250‑1 (MQ=255)
                                                                  cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTc                                                                                                                                                                       >  1:589189/1‑250 (MQ=255)
                                                                  cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCa                                                                                                                                                                      <  1:324194/251‑1 (MQ=255)
                                                                    cacaCGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGTAGTGGTAGGTACGATGGGCCAGTACAGTAACCCGGGAGGGGGAATGTATCTGGTCGAATCGCCGGACGAGTTGGGGGAACTAACGGTCaa                                                                                                                                                                     >  2:253061/1‑250 (MQ=255)
                                                                       acGGTGGTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTGTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAc                                                                                                                                                                 <  2:58327/251‑1 (MQ=255)
                                                                           tgtTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGa                                                                                                                                                               >  1:78523/1‑249 (MQ=255)
                                                                                   aGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGAAGctc                                                                                                                                                     >  1:483803/1‑251 (MQ=255)
                                                                                    ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctct                                                                                                                                                    <  1:382809/251‑1 (MQ=255)
                                                                                      gtgtatgTAACCAAGCAAAAAGCCGTGAGTTGAAGTTGTTTGAGGCACCCTGTCCGCTGGTGCCCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGTAGTTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGGGTGGAAACTGACGGTCAAAATCGAAGAGAAGctctac                                                                                                                                                  <  2:479171/245‑3 (MQ=255)
                                                                                          accgaaCGAAGCAAAAAGCCGTGATTTGACGGTTTTTGATGCCGCCTTGCCGCTGGTGTCCAAATTGCATTTGGAAGTCGCCCGCGCCAGCCGTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                                                                                               <  2:430311/246‑1 (MQ=255)
                                                                                          aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTa                                                                                                                                              <  1:146492/251‑1 (MQ=255)
                                                                                              aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGAc                                                                                                                                          >  1:397084/1‑251 (MQ=255)
                                                                                                       aaaaGCCGCGATTGGACTGTGTTTGTTGCTACCTGGCCGCTGGTGACCAAAGTGCTTATGGAAGTCGCCCGCGCCAGCCGGCGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTAGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAc                                                                                                                                 <  2:356409/251‑1 (MQ=255)
                                                                                                        aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCAACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGAcaa                                                                                                                                  >  1:389067/1‑247 (MQ=255)
                                                                                                                gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTc                                                                                                                           <  1:103076/248‑1 (MQ=255)
                                                                                                                gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCCCTCCTTTCTGACCCAGACCAcactttc                                                                                                                           >  2:103076/1‑242 (MQ=255)
                                                                                                                     gACGGTGCTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGACGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGCAGGTACGATGGGCCAGTACCGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCAGGACGATGCGGGGAAACTGACGGTCAAAAACGAAGAAAAGCTCTCCTTTATGACCCAGACCACGCTTTCGGTGGATg                                                                                                                   >  2:252639/1‑251 (MQ=255)
                                                                                                                       cGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGac                                                                                                                 >  1:645707/1‑251 (MQ=255)
                                                                                                                         gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                                                                                          <  2:158486/144‑1 (MQ=255)
                                                                                                                         gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                                                                                          >  1:158486/1‑144 (MQ=255)
                                                                                                                              ttaTGCCCCCTGTCCGGTGGTGACCAAATTGATTATGTAAGCCGCCCCCGCCGTCCGTCGTGGCGAAGAGCCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAATGTGCGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGTCGATGTGGGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACTCGTCt                                                                                                           <  2:122761/248‑1 (MQ=255)
                                                                                                                               gATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGa                                                                                                         <  1:480145/251‑1 (MQ=255)
                                                                                                                                             gCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAag                                                                                                                                                             >  1:126007/1‑185 (MQ=255)
                                                                                                                                             gCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGACACGCCGGGCATCCGGACGTGGAAGGAACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAAGAAGAGAAGCTCTCCTTTATGAACCAGACCACGCTGTCGGTGGAAGACACGtaggttgtgaacgccgggc                                                                                           >  1:669822/1‑233 (MQ=255)
                                                                                                                                                        aaGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCAGTGGAAGGTACGCTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGGGTGGAAACTGACGGTCAAAAACGAAGAGACGCTCTCCTTTTTGACCCAGAACACGCTGTCGGTGGATGACACGTCTGATGTGATAGACGCGCTGCGGAAACGc                                                                                >  2:257378/1‑251 (MQ=255)
                                                                                                                                                         aGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCt                                                                               <  1:38673/251‑1 (MQ=255)
                                                                                                                                                          gTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGGTGGGCCAGTACAGGAACCCGGAAGGGGGAATGGATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCATTTATGAACCAGACCACGCGATCGGAGGGTGACCCTTCTGGTGTGATCGAAGGGCGGCGTAGACGCtt                                                                              >  2:465199/1‑251 (MQ=255)
                                                                                                                                                                  ggAAGTCGGCCCCGCCAGCCGTCGGTGCGGAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGTTGGGGCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTTGAAACTGACGGTCAAAAACGAAGAGAATCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGaaaa                                                                      <  2:558736/251‑1 (MQ=255)
                                                                                                                                                                                   cagccgtcGTGGCGAAGCGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGAACGATGGGCCAGTACAGTCAACCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAAGCAGTCGGTGGATGACACGTCTGATGTGATCGCAGCACTGCGTTAACGCTTCCCGAAACTTGTCGGACCgcgc                                                        >  2:576564/8‑250 (MQ=255)
                                                                                                                                                                                   cagccgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCgcgc                                                        >  1:546223/8‑250 (MQ=255)
                                                                                                                                                                                    ggccgtcggggCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGGGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCaa                                                      <  2:546223/242‑1 (MQ=255)
                                                                                                                                                                                       cgtcGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGCGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGAAACGTCTGATGGGATCGACGCGCTGCGTAAACGCTTCCCGCAAATTGTAGGTCCGCGCaa                                                      >  2:337601/4‑248 (MQ=255)
                                                                                                                                                                                          cGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGGCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCGTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATg                                                 <  2:397084/250‑1 (MQ=255)
                                                                                                                                                                                            gggCGAAGAGTCGATTCTCATCGGCCACGCCGGGCAGCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGTAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACa                                              <  2:496297/250‑1 (MQ=255)
                                                                                                                                                                                                gaagagTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCGGATGTGATCGACGCGCTGCGTAAACGCTTGCCGAAAATTGTCGGTCCTCGCAATGATGGCATCTg                                          <  2:287409/251‑1 (MQ=255)
                                                                                                                                                                                                  agagtcgaTTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGCTCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCt                                        >  1:100858/8‑251 (MQ=255)
                                                                                                                                                                                                        gaTTCTCATCGGCCACGCCGGGAATCCGGAAGAGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTACGCGCAAAGATGACATCTGCTACGccc                                  >  2:455550/2‑250 (MQ=255)
                                                                                                                                                                                                              cATCGGCCACGCCGGGCATCCGGAAGTGGCAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACt                             >  1:529457/1‑250 (MQ=255)
                                                                                                                                                                                                              cATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACt                             <  1:337601/250‑1 (MQ=255)
                                                                                                                                                                                                               aTCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTa                            >  1:398895/1‑250 (MQ=255)
                                                                                                                                                                                                                   gTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACTTCTTCTCCGCCACTACTAAccat                       >  2:215164/1‑249 (MQ=255)
                                                                                                                                                                                                                        gCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAgg                   >  1:571500/1‑250 (MQ=255)
                                                                                                                                                                                                                            ggCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCg              >  1:89839/1‑251 (MQ=255)
                                                                                                                                                                                                                             gcatccGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGAGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCgg             >  1:22774/5‑251 (MQ=255)
                                                                                                                                                                                                                               atccGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGt            >  1:36942/3‑250 (MQ=255)
                                                                                                                                                                                                                               atccGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGt            >  2:226569/3‑250 (MQ=255)
                                                                                                                                                                                                                                 ccGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGTAcg         <  1:4780/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGCCGCCCTGCGTACACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTACCCGTCAGGAAGCGGTagcc        >  2:204957/1‑248 (MQ=255)
                                                                                                                                                                                                                                  cGGAAGTGGAAGGAACGATAGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGAACCAGACCACGCTGTCGGTGAATGACACGTCTGATGTGATCGACTCTCTGCGTAAACGCTTCCCGAAACTTGGCGGTCCGCGCAAAGAAGACATCTGCATCGCCACGACTACCCGTCAGGAAGCGGTAcgc        >  2:574968/1‑251 (MQ=255)
                                                                                                                                                                                                                                        tGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGTACGCGCCCTg  >  1:658566/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                   
TGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAATTGTCGGTCCGCGCAAAGATGACATCTGCTACGCCACGACTAACCGTCAGGAAGCGGTACGCGCCCTG  >  NC_000913/26426‑26906

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: