Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,587 A→C 5.7% intergenic (‑58/+35) 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,5870AC 5.7% 189.9 / 2.3 74intergenic (‑58/+35)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (25/43);  new base C (3/2);  total (29/45)
Fisher's exact test for biased strand distribution p-value = 3.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.43e-01

GCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGA  >  NC_000913/83442‑83732
                                                                                                                                                 |                                                                                                                                                  
gtaGTTTGTCTTACACACTGAGGCTTGCCTGTAACGCCTGTTCCCCGTCGCGCAAGGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTGGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                                                                 <  1:210197/147‑1 (MQ=255)
   gTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgcg                                                                                                                                              <  1:689408/149‑1 (MQ=255)
    tttttCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgcgc                                                                                                                                             >  1:1086706/1‑149 (MQ=255)
          tttcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCGCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATCGAAAAACCCGCGCCATGGCGCGGGttt                                                                                                                                       >  2:521941/1‑149 (MQ=255)
           ttcACACTCAAGCTTGCCTGTAACGCCTGCTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttt                                                                                                                                      <  2:954735/149‑1 (MQ=255)
            tcACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGttttt                                                                                                                                     >  2:248081/1‑149 (MQ=255)
             cccccTCAAGCTTGCCTGTAACGCCTGTTCCCAGTCGCGCAATGTGGTATCGCAGATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAGAAAAACCCGCGCAATGGCGCGGGtttttt                                                                                                                                    <  2:634851/145‑1 (MQ=255)
             cacaCTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGtttttt                                                                                                                                    <  2:201738/149‑1 (MQ=255)
              cccccccaGCTTACCTGTAACGCCTGTTCACCGCCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATTGCGCGGGTTCtttg                                                                                                                                   <  2:202253/142‑1 (MQ=255)
              acacTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttg                                                                                                                                   >  2:366929/1‑149 (MQ=255)
              acacTCAAGCTTGCCTGTAACGCCGGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGGttttttt                                                                                                                                   <  2:283765/149‑2 (MQ=255)
               cccTCAAGCTTGCCTGTAACGCCTGTTCCCCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgt                                                                                                                                  <  1:566194/147‑1 (MQ=255)
               cacTCAAGCTTGCCTGTAACGCCTGTTCACCGCCGCGCAATGTGGGATCGAAAATAATGACTTGCTGGCTTATGGCTTGGGTCCTTGACTCTTATAGACCGCCCCGCTTCGGGCACAAAAAAACCCGCGCACTGGCGCGGGTTTtttgt                                                                                                                                  >  2:125727/1‑149 (MQ=255)
                  tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTtttgtttg                                                                                                                               >  1:568035/1‑149 (MQ=255)
                  tCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTCCTTGCCGCTGGCTTGGGTTCTTGACTCTTTGAGATGTCCTCGCATCGGGCTAAAAAAAACCCGGCCAATGCGGCGGGGTCtttgtttg                                                                                                                               >  2:389298/1‑149 (MQ=255)
                   gcaGCTTGCCGGTAACGCCTGTTCACCGGCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGa                                                                                                                              <  2:846800/147‑1 (MQ=255)
                       cTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGc                                                                                                                          <  2:160497/149‑1 (MQ=255)
                            cTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTTTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCt                                                                                                                     >  1:399242/1‑149 (MQ=255)
                            cTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCt                                                                                                                     >  1:962474/1‑149 (MQ=255)
                             tGTAACGCCTGTTCACCGTCGCGCCATGTGGTATCGAAAATAAGGAGTTGCTGGGTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTg                                                                                                                    >  2:104916/1‑149 (MQ=255)
                             tGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTg                                                                                                                    >  2:280342/1‑149 (MQ=255)
                                   gCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGCCTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTaa                                                                                                              <  2:920949/149‑1 (MQ=255)
                                   gCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTaa                                                                                                              >  2:688730/1‑149 (MQ=255)
                                        ttCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATgtgg                                                                                                         >  1:819335/1‑146 (MQ=255)
                                                 gtgcAATGCGGTACCGAAAATAATGACTTGCCGGCTCATGGTTTCGGTCCTTGCCTCTTTTAGAGCGCCTCGCGTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCt                                                                                                <  2:740496/147‑1 (MQ=255)
                                                 gcgcAATGTGGTACCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCt                                                                                                <  1:1120802/149‑1 (MQ=255)
                                                  cgcAATGTGGTGTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTGCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGAGTGCGTGCTGGCTTAATGCTGGATGcaggtc                                                                                               >  2:280449/1‑144 (MQ=255)
                                                    cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                                                                             <  2:1013980/149‑1 (MQ=255)
                                                         tgGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTc                                                                                        <  1:511588/149‑1 (MQ=255)
                                                           gTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTTTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTa                                                                                      <  2:399242/149‑1 (MQ=255)
                                                            tATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAACCCCGCGGGATGGCGCGGGTTTTTTGTTTGACTGCGTCCTGGCTTAATGCAGTAGGCAGCTCACCCGTCTat                                                                                     >  1:24464/1‑148 (MQ=255)
                                                                gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCATTCCCGCTGTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCCCTCACTCGACTACcgcg                                                                                 >  1:561887/1‑149 (MQ=255)
                                                                 aaaaTAACGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCgcgc                                                                                >  2:464046/1‑149 (MQ=255)
                                                                     taatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaag                                                                            <  1:751407/149‑1 (MQ=255)
                                                                      aatGACTGGCGGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCCCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                                                                          <  1:493089/149‑1 (MQ=255)
                                                                       ctcaCTTCCGGGCTCCTGGTTTGCTGCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAATCCCGCGCCATGGTGCGGGTTTTGTGTTTGACTGCGTGCTGGCTGAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                                                                         <  2:214323/146‑1 (MQ=255)
                                                                           cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATg                                                                     <  1:626075/149‑1 (MQ=255)
                                                                              gCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                                                                  <  1:280579/149‑1 (MQ=255)
                                                                                tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                                                                >  1:533477/1‑149 (MQ=255)
                                                                                 ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                               <  1:518278/149‑1 (MQ=255)
                                                                                 ggCTCATGGTTTGGGTCCTTGTCCCGTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTta                                                               <  1:521941/149‑1 (MQ=255)
                                                                                  gCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                                                              >  2:409511/1‑149 (MQ=255)
                                                                                     cATGGTTGGCGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag                                                           <  1:156155/149‑1 (MQ=255)
                                                                                     cATGGTCTGGGTCCTTGTCTCTTTTAGAGCGCTCCGCTTCGGGCATAAAAAACCACGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGATGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag                                                           <  2:637625/149‑1 (MQ=255)
                                                                                       tgttttGGGTCCTTTTCTCTTTTAGGGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                         <  2:189258/146‑1 (MQ=255)
                                                                                       tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta                                                         <  2:1086706/149‑1 (MQ=255)
                                                                                        ggTTTGGGTCCCTGTCTCTTTTAGGGCGCCTCGCTTCGGGCATAAAACAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag                                                        <  2:846027/149‑1 (MQ=255)
                                                                                          tttGGGTCTTTGCCTCTTTTAGAGCGCCTCTCTTTCGCCATCAACACACCCGCGCAATGGCGCGGGTTTTTTTTTTGACTGCGTGTTGGCTTAATGCTGGATGCCGCTCACGCGTCTTCCGCGCAAAGAAGTTGCGTTtagtagtagta                                                      <  2:544026/149‑1 (MQ=255)
                                                                                            tGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGa                                                    <  2:453259/149‑1 (MQ=255)
                                                                                            tGGCTCCTTTCCTCTTTTAGGCCTCCTCGCTTCCGGCATAAAAAACCCCGAGCAATGGTGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTTATGCTGGATGCCGCTCACTCGTCTAACGCGCAAAGAAGATGCGTTTAGTAGTAGTAGa                                                    <  2:1007632/149‑1 (MQ=255)
                                                                                                tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTGGACCGa                                                >  1:16276/1‑149 (MQ=255)
                                                                                                tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa                                                >  2:1132852/1‑149 (MQ=255)
                                                                                                      tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                          >  1:800819/1‑149 (MQ=255)
                                                                                                       ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCg                                         <  1:1132852/149‑1 (MQ=255)
                                                                                                         ctTTTAGAGCCCCCCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCCACAGCAATATGCGTTTAGTAGTAGTAGACCGATAAAGCGaa                                       <  1:562653/149‑1 (MQ=255)
                                                                                                            ttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGa                                    >  2:86835/1‑149 (MQ=255)
                                                                                                             tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                   <  2:568035/149‑1 (MQ=255)
                                                                                                             cagagCGCCTCGCATCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCGCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                   <  2:112093/148‑1 (MQ=255)
                                                                                                              agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                  <  1:104916/149‑1 (MQ=255)
                                                                                                                          ttCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTgaa                      >  2:776886/1‑149 (MQ=255)
                                                                                                                           tCGGGCATAAAAAATCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATACAGCGAACGATGTGAGTCATTAAAt                     <  2:24464/149‑1 (MQ=255)
                                                                                                                           tCGGGCATAAAAAAACCCGCGCAATGCCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                     <  1:926289/149‑1 (MQ=255)
                                                                                                                            cGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc                    <  1:1063649/149‑1 (MQ=255)
                                                                                                                              ggCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                  <  1:366929/149‑1 (MQ=255)
                                                                                                                               gCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGc                 <  1:409511/149‑1 (MQ=255)
                                                                                                                                cATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCt                >  1:478556/1‑149 (MQ=255)
                                                                                                                                  acaaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAAAGCTGGATTCTGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc              >  1:1023388/3‑149 (MQ=255)
                                                                                                                                   aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAGGATGAGAGTCATTAAATCAGCACCa             >  2:501835/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            >  1:363248/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCAATGGCGCTGGTTTGTTGTTTGTCTTCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            <  2:232331/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTAAATCGTCTACCGCGCAGAGTAGAGTCGTTGAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            <  1:1133543/149‑1 (MQ=255)
                                                                                                                                       aagcccGCGCAATGGCCCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg         <  2:481420/146‑1 (MQ=255)
                                                                                                                                           ccGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg     >  2:844108/1‑149 (MQ=255)
                                                                                                                                              cgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa  <  1:86835/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGA  >  NC_000913/83442‑83732

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