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

ACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGC  >  NC_000913/3333672‑3334266
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
acaGTAAGGTAATCTGGGCGAAAAGCCCTGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGa                                                                                                                                                                                                                                                                                                        <  1:41223/301‑1 (MQ=255)
  aGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCCACTGATCAATCAGTTAGGCTATTACATCGTCCGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCCCGCGCTGTTTGCTATGAGGGAAATGCGCATGTCCCCATTGCTGCCCTGATCC‑AGTTTATCCA‑CCCCGCCACCCCGCCACACCACCACCGTGTGGaag                                                                                                                                                                                                                                                                                                      >  1:410726/1‑299 (MQ=255)
             cGGGGCGAAAAGCCTGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATg                                                                                                                                                                                                                                                                                           <  2:411289/301‑1 (MQ=255)
                 gCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGAAGTCCAACTGATCAATCAGTTAGGCTATTTCATCGTCCGCGGCGGCGGGAAACGTATTCGTCCGATGATTTCTGTACTGGCTTCCCGAGCTGTTGGCTATGGGGGGAATGCCCATGTCACCCTTTCTGCCCTGGCAGCAGTTTA‑CCA‑CCCCGCGACTCCACTACACCAAGACGCTTTGGAAGCAACAAATATgcacc                                                                                                                                                                                                                                                                                       >  1:439853/1‑298 (MQ=255)
                    aaatgaacGGATTTTGCGTTGAATTGAGAATAAATCTATGAGTTTGTGGCGCAAGAAATTGC‑TGTGTGTAAGGCGTCAAGCCGTGAGCAGGTTCATTGCGAGGTCCAACGTAGCAATCAGGTAGGGTAGTACAGCGTCAGCGGCGGCGGTAAACGGATGAGTCCGATGATGGCTGTACTGGGTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATTCA‑CACGGCGACTCGGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAggg                                                                                                                                                                                                                                                                                    <  1:339400/294‑1 (MQ=255)
                    aaaaGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAAATGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACATATTCGTCCGATGATTGCTGTACTGGCTGCACGCGCTGTTGGCTATGAGGGAAATGCGAATGTCAACATTGCTGCCCTGCACG‑AGTTTATCCACCACGCCAAC‑CTGCTACACCACGCCAGTGTGGCAGAAACAGaaaagcccacgg                                                                                                                                                                                                                                                                                    >  1:518338/1‑290 (MQ=255)
                                     gATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGAATCTGCTACACGACGACGTTGAGGATGAAACAGATAAGCGCAGGGGTAAAACTAACGCaaaa                                                                                                                                                                                                                                                                   >  2:467734/1‑300 (MQ=255)
                                                                        aaGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAgg                                                                                                                                                                                                                                <  2:519010/301‑1 (MQ=255)
                                                                             tgggggggTTT‑TAATGCGGCAAGCCTTGAGCAGCTTATTTCGGACGTCCAACTGTTCATTCAGTTAGTCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGGTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATgcctaccagcacgctggcggcattg                                                                                                                                                                                                                           <  1:467170/296‑23 (MQ=255)
                                                                                  gggTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATtttatt                                                                                                                                                                                                                      <  1:45505/301‑1 (MQ=255)
                                                                                   ggTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTAGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATtttattt                                                                                                                                                                                                                     >  2:227881/1‑301 (MQ=255)
                                                                                           ttGCGGCAATCCTTGAGCAGCTTAATTACGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGGCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCcg                                                                                                                                                                                                              <  2:117991/300‑1 (MQ=255)
                                                                                                     ttccgatctgcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAg                                                                                                                                                                                                    <  2:271275/292‑1 (MQ=255)
                                                                                                              gcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAgagatccgaa                                                                                                                                                                                           >  1:270975/1‑293 (MQ=255)
                                                                                                              gcTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGAcc                                                                                                                                                                                           >  2:529007/1‑301 (MQ=255)
                                                                                                                       cGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGAATTTATTTATACCCGCGCTTTCCAGATGATAACCCGccccagt                                                                                                                                                                                  >  1:234697/1‑296 (MQ=255)
                                                                                                                             ccAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTCTTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCCGGGGTAAAGCTACCGCCCACGCCGCCTTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATCCCCGCGCTTTACAGCTGAAGACCAGCCTCGGGTCGCtt                                                                                                                                                                            >  1:470501/1‑300 (MQ=255)
                                                                                                                             caaattGATCTATCAGTTAGGTTATTAAATCGTTATCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTc                                                                                                                                                                            <  1:638845/296‑1 (MQ=255)
                                                                                                                                      caatcaGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTg                                                                                                                                                                   >  2:440514/1‑301 (MQ=255)
                                                                                                                                      caatcaGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGATGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGCAAAGCTACCGCCAACGACGCATTTGGCAATGCCGCCAGCCGCCTGGTAAGCAATTTTATTTATACCCGCCCTTTCCAGCTAACGACCAGCCACGCTACGACCAAAgaccgt                                                                                                                                                                   >  1:571426/1‑296 (MQ=255)
                                                                                                                                           aGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGt                                                                                                                                                              <  2:440417/301‑1 (MQ=255)
                                                                                                                                             ttgtgtgtttACATGGGCAGCGGCGGCGGTAAACGTATTCGTCCGTTGATTGCTGTACTGTCTGCACGAGCTGTTGTCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCa                                                                                                                                                            <  1:439950/293‑1 (MQ=255)
                                                                                                                                             tgcttttccgatctCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCTCCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCa                                                                                                                                                            <  2:239489/288‑1 (MQ=255)
                                                                                                                                                           cGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCTCCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAATCaagatcggaaaagcg                                                                                                                                              >  1:239222/1‑287 (MQ=255)
                                                                                                                                                              cAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGAGATTTTATTTATACCCGAGCTTTCAAGATGCTGACCAGCCTCAGTTTCCTCAAAGTGCTGGCAATCATTTCAGAAGCGGTATAc                                                                                                                                           >  1:329950/1‑301 (MQ=255)
                                                                                                                                                                      gcggTAAACGGTTTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTGTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGAGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGc                                                                                                                                   <  2:471102/301‑1 (MQ=255)
                                                                                                                                                                       cggTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCa                                                                                                                                  >  2:465435/1‑301 (MQ=255)
                                                                                                                                                                            aaCGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAgg                                                                                                                             >  1:58604/1‑301 (MQ=255)
                                                                                                                                                                                      cctggtggtggCTGTACTGGCTGTACGAGCTGTTGGCTATGAGTGAAATGCGCATGGCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCGGCTACACGACGACGGTGTGGATGAATCAGATATGCGCAGGGGTAAAGCGACGGCCAACGCCGCATTTGGCAAGGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCGCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTg                                                                                                                   <  1:230881/291‑1 (MQ=255)
                                                                                                                                                                                         gatgatTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGAGACGTTCTGCaa                                                                                                                >  1:382882/1‑301 (MQ=255)
                                                                                                                                                                                                   ttACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACg                                                                                                      <  2:58672/300‑1 (MQ=255)
                                                                                                                                                                                                        ggCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTaaa                                                                                                 >  2:36885/1‑300 (MQ=255)
                                                                                                                                                                                                         gcgctattccgatcttgCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACg                                                                                                <  1:93389/285‑1 (MQ=255)
                                                                                                                                                                                                               cGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTAAACGATACga                                                                                          >  2:502167/1‑300 (MQ=255)
                                                                                                                                                                                                                       gggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGAGCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACt                                                                                  <  1:528334/300‑1 (MQ=255)
                                                                                                                                                                                                                        ggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGagatcggaagagcaa                                                                                 >  2:93513/1‑287 (MQ=255)
                                                                                                                                                                                                                        ggCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGAACAGGACATCACta                                                                                 >  2:370102/1‑300 (MQ=255)
                                                                                                                                                                                                                           tgtGAGGGTAATGCTGATGTCACCATTGTTGCTCTGTGCG‑AGTTTATCCA‑CACGGTGATGCTGCTACACGACGTCTTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaag                                                                              <  1:369606/299‑1 (MQ=255)
                                                                                                                                                                                                                           tATGAGGGAAATGCGCATGTTACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaag                                                                              <  2:330355/301‑1 (MQ=255)
                                                                                                                                                                                                                           tATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTgaaa                                                                              >  2:499561/1‑300 (MQ=255)
                                                                                                                                                                                                                                  gAAATGTGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGCAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGCAGAAAACtc                                                                       >  1:494694/1‑300 (MQ=255)
                                                                                                                                                                                                                                  gAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTaagatcggaagagcaaacgtctgaactcc                                                                       >  2:261842/1‑274 (MQ=255)
                                                                                                                                                                                                                                   aaaTGCGGATGTCACCATTTCTGCCCTGGTCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGAAGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTAc                                                                      <  1:227638/301‑1 (MQ=255)
                                                                                                                                                                                                                                             ttCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTa                                                            <  2:495341/300‑1 (MQ=255)
                                                                                                                                                                                                                                              tCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGAAATCACTGAAGAAAACTAAATgagagtaat                                                           >  2:332887/1‑293 (MQ=255)
                                                                                                                                                                                                                                                cccAGTGCTGCCCTGGTCA‑TGTTTTTCCA‑CTCGGCGACTTTGCTACACGTCTACGTTGTGGATGAATCAGATATGCGCAGGTGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTATGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCt                                                         <  1:332482/300‑1 (MQ=255)
                                                                                                                                                                                                                                                  ctttgtggcCCTGTTCG‑ATTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtat                                                       <  2:572140/294‑1 (MQ=255)
                                                                                                                                                                                                                                                    tttctgcCCGGATCG‑AGTTTATCCA‑CACGGCGACGCTGCTACACGACGACGTTGTGGATGAATCAGAGATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAg                                                     <  1:498914/298‑1 (MQ=255)
                                                                                                                                                                                                                                                    gacgttctttctatcttgTTTATCCT‑CACTGCGACTCTGCTACACGTCGACGTGGTGGATGAATCAGTTATGCGGAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtata                                                      <  1:626823/284‑1 (MQ=255)
                                                                                                                                                                                                                                                       ctgcCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGACAACTACATGCGCGTTATCTATAGCaa                                                  >  2:308902/1‑301 (MQ=255)
                                                                                                                                                                                                                                                          cccTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTAAATGCACGTAACCTAAAGCAAAAc                                               >  1:409486/1‑301 (MQ=255)
                                                                                                                                                                                                                                                               aTCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCCGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACAGAAAAAAACTACAAACGCGATaaccaaagacaaaacaccc                                          >  1:285440/1‑283 (MQ=255)
                                                                                                                                                                                                                                                               aTCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCAGTAAACGTAATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATACGGACAACACTGAAGAAAACTAACTCCACGTTATCTAAAGACAAACcgcgc                                          >  1:173347/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                 cG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAACCTACATGCGCGTTATCTATAGCAAAACCGAGCGt                                        >  2:590831/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                    aGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCtataagatcggaagagcaca                                      >  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)
                                                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                                                                         
ACAGTAAGGTAATCGGGGCGAAAAGCCCGGCTTTTGCGATGAATTTAGAAAAAATCAATGAGTTAACCGCGCAAGATATGGCGGGTGT‑TAATGCGGCAATCCTTGAGCAGCTTAATTCCGACGTCCAACTGATCAATCAGTTAGGCTATTACATCGTCAGCGGCGGCGGTAAACGTATTCGTCCGATGATTGCTGTACTGGCTGCACGAGCTGTTGGCTATGAGGGAAATGCGCATGTCACCATTGCTGCCCTGATCG‑AGTTTATCCA‑CACGGCGACTCTGCTACACGACGACGTTGTGGATGAATCAGATATGCGCAGGGGTAAAGCTACCGCCAACGCCGCATTTGGCAATGCCGCCAGCGTGCTGGTAGGCGATTTTATTTATACCCGCGCTTTCCAGATGATGACCAGCCTCGGTTCGCTCAAAGTGCTGGAAGTCATGTCAGAAGCCGTAAACGTCATCGCAGAAGGTGAAGTTCTGCAACTGATGAACGTTAACGATCCGGACATCACTGAAGAAAACTACATGCGCGTTATCTATAGCAAAACCGCGCGTCTGTTTGAGGCTGCCGCGCAGTGTTCCGGGATTCTGGC  >  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: ‑