Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,333,969 T→A 6.6% V87E (GTG→GAG)  ispB → octaprenyl diphosphate synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,333,9690TA 6.6% 184.1 / 2.6 60V87E (GTG→GAG) ispBoctaprenyl diphosphate synthase
Reads supporting (aligned to +/- strand):  ref base T (30/26);  new base A (2/2);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

ACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGC  >  NC_000913/3333672‑3334266
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
acaGTAAGGTAATCTGGGCGAAAAGCCCTGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGa                                                                                                                                                                                                                                                                                                        <  1:41223/301‑1 (MQ=255)
  aGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCCACTGATCAATCAGTTAGGCTATTACATCGTCCGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCCCGCGCTGTTTGCTATGAGGGAAATGCGCATGTCCCCATTGCTGCCCTGATCCAGTTTATCCACCCCGCCACCCCGCCACACCACCACCGTGTGGaag                                                                                                                                                                                                                                                                                                      >  1:410726/1‑299 (MQ=255)
             cGGGGCGAAAAGCCTGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATg                                                                                                                                                                                                                                                                                           <  2:411289/301‑1 (MQ=255)
                 gCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGAAGTCCAACTGATCAATCAGTTAGGCTATTTCATCGTCCGCGGCGGCGGGAAACGTATTCGTCCGATGATTTCTGTACTGGCTTCCCGAGCTGTTGGCTATGGGGGGAATGCCCATGTCACCCTTTCTGCCCTGGCAGCAGTTTACCACCCCGCGACTCCACTACACCAAGACGCTTTGGAAGCAACAAATATgcacc                                                                                                                                                                                                                                                                                       >  1:439853/1‑298 (MQ=255)
                    aaatgaacGGATTTTGCGTTGAATTGAGAATAAATCTATGAGTTTGTGGCGCAAGAAATTGCTGTGTGTAAGGCGTCAAGCCGTGAGCAGGTTCATTGCGAGGTCCAACGTAGCAATCAGGTAGGGTAGTACAGCGTCAGCGGCGGCGGTAAACGGATGAGTCCGATGATGGCTGTACTGGGTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATTCACACGGCGACTCGGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAggg                                                                                                                                                                                                                                                                                    <  1:339400/294‑1 (MQ=255)
                    aaaaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAAATGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACATATTCGTCCGATGATTGCTGTACTGGCTGCACGCGCTGTTGGCTATGAGGGAAATGCGAATGTCAACATTGCTGCCCTGCACGAGTTTATCCACCACGCCAACCTGCTACACCACGCCAGTGTGGCAGAAACAGaaaagcccacgg                                                                                                                                                                                                                                                                                    >  1:518338/1‑290 (MQ=255)
                                     gATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGAATCTGCTACACGACGACGTTGAGGATGAAACAGATAAGCGCAGGGGTAAAACTAACGCaaaa                                                                                                                                                                                                                                                                   >  2:467734/1‑300 (MQ=255)
                                                                        aaGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAgg                                                                                                                                                                                                                                <  2:519010/301‑1 (MQ=255)
                                                                             tgggggggTTTTAATGCGGCAAGCCTTGAGCAGCTTATTTCGGACGTCCAACTGTTCATTCAGTTAGTCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGGTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATgcctaccagcacgctggcggcattg                                                                                                                                                                                                                           <  1:467170/296‑23 (MQ=255)
                                                                                  gggTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATtttatt                                                                                                                                                                                                                      <  1:45505/301‑1 (MQ=255)
                                                                                   ggTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTAGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATtttattt                                                                                                                                                                                                                     >  2:227881/1‑301 (MQ=255)
                                                                                           ttGCGGCAATCCTTGAGCAGCTTAATTACGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGGCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCcg                                                                                                                                                                                                              <  2:117991/300‑1 (MQ=255)
                                                                                                     ttccgatctgcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAg                                                                                                                                                                                                    <  2:271275/292‑1 (MQ=255)
                                                                                                              gcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAgagatccgaa                                                                                                                                                                                           >  1:270975/1‑293 (MQ=255)
                                                                                                              gcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGAcc                                                                                                                                                                                           >  2:529007/1‑301 (MQ=255)
                                                                                                                       cGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGAATTTATTTATACCCGCGCTTTCCAGATGATAACCCGccccagt                                                                                                                                                                                  >  1:234697/1‑296 (MQ=255)
                                                                                                                             ccAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTCTTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCCGGGGTAAAGCTACCGCCCACGCCGCCTTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATCCCCGCGCTTTACAGCTGAAGACCAGCCTCGGGTCGCtt                                                                                                                                                                            >  1:470501/1‑300 (MQ=255)
                                                                                                                             caaattGATCTATCAGTTAGGTTATTAAATCGTTATCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTc                                                                                                                                                                            <  1:638845/296‑1 (MQ=255)
                                                                                                                                      caatcaGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTg                                                                                                                                                                   >  2:440514/1‑301 (MQ=255)
                                                                                                                                      caatcaGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGATGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGCAAAGCTACCGCCAACGACGCATTTGGCAATGCCGCCAGCCGCCTGGTAAGCAATTTTATTTATACCCGCCCTTTCCAGCTAACGACCAGCCACGCTACGACCAAAgaccgt                                                                                                                                                                   >  1:571426/1‑296 (MQ=255)
                                                                                                                                           aGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGt                                                                                                                                                              <  2:440417/301‑1 (MQ=255)
                                                                                                                                             ttgtgtgtttACATGGGCAGCGGCGGCGGTAAACGTATTCGTCCGTTGATTGCTGTACTGTCTGCACGAGCTGTTGTCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCa                                                                                                                                                            <  1:439950/293‑1 (MQ=255)
                                                                                                                                             tgcttttccgatctCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCTCCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCa                                                                                                                                                            <  2:239489/288‑1 (MQ=255)
                                                                                                                                                           cGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCTCCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAATCaagatcggaaaagcg                                                                                                                                              >  1:239222/1‑287 (MQ=255)
                                                                                                                                                              cAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGAGATTTTATTTATACCCGAGCTTTCAAGATGCTGACCAGCCTCAGTTTCCTCAAAGTGCTGGCAATCATTTCAGAAGCGGTATAc                                                                                                                                           >  1:329950/1‑301 (MQ=255)
                                                                                                                                                                      gcggTAAACGGTTTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTGTATCCACACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGc                                                                                                                                   <  2:471102/301‑1 (MQ=255)
                                                                                                                                                                       cggTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCa                                                                                                                                  >  2:465435/1‑301 (MQ=255)
                                                                                                                                                                            aaCGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAgg                                                                                                                             >  1:58604/1‑301 (MQ=255)
                                                                                                                                                                                      cctggtggtggCTGTACTGGCTGTACGAGCTGTTGGCTATGAGTGAAATGCGCATGGCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCGGCTACACGACGACGGTGTGGATGAATCAGATATGCGCAGGGGTAAAGCGACGGCCAACGCCGCATTTGGCAAGGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCGCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTg                                                                                                                   <  1:230881/291‑1 (MQ=255)
                                                                                                                                                                                         gatgatTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGAGACGTTCTGCaa                                                                                                                >  1:382882/1‑301 (MQ=255)
                                                                                                                                                                                                   ttACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACg                                                                                                      <  2:58672/300‑1 (MQ=255)
                                                                                                                                                                                                        ggCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTaaa                                                                                                 >  2:36885/1‑300 (MQ=255)
                                                                                                                                                                                                         gcgctattccgatcttgCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACg                                                                                                <  1:93389/285‑1 (MQ=255)
                                                                                                                                                                                                               cGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTAAACGATACga                                                                                          >  2:502167/1‑300 (MQ=255)
                                                                                                                                                                                                                       gggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGAGCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACt                                                                                  <  1:528334/300‑1 (MQ=255)
                                                                                                                                                                                                                        ggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGagatcggaagagcaa                                                                                 >  2:93513/1‑287 (MQ=255)
                                                                                                                                                                                                                        ggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGAACAGGACATCACta                                                                                 >  2:370102/1‑300 (MQ=255)
                                                                                                                                                                                                                           tgtGAGGGTAATGCTGATGTCACCATTGTTGCTCTGTGCGAGTTTATCCACACGGTGATGCTGCTACACGACGTCTTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaag                                                                              <  1:369606/299‑1 (MQ=255)
                                                                                                                                                                                                                           tATGAGGGAAATGCGCATGTTACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaag                                                                              <  2:330355/301‑1 (MQ=255)
                                                                                                                                                                                                                           tATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaaa                                                                              >  2:499561/1‑300 (MQ=255)
                                                                                                                                                                                                                                  gAAATGTGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGCAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGCAGAAAACtc                                                                       >  1:494694/1‑300 (MQ=255)
                                                                                                                                                                                                                                  gAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTaagatcggaagagcaaacgtctgaactcc                                                                       >  2:261842/1‑274 (MQ=255)
                                                                                                                                                                                                                                   aaaTGCGGATGTCACCATTTCTGCCCTGGTCGAGTTTATCCACACGGCGACTCTGCTACACGACGAAGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTAc                                                                      <  1:227638/301‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTa                                                            <  2:495341/300‑1 (MQ=255)
                                                                                                                                                                                                                                              tCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGAAATCACTGAAGAAAACTAAATgagagtaat                                                           >  2:332887/1‑293 (MQ=255)
                                                                                                                                                                                                                                                cccAGTGCTGCCCTGGTCATGTTTTTCCACTCGGCGACTTTGCTACACGTCTACGTTGTGGATGAATCAGATATGCGCAGGTGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTATGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCt                                                         <  1:332482/300‑1 (MQ=255)
                                                                                                                                                                                                                                                  ctttgtggcCCTGTTCGATTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtat                                                       <  2:572140/294‑1 (MQ=255)
                                                                                                                                                                                                                                                    tttctgcCCGGATCGAGTTTATCCACACGGCGACGCTGCTACACGACGACGTTGTGGATGAATCAGAGATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAg                                                     <  1:498914/298‑1 (MQ=255)
                                                                                                                                                                                                                                                    gacgttctttctatcttgTTTATCCTCACTGCGACTCTGCTACACGTCGACGTGGTGGATGAATCAGTTATGCGGAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtata                                                      <  1:626823/284‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctgcCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGACAACTACATGCGCGTTATCTATAGCaa                                                  >  2:308902/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          cccTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTAAATGCACGTAACCTAAAGCAAAAc                                               >  1:409486/1‑301 (MQ=255)
                                                                                                                                                                                                                                                               aTCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCCGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACAGAAAAAAACTACAAACGCGATaaccaaagacaaaacaccc                                          >  1:285440/1‑283 (MQ=255)
                                                                                                                                                                                                                                                               aTCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCAGTAAACGTAATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATACGGACAACACTGAAGAAAACTAACTCCACGTTATCTAAAGACAAACcgcgc                                          >  1:173347/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 cGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAACCTACATGCGCGTTATCTATAGCAAAACCGAGCGt                                        >  2:590831/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    aGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtataagatcggaagagcaca                                      >  2:627700/1‑285 (MQ=255)
                                                                                                                                                                                                                                                                             cactcgtctaCTCTGCTAAACGACGACTTTGTGTATGAATCATATATGCGCAGGGGTAAATCTTCCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGTAAGTCATGTCATAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAgg                              <  2:173544/292‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          ggCTACACGGCGACGTTGTGGATTAATCAGATATGCGCAGGGGTAATGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAgtg                 <  1:464871/300‑1 (MQ=255)
                                                                                                                                                                                                                                                                                           gCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAgtgt                >  2:360926/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                acgacgacgTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCgg           >  2:258833/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                         ttGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTctggc  <  2:410049/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
ACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGTTAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCGAGTTTATCCACACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGC  >  NC_000913/3333672‑3334266

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