Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,590 G→T 10.0% intergenic (‑61/+32) leuA ← / ← leuL 2‑isopropylmalate synthase/leu operon leader peptide

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091383,5900GT10.0% 418.5 / 25.7 126intergenic (‑61/+32)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base G (51/62);  new base T (9/4);  total (60/66)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.86e-01

Only 100 of 127 total aligned reads displayed.

CAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATG  >  NC_000913/83443‑83736
                                                                                                                                                   |                                                                                                                                                   
cAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATgg                                                                                                                                                    >  4:132384/1‑149 (MQ=255)
    ttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgcgc                                                                                                                                                 >  2:751734/1‑148 (MQ=255)
      ttctttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCgg                                                                                                                                               <  1:751734/148‑1 (MQ=255)
          ttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttt                                                                                                                                           >  1:966294/1‑148 (MQ=255)
           tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt                                                                                                                                         >  3:881767/1‑149 (MQ=255)
            cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAATCCGCGCAATGGCGCGGGtttttt                                                                                                                                        <  4:60972/149‑1 (MQ=255)
            cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt                                                                                                                                         >  2:381981/1‑148 (MQ=255)
              cacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttg                                                                                                                                       <  1:443107/148‑1 (MQ=255)
               acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgc                                                                                                                                                   >  3:822006/1‑135 (MQ=255)
               acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgc                                                                                                                                                   <  4:822006/135‑1 (MQ=255)
               acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgtt                                                                                                                                     >  3:326026/1‑149 (MQ=255)
               acTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgt                                                                                                                                      <  1:121705/148‑1 (MQ=255)
                cTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgttt                                                                                                                                    <  3:377842/149‑1 (MQ=255)
                  cAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAATAAACCCGCGCAATGGCGCGGGTTTtttctttg                                                                                                                                   <  1:278628/148‑1 (MQ=255)
                     gCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACt                                                                                                                                >  1:417339/1‑148 (MQ=255)
                      cctGCCTGTAACGCCTTTTCCCCGGCGCGCAATTTGGGATCGAAAATAATGACTTGGTGGGCCATTGTTTGGGTTCTTGTTTCTTTTAGAGCGCCCCGCTTCCGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTg                                                                                                                               <  1:147822/146‑1 (MQ=255)
                        tGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCg                                                                                                                             >  1:614821/1‑148 (MQ=255)
                        tGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGt                                                                                                                            <  4:838506/149‑1 (MQ=255)
                           cTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGGt                                                                                                                                             <  3:556175/129‑1 (MQ=255)
                           cTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGGt                                                                                                                                             >  4:556175/1‑129 (MQ=255)
                               aaCGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGc                                                                                                                     <  1:955103/148‑1 (MQ=255)
                                     tGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAATAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGc                                                                                                              <  3:520152/149‑1 (MQ=255)
                                     tGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATg                                                                                                               <  2:540972/148‑1 (MQ=255)
                                       ttCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTg                                                                                                            <  4:86564/149‑1 (MQ=255)
                                                        tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTc                                                                                           >  3:72550/1‑149 (MQ=255)
                                                        tgGTATCGAAAATAATGACTCGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGGTTCGGGCATAAAAAAACCCGCGCAAGGGCGCAGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATACCGCTCACTCGt                                                                                            <  2:281071/148‑1 (MQ=255)
                                                             tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAGAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcg                                                                                      >  4:1188/1‑149 (MQ=255)
                                                                aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTATCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                                                                    >  1:794947/1‑148 (MQ=255)
                                                                 aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCgcgc                                                                                   >  2:889796/1‑148 (MQ=255)
                                                                  aataatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCa                                                                                  >  1:931134/1‑148 (MQ=255)
                                                                      atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                             <  4:326026/149‑1 (MQ=255)
                                                                      atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                             <  4:881767/149‑1 (MQ=255)
                                                                           ttGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGGGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTTGATGCGTGGCGGGCTTAATGCGGGAGGCGGCTAACTCGCTTACGGGGCAAAGAAGATg                                                                         >  3:666943/1‑149 (MQ=255)
                                                                            tGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCg                                                                       >  1:411235/1‑148 (MQ=255)
                                                                              cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                     <  3:453492/149‑1 (MQ=255)
                                                                              cTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                     <  3:363981/149‑1 (MQ=255)
                                                                                ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                    >  2:593302/1‑148 (MQ=255)
                                                                                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGa                                                                                                          <  3:718807/107‑1 (MQ=255)
                                                                                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGa                                                                                                          >  4:718807/1‑107 (MQ=255)
                                                                                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                                 <  2:580928/148‑1 (MQ=255)
                                                                                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt                                                                 <  2:930424/148‑1 (MQ=255)
                                                                                   tCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑ACTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagta                                                                >  3:815167/1‑149 (MQ=255)
                                                                                    cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                    <  2:955972/144‑1 (MQ=255)
                                                                                    cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGACTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                    >  1:955972/1‑144 (MQ=255)
                                                                                      tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagt                                                              <  1:833523/148‑1 (MQ=255)
                                                                                        gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt                                                           <  3:124838/149‑1 (MQ=255)
                                                                                           tGGGTCCTGGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTGGTCTACCGCGCTAAGAAGATGCGTTTAGTAGTAGTAGa                                                        <  3:593002/149‑1 (MQ=255)
                                                                                            gggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc                                                       <  3:6762/149‑1 (MQ=255)
                                                                                             ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAcc                                                      >  3:233450/1‑149 (MQ=255)
                                                                                                ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                    >  2:598506/1‑148 (MQ=255)
                                                                                                ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                    >  2:350099/1‑148 (MQ=255)
                                                                                                ccTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                    >  2:577098/1‑148 (MQ=255)
                                                                                                  ttGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCg                                                                                                     <  3:448975/97‑1 (MQ=255)
                                                                                                  ttGGCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCg                                                                                                     >  4:448975/1‑97 (MQ=255)
                                                                                                   tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaac                                                <  4:823892/149‑2 (MQ=255)
                                                                                                   tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                 >  2:204984/1‑148 (MQ=255)
                                                                                                   tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                 <  1:889796/148‑1 (MQ=255)
                                                                                                   tGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaa                                                 <  1:414507/148‑1 (MQ=255)
                                                                                                    gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATaaa                                                <  1:950605/148‑1 (MQ=255)
                                                                                                    gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAg                                               <  4:233450/149‑1 (MQ=255)
                                                                                                     tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                     >  1:631942/1‑126 (MQ=255)
                                                                                                     tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                              <  4:339984/149‑1 (MQ=255)
                                                                                                     tcgctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGtt                                                                     <  2:631942/123‑1 (MQ=255)
                                                                                                       tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                            >  3:804803/1‑149 (MQ=255)
                                                                                                       tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGa                                            >  3:261578/1‑149 (MQ=255)
                                                                                                         ttttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                          >  3:343633/1‑149 (MQ=255)
                                                                                                         ttttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAc                                          >  3:478353/1‑149 (MQ=255)
                                                                                                             agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                       <  2:925389/148‑1 (MQ=255)
                                                                                                                gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgtg                                    <  1:204984/148‑1 (MQ=255)
                                                                                                                gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtgtg                                    <  1:577098/148‑1 (MQ=255)
                                                                                                                  gcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAg                                  <  2:966294/148‑1 (MQ=255)
                                                                                                                   ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTc                                >  4:861578/1‑149 (MQ=255)
                                                                                                                     tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                   >  4:470827/1‑144 (MQ=255)
                                                                                                                     tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                   <  3:470827/144‑1 (MQ=255)
                                                                                                                     tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                               >  1:612925/1‑148 (MQ=255)
                                                                                                                     tCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                               >  2:355096/1‑148 (MQ=255)
                                                                                                                      cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                               >  1:694386/1‑147 (MQ=255)
                                                                                                                      cGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCa                               <  2:694386/147‑1 (MQ=255)
                                                                                                                       gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCAtt                             <  2:479788/148‑1 (MQ=255)
                                                                                                                         ttCGGGCATAAAAAAACCCGCGCAATGGCGCGGGCTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGGCATTaa                           >  2:897743/1‑148 (MQ=255)
                                                                                                                             ggCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                      <  3:557470/149‑1 (MQ=255)
                                                                                                                                 tAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc                  >  3:432869/1‑149 (MQ=255)
                                                                                                                                 ataaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAATAAGATGCGTTTAATAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc                  <  4:700346/147‑1 (MQ=255)
                                                                                                                                  aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCa                 >  3:896236/1‑149 (MQ=255)
                                                                                                                                   aaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg                >  4:150583/1‑149 (MQ=255)
                                                                                                                                    aaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGa               <  4:615694/149‑1 (MQ=255)
                                                                                                                                      aaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg             >  3:342727/1‑149 (MQ=255)
                                                                                                                                        aCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGa                                   <  3:601141/125‑1 (MQ=255)
                                                                                                                                        aCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCACAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                         >  1:892632/1‑135 (MQ=255)
                                                                                                                                        aCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCACAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                         <  2:892632/135‑1 (MQ=255)
                                                                                                                                        aCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGa            <  1:714887/148‑1 (MQ=255)
                                                                                                                                         cccGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGaa           >  2:549162/1‑148 (MQ=255)
                                                                                                                                         cccGCGCCATTGCGCGGGTTTTTGTTTTGCCTGCGT‑GCTGGCTTAATCCTGAATCCCGCTCACTGGTCTACCGCGCAAAGAAGATGCGTTAGGTAGTAGTAGACCGATAAAGCGAACGATGGGAGTCATTAAATAAGTTCCTGATGaag          >  3:821826/1‑148 (MQ=255)
                                                                                                                                         cccGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAAt          <  4:261578/149‑1 (MQ=255)
                                                                                                                                         cccGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAAt          <  3:45812/149‑1 (MQ=255)
                                                                                                                                           cgcgCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg         <  2:113859/148‑1 (MQ=255)
                                                                                                                                             cgcAACGGGGCGGGTGTTTTGTTTGACTGCGT‑GCGGGCGTAATGCTGGATGCCGCTCACGCGTCCACCGCGCAAAGAAGATGCGAGTAGTAGGAGGAGACCGATAAAGCGAACGAGGGGAGTCATTAAATCAGCTCCAGATGAATGCg       <  2:561745/148‑1 (MQ=255)
                                                                                                                                               cAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGata    <  4:343633/149‑1 (MQ=255)
                                                                                                                                               cAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGata    <  4:170880/149‑1 (MQ=255)
                                                                                                                                                 aTGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATg  <  4:804803/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGT‑GCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGATATG  >  NC_000913/83443‑83736

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