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

GCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGA  >  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)
                                                                      aatGACTGGCGGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCCCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaag‑a                                                                          <  1:493089/149‑1 (MQ=255)
                                                                       ctcaCTTCCGGGCTCCTGGTTTGCTGCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAATCCCGCGCCATGGTGCGGGTTTTGTGTTTGACTGCGTGCTGGCTGAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaag‑aa                                                                         <  2:214323/146‑1 (MQ=255)
                                                                           cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATg                                                                     <  1:626075/149‑1 (MQ=255)
                                                                              gCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGt                                                                  <  1:280579/149‑1 (MQ=255)
                                                                                tGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGttt                                                                >  1:533477/1‑149 (MQ=255)
                                                                                 ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTta                                                               <  1:518278/149‑1 (MQ=255)
                                                                                 ggCTCATGGTTTGGGTCCTTGTCCCGTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTta                                                               <  1:521941/149‑1 (MQ=255)
                                                                                  gCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtag                                                              >  2:409511/1‑149 (MQ=255)
                                                                                     cATGGTTGGCGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtagtag                                                           <  1:156155/149‑1 (MQ=255)
                                                                                     cATGGTCTGGGTCCTTGTCTCTTTTAGAGCGCTCCGCTTCGGGCATAAAAAACCACGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGATGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtagtag                                                           <  2:637625/149‑1 (MQ=255)
                                                                                       tgttttGGGTCCTTTTCTCTTTTAGGGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtagtagta                                                         <  2:189258/146‑1 (MQ=255)
                                                                                       tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtagtagta                                                         <  2:1086706/149‑1 (MQ=255)
                                                                                        ggTTTGGGTCCCTGTCTCTTTTAGGGCGCCTCGCTTCGGGCATAAAACAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTtagtagtag                                                        <  2:846027/149‑1 (MQ=255)
                                                                                          tttGGGTCTTTGCCTCTTTTAGAGCGCCTCTCTTTCGCCATCAACACACCCGCGCAATGGCGCGGGTTTTTTTTTTGACTGCGTGTTGGCTTAATGCTGGATGCCGCTCACGCGTCTTCCGCGCAAAG‑AAGTTGCGTTtagtagtagta                                                      <  2:544026/149‑1 (MQ=255)
                                                                                            tGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGa                                                    <  2:453259/149‑1 (MQ=255)
                                                                                            tGGCTCCTTTCCTCTTTTAGGCCTCCTCGCTTCCGGCATAAAAAACCCCGAGCAATGGTGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTTATGCTGGATGCCGCTCACTCGTCTAACGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGa                                                    <  2:1007632/149‑1 (MQ=255)
                                                                                                tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTGGACCGa                                                >  1:16276/1‑149 (MQ=255)
                                                                                                tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGa                                                >  2:1132852/1‑149 (MQ=255)
                                                                                                      tctcTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGc                                          >  1:800819/1‑149 (MQ=255)
                                                                                                       ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCg                                         <  1:1132852/149‑1 (MQ=255)
                                                                                                         ctTTTAGAGCCCCCCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGC‑CACAGCAATATGCGTTTAGTAGTAGTAGACCGATAAAGCGaa                                       <  1:562653/149‑1 (MQ=255)
                                                                                                            ttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGa                                    >  2:86835/1‑149 (MQ=255)
                                                                                                             tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                   <  2:568035/149‑1 (MQ=255)
                                                                                                             cagagCGCCTCGCATCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCGCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAt                                   <  2:112093/148‑1 (MQ=255)
                                                                                                              agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGAtg                                  <  1:104916/149‑1 (MQ=255)
                                                                                                                          ttCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTgaa                      >  2:776886/1‑149 (MQ=255)
                                                                                                                           tCGGGCATAAAAAATCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATACAGCGAACGATGTGAGTCATTAAAt                     <  2:24464/149‑1 (MQ=255)
                                                                                                                           tCGGGCATAAAAAAACCCGCGCAATGCCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAAt                     <  1:926289/149‑1 (MQ=255)
                                                                                                                            cGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATc                    <  1:1063649/149‑1 (MQ=255)
                                                                                                                              ggCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAg                  <  1:366929/149‑1 (MQ=255)
                                                                                                                               gCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGc                 <  1:409511/149‑1 (MQ=255)
                                                                                                                                cATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCt                >  1:478556/1‑149 (MQ=255)
                                                                                                                                  acaaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAAAGCTGGATTCTGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTcc              >  1:1023388/3‑149 (MQ=255)
                                                                                                                                   aaaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAAGGATGAGAGTCATTAAATCAGCACCa             >  2:501835/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCCATTGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            >  1:363248/1‑149 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCAATGGCGCTGGTTTGTTGTTTGTCTTCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            <  2:232331/149‑1 (MQ=255)
                                                                                                                                    aaaaaaCCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTAAATCGTCTACCGCGCAGAG‑TAGAGTCGTTGAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAg            <  1:1133543/149‑1 (MQ=255)
                                                                                                                                       aagcccGCGCAATGGCCCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATg         <  2:481420/146‑1 (MQ=255)
                                                                                                                                           ccGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATg     >  2:844108/1‑149 (MQ=255)
                                                                                                                                              cgcAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGa  <  1:86835/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GCAGTTTTTCTTTCACACTCAAGCTTGCCTGTAACGCCTGTTCACCGTCGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGGTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAG‑AAGATGCGTTTAGTAGTAGTAGACCGATAAAGCGAACGATGTGAGTCATTAAATCAGCTCCAGATGAATGCGA  >  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: ‑