Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 26,609 C→T R111R (CGC→CGT 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,6090CT96.6% 170.2 / ‑3.4 60R111R (CGC→CGTispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base T (26/31);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 4.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.65e-01

CGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAA  >  NC_000913/26371‑26838
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgaagag                                                                                                                                                                                                                           <  2:305220/251‑2 (MQ=255)
cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgaagag                                                                                                                                                                                                                           <  1:31484/251‑2 (MQ=255)
      atGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCATGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCgaagagtcgatt                                                                                                                                                                                                                     >  1:56617/1‑244 (MQ=255)
               aCGAAGTGGTACATAACCGCTATGTGGGCGATAGGTTGCGTGAGCGAGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGACTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCCTATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCAtaggc                                                                                                                                                                                                            >  1:272532/1‑247 (MQ=255)
                               ccGCTATGAGGACGATAGCTCGCGTGAGCGTGGGGCTATCTTTCTTGAGCAGATTAGCGAAGTACAGGACGGCGCGATCCTGATTGTCTTCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGAGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGGTGCGAAAAGTCGATTCTCATCGGCCACGCCGGGCATCCgg                                                                                                                                                                                            >  2:137683/1‑251 (MQ=255)
                                    attgggTCGATAGCTTTCGTGAGCGTGGGGTTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTGCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCCCCCGGAAg                                                                                                                                                                                         <  2:479571/245‑1 (MQ=255)
                                      gtgGTCGAGATCTTGCTTTAGCGTGGGGCTATCTTTATTGAGCAGATTATCGAAGTTCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTTCGTAACGAAGCAAAAAGCCTCGATTTGACGGTGTGTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGa                                                                                                                                                                                     <  2:541023/251‑1 (MQ=255)
                                             aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGCAGTACCGGCCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGCGACCAAAGGGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCGGTGGAAggtacg                                                                                                                                                                              >  1:651350/1‑247 (MQ=255)
                                               aGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtacgat                                                                                                                                                                            <  1:504640/251‑7 (MQ=255)
                                                       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)
                                                                                                                                                        gAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                                                                                >  1:392346/1‑125 (MQ=255)
                                                                                                                                                        gAAGCAAAAAGCCGCGATTTGACGGGGTTTGATGCCACCTTTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGCCGTCGTGGCTAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                                                                                <  2:392346/126‑2 (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)
                                                                                                                                                                                                                                                |                                                                                                                                                                                                                                     
CGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAGTCGCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAA  >  NC_000913/26371‑26838

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