Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,272 A→C 82.7% intergenic (+7/‑79) araC → / → yabI DNA‑binding transcriptional dual regulator AraC/DedA family protein YabI

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091371,2720AC82.7% 2162.4 / inf 648intergenic (+7/‑79)araC/yabIDNA‑binding transcriptional dual regulator AraC/DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (58/54);  new base C (276/260);  total (334/314)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

Only 100 of 654 total aligned reads displayed.

GAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACA  >  NC_000913/71174‑71370
                                                                                                  |                                                                                                  
   tatataAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                 >  2:301619/1‑99 (MQ=255)
      ttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgct                                                                                              <  1:461799‑M1/99‑4 (MQ=255)
       tAAAAAATGCACCCGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctc                                                                                             <  2:1581564‑M1/99‑5 (MQ=255)
         aaaaaTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGa                                                                                           <  2:224968/99‑1 (MQ=255)
          aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaa                                                                                          <  2:701071‑M1/99‑8 (MQ=255)
          aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaa                                                                                          >  2:2447994‑M1/1‑92 (MQ=255)
           aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaat                                                                                         <  1:1769780‑M1/99‑9 (MQ=255)
            aaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaatt                                                                                        <  2:2202767‑M1/99‑10 (MQ=255)
             aTGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAAt                                                                                       >  2:500474/1‑99 (MQ=255)
              tGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcg                                                                                      >  1:1392220‑M1/1‑88 (MQ=255)
              tGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcg                                                                                      >  2:1245030‑M1/1‑88 (MQ=255)
                cACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgta                                                                                    <  2:930112‑M1/99‑14 (MQ=255)
                 aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaa                                                                                   <  1:1147524‑M1/99‑15 (MQ=255)
                 aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaa                                                                                   <  1:1632501‑M1/99‑15 (MQ=255)
                  ccGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaat                                                                                  >  1:1909759‑M1/1‑84 (MQ=255)
                     gggCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCt                                                                               <  1:21698/99‑1 (MQ=255)
                      ggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatg                                                                              <  1:951955‑M1/99‑20 (MQ=255)
                      ggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatg                                                                              <  1:334025‑M1/99‑20 (MQ=255)
                         cAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGAc                                                                           >  2:1818257/1‑99 (MQ=255)
                         cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtca                                                                           <  2:2103293‑M1/99‑23 (MQ=255)
                          aGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACg                                                                          >  1:875590/1‑99 (MQ=255)
                          aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcat                                                                          <  1:269802‑M1/99‑24 (MQ=255)
                           gCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcata                                                                         >  2:590054‑M1/1‑75 (MQ=255)
                           gCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcata                                                                         >  1:1330528‑M1/1‑75 (MQ=255)
                             ccGGGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAg                                                                       <  2:438958/99‑1 (MQ=255)
                             ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagc                                                                       <  2:853670‑M1/99‑27 (MQ=255)
                              cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagct                                                                      >  2:815159‑M1/1‑72 (MQ=255)
                               gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctg                                                                     >  1:1737488‑M1/1‑71 (MQ=255)
                               gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctg                                                                     <  2:2307026‑M1/99‑29 (MQ=255)
                                agcgagTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAg                                                                    >  2:734815/1‑99 (MQ=255)
                                agcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgt                                                                    >  2:1951227‑M1/1‑70 (MQ=255)
                                agcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgt                                                                    >  2:2533529‑M1/1‑70 (MQ=255)
                                 gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtt                                                                   <  1:1582291‑M1/99‑31 (MQ=255)
                                 gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtt                                                                   <  1:2493605‑M1/99‑31 (MQ=255)
                                 gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtt                                                                   <  1:228295‑M1/99‑31 (MQ=255)
                                 gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtt                                                                   <  1:1001581‑M1/99‑31 (MQ=255)
                                  cgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgttt                                                                  <  1:1104645‑M1/99‑32 (MQ=255)
                                      tttCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAgg                                                              >  2:1707770/1‑99 (MQ=255)
                                      ttccgtgCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctg                                                              >  1:1069624‑M1/4‑64 (MQ=255)
                                      ttccgtgCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctg                                                              >  1:689087‑M1/4‑64 (MQ=255)
                                       tccgtGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgt                                                             >  1:627753‑M1/3‑63 (MQ=255)
                                        ccgTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtg                                                            <  1:1945624‑M1/98‑38 (MQ=255)
                                         cGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgt                                                           <  1:1834815‑M1/99‑39 (MQ=255)
                                          gTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtg                                                          >  1:800782‑M1/1‑60 (MQ=255)
                                          gTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtg                                                          >  2:1407280‑M1/1‑60 (MQ=255)
                                           tGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTg                                                         >  1:1520823/1‑99 (MQ=255)
                                                gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTTCCGAgctcgaattcgtgatcatgtcatagctgtttcctgtgtgaaattg                                                    <  2:1930206‑M1/99‑46 (MQ=255)
                                                gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattg                                                    >  2:2357901‑M1/1‑54 (MQ=255)
                                                 ttGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgt                                                   >  2:1028686‑M1/1‑53 (MQ=255)
                                                 ttGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgt                                                   >  1:1272255‑M1/1‑53 (MQ=255)
                                                   gtgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgtta                                                 <  1:1869261‑M1/99‑49 (MQ=255)
                                                   gtgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgtta                                                 <  2:1357306‑M1/99‑49 (MQ=255)
                                                      aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                              <  2:120307‑M1/99‑52 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             <  1:1999790‑M1/99‑53 (MQ=255)
                                                        gaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgc                                            >  2:1500638‑M1/1‑46 (MQ=255)
                                                        gaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgc                                            >  1:240770‑M1/1‑46 (MQ=255)
                                                          aaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctc                                          >  1:503494‑M1/1‑44 (MQ=255)
                                                          aaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctc                                          >  1:1600561‑M1/1‑44 (MQ=255)
                                                          aaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTc                                          <  2:2322666/99‑1 (MQ=255)
                                                           aaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctca                                         <  2:883144‑M1/99‑57 (MQ=255)
                                                           aaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctca                                         >  1:570655‑M1/1‑43 (MQ=255)
                                                            aaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcac                                        >  1:2466936‑M1/1‑42 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      >  1:2323889‑M1/1‑40 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      >  1:777624‑M1/1‑40 (MQ=255)
                                                                gAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaatt                                    >  1:1462163‑M1/1‑38 (MQ=255)
                                                                 aaTGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattc                                   >  1:2505249‑M1/1‑37 (MQ=255)
                                                                  atgatgTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCAttt                                  <  2:1068092/99‑1 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                 <  1:725767‑M1/99‑65 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                 <  1:587267‑M1/99‑65 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                 <  1:2367959‑M1/99‑65 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                 <  1:372780‑M1/99‑65 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTg                                 >  1:2249709/1‑99 (MQ=255)
                                                                    gatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccac                                >  2:1455530‑M1/1‑34 (MQ=255)
                                                                     atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                               <  1:2535826‑M1/99‑67 (MQ=255)
                                                                     atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                               <  2:764070‑M1/99‑67 (MQ=255)
                                                                      tgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacac                              <  2:1545892‑M1/99‑68 (MQ=255)
                                                                        tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaa                            >  2:2474795‑M1/1‑30 (MQ=255)
                                                                        tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaa                            <  2:2052325‑M1/99‑70 (MQ=255)
                                                                        tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaataccacacaa                            <  1:1029825‑M1/99‑70 (MQ=255)
                                                                         aGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaac                           <  1:2049723‑M1/99‑71 (MQ=255)
                                                                         aGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaac                           <  1:1242932‑M1/99‑71 (MQ=255)
                                                                          gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                          <  2:474144‑M1/99‑72 (MQ=255)
                                                                          gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                          <  2:2377367‑M1/99‑72 (MQ=255)
                                                                            cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                        >  2:1162114‑M1/1‑26 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACAtt                      >  2:788437/1‑99 (MQ=255)
                                                                               cAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTa                     <  1:1809175/99‑1 (MQ=255)
                                                                                 aGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagc                   <  1:153436‑M1/99‑79 (MQ=255)
                                                                                 aGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagc                   <  2:1745221‑M1/99‑79 (MQ=255)
                                                                                 aGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagc                   <  1:211102‑M1/99‑79 (MQ=255)
                                                                                  gTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagcc                  >  2:1309766‑M1/1‑20 (MQ=255)
                                                                                   ttGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttagccgctcacaattccacacaacatacgagccg                 <  2:1640040‑M1/99‑81 (MQ=255)
                                                                                      tCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaa              >  2:55791‑M1/1‑16 (MQ=255)
                                                                                      tCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGc              >  2:1891973/1‑99 (MQ=255)
                                                                                       cATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaag             >  2:406425‑M1/1‑15 (MQ=255)
                                                                                        aTAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagc            >  1:1659173‑M1/1‑14 (MQ=255)
                                                                                          aaTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcat          >  2:1129598‑M1/1‑12 (MQ=255)
                                                                                          aaTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcat          <  1:1170840‑M1/99‑88 (MQ=255)
                                                                                           aTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcata         >  2:1678952‑M1/1‑11 (MQ=255)
                                                                                           aTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcata         >  2:2271598‑M1/1‑11 (MQ=255)
                                                                                             tGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaa       <  1:2228984‑M1/99‑91 (MQ=255)
                                                                                                  |                                                                                                  
GAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACA  >  NC_000913/71174‑71370

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