Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 71,272 A→C 100% 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,2720AC99.0% 30639.9 / 138.7 9018intergenic (+7/‑79)araC/yabIDNA‑binding transcriptional dual regulator AraC/DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (59/32);  new base C (4510/4414);  total (4571/4447)
Fisher's exact test for biased strand distribution p-value = 8.10e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 9019 total aligned reads displayed.

GTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAAC  >  NC_000913/71176‑71369
                                                                                                |                                                                                                 
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGCTGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                <  1:961098/99‑1 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                >  1:360105/1‑99 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                >  2:1490836/1‑99 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                >  2:368606/1‑99 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                <  1:1655934/99‑1 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                <  1:1543242/99‑1 (MQ=255)
 tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                >  2:847454/1‑99 (MQ=255)
  aTTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAg                                                                                               >  2:1772221‑M1/1‑98 (MQ=255)
   tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                              >  1:1033181‑M1/1‑97 (MQ=255)
   tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                              >  1:769033‑M1/1‑97 (MQ=255)
   tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                              >  2:483887‑M1/1‑97 (MQ=255)
      aaaaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcg                                                                                           >  1:214185‑M1/1‑94 (MQ=255)
       aaaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcga                                                                                          >  1:620704‑M1/1‑93 (MQ=255)
        aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaa                                                                                         <  2:932836‑M1/99‑8 (MQ=255)
         aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaat                                                                                        <  2:1992805‑M1/99‑9 (MQ=255)
          aaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaatt                                                                                       >  2:219805‑M1/1‑90 (MQ=255)
          aaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaatt                                                                                       >  1:1576392‑M1/1‑90 (MQ=255)
           aTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattc                                                                                      <  1:11462‑M1/99‑11 (MQ=255)
            tGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcg                                                                                     <  1:193694‑M1/99‑12 (MQ=255)
             gCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgt                                                                                    >  1:698516‑M1/1‑87 (MQ=255)
             gCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgt                                                                                    <  2:1068449‑M1/99‑13 (MQ=255)
               aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaa                                                                                  <  1:1110660‑M1/99‑15 (MQ=255)
               aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaa                                                                                  <  1:320701‑M1/99‑15 (MQ=255)
               aCCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaa                                                                                  <  2:1957449‑M1/99‑15 (MQ=255)
                ccGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaat                                                                                 <  1:581654‑M1/99‑16 (MQ=255)
                ccGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaat                                                                                 <  1:129564‑M1/99‑16 (MQ=255)
                 cGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatc                                                                                >  2:1887974‑M1/1‑83 (MQ=255)
                 cGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatc                                                                                >  2:1033264‑M1/1‑83 (MQ=255)
                  ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               <  1:1435646‑M1/99‑18 (MQ=255)
                  ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               <  1:1181919‑M1/99‑18 (MQ=255)
                  ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               >  1:2018368‑M1/1‑82 (MQ=255)
                  ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               <  1:1904101‑M1/99‑18 (MQ=255)
                  ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               <  2:2026791‑M1/99‑18 (MQ=255)
                  cgggcCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                               <  1:851210‑M1/98‑18 (MQ=255)
                     gCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgt                                                                            <  1:1399295‑M1/99‑21 (MQ=255)
                       cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtca                                                                          <  1:1722565‑M1/99‑23 (MQ=255)
                        aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcat                                                                         >  2:565199‑M1/1‑76 (MQ=255)
                         gCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcata                                                                        >  1:1224672‑M1/1‑75 (MQ=255)
                          cccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatag                                                                       <  1:388461‑M1/99‑26 (MQ=255)
                            cgagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagct                                                                     >  2:636926‑M1/1‑72 (MQ=255)
                             gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctg                                                                    <  2:1196226‑M1/99‑29 (MQ=255)
                             gagcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctg                                                                    >  2:250635‑M1/1‑71 (MQ=255)
                              agcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgt                                                                   >  2:137403‑M1/1‑70 (MQ=255)
                              agcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgt                                                                   <  2:441238‑M1/99‑30 (MQ=255)
                                 gagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttc                                                                <  1:1868621‑M1/99‑33 (MQ=255)
                                       cGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgt                                                          <  2:141774‑M1/99‑39 (MQ=255)
                                       cGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgt                                                          <  1:1829646‑M1/99‑39 (MQ=255)
                                         tGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtga                                                        >  1:659854‑M1/1‑59 (MQ=255)
                                           ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                      <  2:29030‑M1/99‑43 (MQ=255)
                                           ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                      <  2:771405‑M1/99‑43 (MQ=255)
                                           ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                      <  2:1850248‑M1/99‑43 (MQ=255)
                                           ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                      <  2:1294869‑M1/99‑43 (MQ=255)
                                             ggTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaatt                                                    >  2:407167‑M1/1‑55 (MQ=255)
                                              gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattg                                                   >  2:1569707‑M1/1‑54 (MQ=255)
                                              gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattg                                                   >  2:1632808‑M1/1‑54 (MQ=255)
                                                 gtgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgtta                                                <  2:706083‑M1/99‑49 (MQ=255)
                                                  tgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttat                                               >  1:1364115‑M1/1‑50 (MQ=255)
                                                  tgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttat                                               >  1:1000540‑M1/1‑50 (MQ=255)
                                                  tgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttat                                               <  2:1917707‑M1/99‑50 (MQ=255)
                                                   gaagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatc                                              <  1:1975451‑M1/99‑51 (MQ=255)
                                                    aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                             >  1:890234‑M1/1‑48 (MQ=255)
                                                    aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                             >  1:1503348‑M1/1‑48 (MQ=255)
                                                    aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                             <  1:807814‑M1/99‑52 (MQ=255)
                                                    aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                             >  2:1226508‑M1/1‑48 (MQ=255)
                                                    aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                             >  1:43001‑M1/1‑48 (MQ=255)
                                                     agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                            >  1:465450‑M1/1‑47 (MQ=255)
                                                     agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                            >  2:1602452‑M1/1‑47 (MQ=255)
                                                      gaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgc                                           >  1:920207‑M1/1‑46 (MQ=255)
                                                        aaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctc                                         >  2:1701574‑M1/1‑44 (MQ=255)
                                                         aaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctca                                        >  2:152848‑M1/1‑43 (MQ=255)
                                                         aaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctca                                        <  1:249223‑M1/99‑57 (MQ=255)
                                                           aGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcaca                                      >  2:1453736‑M1/1‑41 (MQ=255)
                                                            gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                     >  2:671383‑M1/1‑40 (MQ=255)
                                                             tGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaat                                    >  2:327255‑M1/1‑39 (MQ=255)
                                                               aaTGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCAtt                                  <  2:526480/99‑1 (MQ=255)
                                                                atgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcc                                 >  2:1334175‑M1/1‑36 (MQ=255)
                                                                atgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcc                                 >  2:892750‑M1/1‑36 (MQ=255)
                                                                  gatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccac                               >  1:281184‑M1/1‑34 (MQ=255)
                                                                   atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                              >  1:1331172‑M1/1‑33 (MQ=255)
                                                                    tgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacac                             >  2:1736007‑M1/1‑32 (MQ=255)
                                                                      tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaa                           <  2:110850‑M1/99‑70 (MQ=255)
                                                                        gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                         >  1:993824‑M1/1‑28 (MQ=255)
                                                                        gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                         >  2:1812467‑M1/1‑28 (MQ=255)
                                                                        gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                         >  2:1153534‑M1/1‑28 (MQ=255)
                                                                        gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                         <  1:1472455‑M1/99‑72 (MQ=255)
                                                                        gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                         <  1:501591‑M1/99‑72 (MQ=255)
                                                                          cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                       <  1:543977‑M1/99‑74 (MQ=255)
                                                                          cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                       <  1:175897‑M1/99‑74 (MQ=255)
                                                                           gTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatac                      >  2:597550‑M1/1‑25 (MQ=255)
                                                                            tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                     >  1:1792482‑M1/1‑24 (MQ=255)
                                                                            tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaatcccacacaacatacg                     >  2:1257387‑M1/1‑24 (MQ=255)
                                                                              aaGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgag                   <  1:1621614‑M1/99‑78 (MQ=255)
                                                                                gTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagcc                 <  2:1665714‑M1/99‑80 (MQ=255)
                                                                                gTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagcc                 >  2:810662‑M1/1‑20 (MQ=255)
                                                                                 ttGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccg                >  2:740090‑M1/1‑19 (MQ=255)
                                                                                       tAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagca          >  1:1067318‑M1/1‑13 (MQ=255)
                                                                                        aaTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcat         >  1:1691501‑M1/1‑12 (MQ=255)
                                                                                          ttGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataa       <  1:1259328‑M1/99‑90 (MQ=255)
                                                                                          ttGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataa       >  2:967650‑M1/1‑10 (MQ=255)
                                                                                            ggTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaag     <  1:734361‑M1/99‑92 (MQ=255)
                                                                                                |                                                                                                 
GTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAAC  >  NC_000913/71176‑71369

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