Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 83,587 A→C 17.9% 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,5870AC17.9% 329.4 / 69.2 150intergenic (‑58/+35)leuA/leuL2‑isopropylmalate synthase/leu operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (55/68);  new base C (13/14);  total (68/82)
Fisher's exact test for biased strand distribution p-value = 8.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 151 total aligned reads displayed.

CGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGA  >  NC_000913/83490‑83686
                                                                                                 |                                                                                                    
cgcgcgATGTGGTATCGAAACTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGGGCCTCGCTTGGGGAATAAAAAAACCCGCGTAAt                                                                                                    >  1:2388152/1‑100 (MQ=255)
 gcgcAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                   <  1:704328/100‑1 (MQ=255)
 gcgcAATGTGGTATCGAAAATAATGACTTGCTGGCGCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCAGGCATAAAAAAACCCGCGCAATg                                                                                                   >  1:595971/1‑100 (MQ=255)
 gcgcAATGTGGTATCGAAAAAAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATg                                                                                                   >  1:1586199/1‑100 (MQ=255)
    cAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGgc                                                                                                 <  1:2503466/99‑1 (MQ=255)
     aaTGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGGCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTgcg                                                                                                <  1:136369/99‑1 (MQ=255)
          ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tt                                                                                         <  2:2203185/100‑1 (MQ=255)
          ggTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tt                                                                                         >  2:1656912/1‑100 (MQ=255)
            tATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCAT‑AAAAAACCCGCGCCATTGCGCGGG‑ttttt                                                                                      <  2:268941/100‑1 (MQ=38)
             aTCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttt                                                                                      <  1:2283524/100‑1 (MQ=255)
              tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttt                                                                                     <  2:1172844/100‑1 (MQ=255)
              tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttt                                                                                     <  1:2583637/100‑1 (MQ=255)
              tCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttt                                                                                     <  1:1042365/100‑1 (MQ=255)
              tCGAAAATAATGACTTGCGGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑tttttt                                                                                     <  1:1877508/100‑1 (MQ=255)
               cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttg                                                                                    >  1:166937/1‑100 (MQ=255)
               cGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttg                                                                                    <  2:1031294/100‑1 (MQ=255)
                tccgatcttGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttttta                                                                                   <  1:2592531/92‑3 (MQ=37)
                gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttttta                                                                                   <  2:2573414/100‑3 (MQ=255)
                gAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑ttttttt                                                                                    <  1:137797/99‑2 (MQ=255)
                 aaaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCCATTGCGCGGG‑tttttttat                                                                                  <  2:1977912/100‑4 (MQ=38)
                 aaaaCAATGACTTGCTCGCTCATGGTTTGGGTCCTTCTCTCTTTTAGAGCGCCTCGCTTCGTGCATAAAATAACCCTCGCCAGTGCGCGGG‑tttcttt                                                                                    <  1:1272810/98‑2 (MQ=255)
                  aaaTAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTtttgttt                                                                                 <  1:1297894/100‑1 (MQ=255)
                      aatGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACt                                                                             >  2:2362094/1‑100 (MQ=255)
                       atGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTATTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGAGCGGG‑TTTTTTGTTTGACTg                                                                            >  2:1527053/1‑100 (MQ=255)
                           cTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTg                                                                        >  2:73129/1‑100 (MQ=255)
                             tGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCt                                                                      <  1:1432029/100‑1 (MQ=255)
                                 ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCt                                                                  >  1:2304205/1‑100 (MQ=255)
                                 ggCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCt                                                                  <  2:402963/100‑1 (MQ=255)
                                     cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAAt                                                              >  2:523990/1‑100 (MQ=255)
                                     cATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAAt                                                              <  2:1905434/100‑1 (MQ=255)
                                      aTGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                             >  1:1514821/1‑100 (MQ=255)
                                      aTGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATg                                                             <  1:2370438/100‑1 (MQ=255)
                                       tGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGc                                                            >  2:2139094/1‑100 (MQ=255)
                                        tgTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCt                                                           <  1:344348/99‑1 (MQ=255)
                                         gTTTGGGTCCTTGTCTCTTTTAGAGCTCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCtt                                                          <  2:1277212/100‑2 (MQ=255)
                                         gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                          <  2:474959/100‑1 (MQ=255)
                                         gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                          <  2:2342795/100‑1 (MQ=255)
                                         gTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTg                                                          <  1:1723082/100‑1 (MQ=255)
                                          tttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTgg                                                         >  2:1426356/1‑100 (MQ=255)
                                          tttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTgg                                                         <  2:336128/100‑1 (MQ=255)
                                          tctGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTgg                                                         <  2:878232/98‑1 (MQ=255)
                                           ttGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGa                                                        <  1:228143/100‑1 (MQ=255)
                                            tGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGAt                                                       <  2:1237945/100‑1 (MQ=255)
                                             gggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGaga                                                      >  1:878229/1‑98 (MQ=255)
                                              ggTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGc                                                     <  2:1127845/100‑1 (MQ=255)
                                               gTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGcc                                                    <  2:2042943/100‑1 (MQ=255)
                                                tCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCg                                                   <  1:118684/100‑1 (MQ=255)
                                                  cTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCt                                                 <  2:1246188/100‑1 (MQ=255)
                                                    tGTCTCTTTTTGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTc                                                <  1:339127/99‑1 (MQ=255)
                                                     gTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                              >  2:2566113/1‑100 (MQ=255)
                                                     gTCTATTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGAGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                              <  1:1527044/100‑1 (MQ=255)
                                                       ctctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCAc                                              <  1:1568317/98‑1 (MQ=255)
                                                        tctTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCg                                           >  2:1642520/1‑100 (MQ=255)
                                                           tttAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCt                                        >  1:2041128/1‑100 (MQ=255)
                                                             tAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAc                                      >  1:2326598/1‑100 (MQ=255)
                                                              agagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGCCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTAcc                                     >  1:2110830/1‑100 (MQ=255)
                                                               gagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcg                                    <  2:842830/100‑1 (MQ=255)
                                                               gagCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcg                                    <  2:1764267/100‑1 (MQ=255)
                                                                 gcgcCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACcgcg                                  >  1:1106666/1‑100 (MQ=255)
                                                                    ccTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaa                               >  2:845702/1‑100 (MQ=255)
                                                                     cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaaa                              <  1:2311342/100‑1 (MQ=255)
                                                                     cTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCaaa                              >  1:143154/1‑100 (MQ=255)
                                                                        gCTTCGGGCATAAAAAATCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaaga                            <  1:2535358/99‑1 (MQ=255)
                                                                        gCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaa                           <  2:702071/100‑1 (MQ=255)
                                                                         cTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAaagaag                          <  1:1743730/100‑1 (MQ=255)
                                                                          ttCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAagaaga                         <  1:98842/100‑1 (MQ=255)
                                                                           tCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGAt                        <  1:819792/100‑1 (MQ=255)
                                                                           tCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGAt                        >  1:760268/1‑100 (MQ=255)
                                                                            cGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGAt                        <  2:799473/99‑1 (MQ=255)
                                                                             gtcacacaaaaaaCCCGCGCCATTGCGCGGG‑CTTTTTGTTTGACTGCGTGCTGGCTTAATGCGGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATg                       >  1:2339140/8‑99 (MQ=37)
                                                                             gggCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGc                      >  1:415013/1‑100 (MQ=255)
                                                                               gCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                    <  2:426424/100‑1 (MQ=255)
                                                                                aaacaaaaaaCCCGCGCCATTGCGCGGGTTTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGt                    >  1:2213181/5‑100 (MQ=37)
                                                                                 tctaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                  <  2:71592/97‑1 (MQ=255)
                                                                                 atagtaaagCCCGCGCAATGGCGCGGG‑TTTTGTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                  <  1:1811174/95‑1 (MQ=255)
                                                                                 aacaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                  >  2:2173264/4‑100 (MQ=38)
                                                                                 aTAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGttt                  <  1:1165689/100‑1 (MQ=255)
                                                                                  acaaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGGCGCTCACTCGTCTACCGCGCAAAGAAGATggagtta                 >  1:484487/3‑94 (MQ=37)
                                                                                   caaaaaaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                >  2:2216835/2‑100 (MQ=255)
                                                                                   aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                <  2:579423/100‑1 (MQ=255)
                                                                                   aaaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtag                <  1:568996/100‑1 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               <  1:2139080/100‑1 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               >  2:1619396/1‑100 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               <  2:955/100‑1 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               >  2:2056405/1‑100 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               <  2:107321/100‑1 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagt               <  2:1315591/100‑1 (MQ=255)
                                                                                    aaaaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtaga               >  1:71591/1‑99 (MQ=255)
                                                                                      aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag             >  1:1965754/1‑100 (MQ=255)
                                                                                      aaaaCCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtag             <  1:1656903/100‑1 (MQ=255)
                                                                                        aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTTCTTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta           >  1:2054122/1‑100 (MQ=39)
                                                                                        aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGTCTTACTGCTGTATGCCGCTCACTCGTCTACCGCGCAAAGATGATCCGATtagtagta           >  2:1264569/1‑100 (MQ=37)
                                                                                        aaCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagta           >  2:1038788/1‑100 (MQ=255)
                                                                                         aCCCGCGCCATTGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtag          >  2:1316231/1‑100 (MQ=255)
                                                                                          cccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         <  2:2507152/100‑1 (MQ=255)
                                                                                          cccGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         <  2:836747/100‑1 (MQ=255)
                                                                                            cgcgCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTtagtagtagt         <  1:69514/98‑1 (MQ=255)
                                                                                              cgcAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGAc     <  2:2460943/100‑1 (MQ=255)
                                                                                                cAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCg   >  2:914363/1‑100 (MQ=255)
                                                                                                 aaTGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGa  <  2:1145106/100‑1 (MQ=255)
                                                                                                 |                                                                                                    
CGCGCAATGTGGTATCGAAAATAATGACTTGCTGGCTCATGGTTTGGGTCCTTGTCTCTTTTAGAGCGCCTCGCTTCGGGCATAAAAAAACCCGCGCAATGGCGCGGG‑TTTTTTGTTTGACTGCGTGCTGGCTTAATGCTGGATGCCGCTCACTCGTCTACCGCGCAAAGAAGATGCGTTTAGTAGTAGTAGACCGA  >  NC_000913/83490‑83686

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