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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091326,6060TC95.0% 160.2 / ‑2.7 60S110S (AGT→AGC) ispH4‑hydroxy‑3‑methylbut‑2‑enyl diphosphate reductase
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (26/31);  total (26/34)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.09e-01

CGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAA  >  NC_000913/26371‑26838
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         
cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCgaagag                                                                                                                                                                                                                           <  2:305220/251‑2 (MQ=255)
cGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCgaagag                                                                                                                                                                                                                           <  1:31484/251‑2 (MQ=255)
      atGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCATGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCgaagagtcgatt                                                                                                                                                                                                                     >  1:56617/1‑244 (MQ=255)
               aCGAAGTGGTACATAACCGCTATGTGGGCGATAGGTTGCGTGAGCGAGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGACTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCCTATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCAtaggc                                                                                                                                                                                                            >  1:272532/1‑247 (MQ=255)
                               ccGCTATGAGGACGATAGCTCGCGTGAGCGTGGGGCTATCTTTCTTGAGCAGATTAGCGAAGTACAGGACGGCGCGATCCTGATTGTCTTCGCACACGGGGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGAGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGGTGCGAAAAGTCGATTCTCATCGGCCACGCCGGGCATCCgg                                                                                                                                                                                            >  2:137683/1‑251 (MQ=255)
                                    attgggTCGATAGCTTTCGTGAGCGTGGGGTTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTGCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCCCCCGGAAg                                                                                                                                                                                         <  2:479571/245‑1 (MQ=255)
                                      gtgGTCGAGATCTTGCTTTAGCGTGGGGCTATCTTTATTGAGCAGATTATCGAAGTTCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTTCGTAACGAAGCAAAAAGCCTCGATTTGACGGTGTGTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGa                                                                                                                                                                                     <  2:541023/251‑1 (MQ=255)
                                             aTAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGCAGTACCGGCCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGCGACCAAAGGGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCGGTGGAAggtacg                                                                                                                                                                              >  1:651350/1‑247 (MQ=255)
                                               aGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAggtacgat                                                                                                                                                                            <  1:504640/251‑7 (MQ=255)
                                                       ggAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagt                                                                                                                                                                    <  1:290244/250‑1 (MQ=255)
                                                         ggCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtac                                                                                                                                                                  <  1:1879/250‑1 (MQ=255)
                                                              ggggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCcagtacagta                                                                                                                                                              >  1:96062/1‑250 (MQ=255)
                                                               gggCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAc                                                                                                                                                            >  1:58327/1‑251 (MQ=255)
                                                                ggCTATCTTTATTGGGCAGATTAGCGAAGTACCGGCCGGCGCGATCCTGATTTTCCCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGGCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAAcc                                                                                                                                                           <  2:43184/251‑1 (MQ=255)
                                                                               gCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACTGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGt                                                                                                                                            <  1:51101/251‑1 (MQ=255)
                                                                                  gATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATc                                                                                                                                         <  1:233382/251‑1 (MQ=255)
                                                                                   tttAGAGAAGAACCGGACGGCGCGATCCTGCTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAAAGAAGCAAAAAGCCGCGATTTGACGGTGTTTGCTGCCACCTGTCCGCTGGTGACCAACGTGAATCTGGAAGTAGCACGCGCCAG‑CC‑GTCGTGGCGGAGAGTCGCTTCTCATAGGCCACGCCGGGCATCCGGAAGTGGAAGGTACAATGGGCCAGCACAGTAACCCGGCAggggggctgattct                                                                                                                                        >  2:239633/2‑242 (MQ=255)
                                                                                        cGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCg                                                                                                                                   <  1:204726/251‑1 (MQ=255)
                                                                                            gTACCGGACGGCGCGATCCTTATTTTCTCCGCACACGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATGTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATc                                                                                                                               <  2:194581/251‑1 (MQ=255)
                                                                                             gaCCGGACGGCGAGATTCTGTTTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTTGCGTTGTGGCGGTGTTTGATGCCACCGGTCCGCTGGGGTCCAAAGGGCATATGGAAGCCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCAGCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTCCAGTAACCCGGAAGGGGGAATGTATCTGGTAGaagcg                                                                                                                              <  2:376735/250‑4 (MQ=255)
                                                                                              aCCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCg                                                                                                                              <  1:17025/250‑1 (MQ=255)
                                                                                                cGGACAGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCg                                                                                                                           >  1:489675/1‑251 (MQ=255)
                                                                                                  gACGGCGCGATCCGTATTTCTTCCGCACCCGGTGTTTCGCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGAGTCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCTCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCAGGCCGGGCATCCGGAAGTGGAAGGTACGATGTGCCAGTACAGTAGCCCGGAAGGGGGAATTTATCTGGTCGAATTGCCGGa                                                                                                                         <  2:640303/251‑1 (MQ=255)
                                                                                                   aCGGCGCGATCCTGATTTTCTCCGCTCACGTTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGGTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                        <  1:524642/251‑1 (MQ=255)
                                                                                                   aCGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGAc                                                                                                                        >  1:511134/1‑251 (MQ=255)
                                                                                                       cgcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtg                                                                                                                    <  1:618366/251‑1 (MQ=255)
                                                                                                        gcgATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGAtgt                                                                                                                   >  1:194208/1‑251 (MQ=255)
                                                                                                          gATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGCAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGGATGTATCTGGTCGAATCGCCGGACGAtgtgt                                                                                                                 >  1:567684/1‑251 (MQ=255)
                                                                                                                 atttctgccGCACACGGTGTTTCTCAGGCGGTCCGTAACGAAGCAAAAAGCCGAGATTTGACGGTGTTTGATGCCCCCTGTCCGCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGGATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAAca                                                                                                          <  2:33752/244‑2 (MQ=255)
                                                                                                                  atttatgccgactaggttgttCTTAGGCGGGACGTAAAGAAGCAAAAAGCAGCGATTTGACGGGGTTTTATGCAACCGGCCTGCTGGTGCCCAAAGTGCATATGGAAGGCGCCCGCGCCAGCTT‑GCAG‑GGCGAAGAGTCGATTCCCAACGGACACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCTGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTGGAAACTg                                                                                                         <  2:126036/232‑1 (MQ=255)
                                                                                                                   tttCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCTCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTg                                                                                                         <  1:16671/250‑1 (MQ=255)
                                                                                                                         cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTc                                                                                                   >  1:589189/1‑250 (MQ=255)
                                                                                                                         cGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCa                                                                                                  <  1:324194/251‑1 (MQ=255)
                                                                                                                           cacaCGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGTAGTGGTAGGTACGATGGGCCAGTACAGTAACCCGGGAGGGGGAATGTATCTGGTCGAATCGCCGGACGAGTTGGGGGAACTAACGGTCaa                                                                                                 >  2:253061/1‑250 (MQ=255)
                                                                                                                              acGGTGGTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGAGTTGACGGTGTGTGAGGCCCCCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAc                                                                                             <  2:58327/251‑1 (MQ=255)
                                                                                                                                  tgtTTCTCAGGCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGa                                                                                           >  1:78523/1‑249 (MQ=255)
                                                                                                                                          aGGCGGTACGTAACGAAGCAAAAAGACGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGACGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGACGAGAAGctc                                                                                 >  1:483803/1‑251 (MQ=255)
                                                                                                                                           ggCGGTACGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAATTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGctct                                                                                <  1:382809/251‑1 (MQ=255)
                                                                                                                                             gtgtatgTAACCAAGCAAAAAGCCGTGAGTTGAAGTTGTTTGAGGCACCCTGTCCGCTGGTGCCCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGTAGTTACGATGGGCCAGTACAGGAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGGGTGGAAACTGACGGTCAAAATCGAAGAGAAGctctac                                                                              <  2:479171/245‑3 (MQ=255)
                                                                                                                                                 accgaaCGAAGCAAAAAGCCGTGATTTGACGGTTTTTGATGCCGCCTTGCCGCTGGTGTCCAAATTGCATTTGGAAGTCGCCCGCGCCAG‑CC‑GTCGGGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCttt                                                                           <  2:430311/246‑1 (MQ=255)
                                                                                                                                                 aCGTAACGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTa                                                                          <  1:146492/251‑1 (MQ=255)
                                                                                                                                                     aaCGAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGAc                                                                      >  1:397084/1‑251 (MQ=255)
                                                                                                                                                        gAAGCAAAAAGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                                                                                >  1:392346/1‑125 (MQ=255)
                                                                                                                                                        gAAGCAAAAAGCCGCGATTTGACGGGGTTTGATGCCACCTTTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCTAAGAGTCGATTCTCATCGGCCACGCCGGGCat                                                                                                                                                                                                <  2:392346/126‑2 (MQ=255)
                                                                                                                                                              aaaaGCCGCGATTGGACTGTGTTTGTTGCTACCTGGCCGCTGGTGACCAAAGTGCTTATGGAAGTCGCCCGCGCCAG‑CC‑GGCGTGGCGAAGAGGCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTATCCCGGAAGGGGGAATGTATCTGGTAGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCAc                                                             <  2:356409/251‑1 (MQ=255)
                                                                                                                                                               aaaGCCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCAACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGAcaa                                                              >  1:389067/1‑247 (MQ=255)
                                                                                                                                                                       gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTc                                                       <  1:103076/248‑1 (MQ=255)
                                                                                                                                                                       gATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCCCTCCTTTCTGACCCAGACCAcactttc                                                       >  2:103076/1‑242 (MQ=255)
                                                                                                                                                                            gACGGTGCTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGACGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGCAGGTACGATGGGCCAGTACCGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCAGGACGATGCGGGGAAACTGACGGTCAAAAACGAAGAAAAGCTCTCCTTTATGACCCAGACCACGCTTTCGGTGGATg                                               >  2:252639/1‑251 (MQ=255)
                                                                                                                                                                              cGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGac                                             >  1:645707/1‑251 (MQ=255)
                                                                                                                                                                                gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                      <  2:158486/144‑1 (MQ=255)
                                                                                                                                                                                gtgtTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGaa                                                                                                                                                      >  1:158486/1‑144 (MQ=255)
                                                                                                                                                                                     ttaTGCCCCCTGTCCGGTGGTGACCAAATTGATTATGTAAGCCGCCCCCGCC‑G‑TCCGTCGTGGCGAAGAGCCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAATGTGCGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGTCGATGTGGGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACTCGTCt                                       <  2:122761/248‑1 (MQ=255)
                                                                                                                                                                                      gATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGa                                     <  1:480145/251‑1 (MQ=255)
                                                                                                                                                                                                    gCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAag                                                                                         >  1:126007/1‑185 (MQ=255)
                                                                                                                                                                                                    gCTGGTGACCAAAGTGAATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGACACGCCGGGCATCCGGACGTGGAAGGAACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAAAGAAGAGAAGCTCTCCTTTATGAACCAGACCACGCTGTCGGTGGAAGACACGtaggttgtgaacgccgggc                       >  1:669822/1‑233 (MQ=255)
                                                                                                                                                                                                               aaGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGCAGTGGAAGGTACGCTGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGGGTGGAAACTGACGGTCAAAAACGAAGAGACGCTCTCCTTTTTGACCCAGAACACGCTGTCGGTGGATGACACGTCTGATGTGATAGACGCGCTGCGGAAACGc            >  2:257378/1‑251 (MQ=255)
                                                                                                                                                                                                                aGTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCt           <  1:38673/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTGCATATGGAAGTCGCCCGCGCCAG‑CC‑GTCGTGGCGAAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGGTGGGCCAGTACAGGAACCCGGAAGGGGGAATGGATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCATTTATGAACCAGACCACGCGATCGGAGGGTGACCCTTCTGGTGTGATCGAAGGGCGGCGTAGACGCtt          >  2:465199/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggAAGTCGGCCCCGCCAG‑CC‑GTCGGTGCGGAGAGTCGATTCTCATCGGCCACGCCGGGCATCCGGAAGTGGAAGGTACGTTGGGGCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCTGACGATGTGTTGAAACTGACGGTCAAAAACGAAGAGAATCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGaaaa  <  2:558736/251‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                         
CGATATATGTCCGTCACGAAGTGGTACATAACCGCTATGTGGTCGATAGCTTGCGTGAGCGTGGGGCTATCTTTATTGAGCAGATTAGCGAAGTACCGGACGGCGCGATCCTGATTTTCTCCGCACACGGTGTTTCTCAGGCGGTACGTAACGAAGCAAAAAGTCGCGATTTGACGGTGTTTGATGCCACCTGTCCGCTGGTGACCAAAGTGCATATGGAAGTCGCCCGCGCCAG‑TC‑GCCGTGGCGAAGAATCTATTCTCATCGGTCACGCCGGGCACCCGGAAGTGGAAGGGACAATGGGCCAGTACAGTAACCCGGAAGGGGGAATGTATCTGGTCGAATCGCCGGACGATGTGTGGAAACTGACGGTCAAAAACGAAGAGAAGCTCTCCTTTATGACCCAGACCACGCTGTCGGTGGATGACACGTCTGATGTGATCGACGCGCTGCGTAAACGCTTCCCGAAAA  >  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: ‑