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,2720AC97.5% 18803.0 / 279.2 5189intergenic (+7/‑79)araC/yabIDNA‑binding transcriptional dual regulator AraC/DedA family protein YabI
Reads supporting (aligned to +/- strand):  ref base A (60/69);  new base C (2557/2500);  total (2619/2570)
Fisher's exact test for biased strand distribution p-value = 3.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 5194 total aligned reads displayed.

GAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACA  >  NC_000913/71174‑71370
                                                                                                  |                                                                                                  
   tATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGa                                                                                                 >  1:2000735/1‑99 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAAt                                                                                               <  2:85551/99‑1 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  1:2655166‑M1/1‑97 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  1:2299984‑M1/1‑97 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  1:2092320‑M1/1‑97 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  2:1184976‑M1/1‑97 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  2:169490‑M1/1‑97 (MQ=255)
     tttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgc                                                                                               >  1:1195694‑M1/1‑97 (MQ=255)
      ttAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgct                                                                                              >  1:973644‑M1/1‑96 (MQ=255)
          aaaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaa                                                                                          <  2:2211460‑M1/99‑8 (MQ=255)
           aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaat                                                                                         <  2:2307146‑M1/99‑9 (MQ=255)
           aaaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaat                                                                                         >  1:715432‑M1/1‑91 (MQ=255)
            aaTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaatt                                                                                        <  1:2707993‑M1/99‑10 (MQ=255)
             aTGCACCGGGGCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattc                                                                                       <  2:1851791‑M1/99‑11 (MQ=255)
                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                                <  2:2357422‑M1/99‑18 (MQ=255)
                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                                <  1:2811439‑M1/99‑18 (MQ=255)
                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                                <  2:1804652‑M1/99‑18 (MQ=255)
                    ggggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatca                                                                                <  2:246009‑M1/99‑18 (MQ=255)
                     gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcat                                                                               <  1:2040641‑M1/99‑19 (MQ=255)
                     gggCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcat                                                                               >  2:598790‑M1/1‑81 (MQ=255)
                       gCCAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgt                                                                             >  2:1029586‑M1/1‑79 (MQ=255)
                        ccAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtc                                                                            <  1:1764066‑M1/99‑22 (MQ=255)
                        ccAGCCCGAGCGAGTTCCGTGCCGGTGGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtc                                                                            >  1:901224‑M1/1‑78 (MQ=255)
                         cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtca                                                                           <  2:1645546‑M1/99‑23 (MQ=255)
                         cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtca                                                                           <  1:2559467‑M1/99‑23 (MQ=255)
                         cAGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtca                                                                           <  1:123790‑M1/99‑23 (MQ=255)
                          aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcat                                                                          >  1:185162‑M1/1‑76 (MQ=255)
                          aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcat                                                                          >  2:1996919‑M1/1‑76 (MQ=255)
                          aGCCCGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcat                                                                          <  2:1354914‑M1/99‑24 (MQ=255)
                             ccGAGCGAGTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagc                                                                       <  2:2404633‑M1/99‑27 (MQ=255)
                                 gcgagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtt                                                                   <  2:2679960‑M1/99‑31 (MQ=255)
                                   gagTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttc                                                                 >  1:2141427‑M1/1‑67 (MQ=255)
                                    agTTCCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcc                                                                <  1:2515434‑M1/99‑34 (MQ=255)
                                     gttccgtgcCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcct                                                               <  1:1053729‑M1/95‑35 (MQ=255)
                                      ttccgtgCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctg                                                              >  1:1284602‑M1/4‑64 (MQ=255)
                                      ttccgtgCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctg                                                              >  1:573941‑M1/4‑64 (MQ=255)
                                       tccgtGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgt                                                             <  2:2251695‑M1/97‑37 (MQ=255)
                                           tGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtga                                                         >  1:1811910‑M1/1‑59 (MQ=255)
                                           tGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtga                                                         >  2:195442‑M1/1‑59 (MQ=255)
                                           tGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtga                                                         >  1:855774‑M1/1‑59 (MQ=255)
                                             ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                       <  2:1595261‑M1/99‑43 (MQ=255)
                                             ccGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaa                                                       <  1:527210‑M1/99‑43 (MQ=255)
                                              cGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaat                                                      <  2:388324‑M1/99‑44 (MQ=255)
                                               ggTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaatt                                                     <  2:2575408‑M1/99‑45 (MQ=255)
                                                gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattg                                                    >  1:1379554‑M1/1‑54 (MQ=255)
                                                gTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattg                                                    >  2:703408‑M1/1‑54 (MQ=255)
                                                 ttGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgt                                                   >  2:649314‑M1/1‑53 (MQ=255)
                                                   gtgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgtta                                                 >  1:1430352‑M1/1‑51 (MQ=255)
                                                    tgAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttat                                                <  1:477923‑M1/99‑50 (MQ=255)
                                                     gaagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatc                                               <  2:2776103‑M1/99‑51 (MQ=255)
                                                      aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                              <  1:326365‑M1/99‑52 (MQ=255)
                                                      aagaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatcc                                              <  1:2352373‑M1/99‑52 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             >  1:1543180‑M1/1‑47 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             >  2:1500616‑M1/1‑47 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             >  2:284342‑M1/1‑47 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             >  1:375186‑M1/1‑47 (MQ=255)
                                                       agaaAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccg                                             >  1:753491‑M1/1‑47 (MQ=255)
                                                         aaaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgct                                           >  1:2455621‑M1/1‑45 (MQ=255)
                                                         aaaaaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgct                                           >  1:2398202‑M1/1‑45 (MQ=255)
                                                            aaGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcac                                        >  1:2246738‑M1/1‑42 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      <  2:2729272‑M1/99‑60 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      <  2:2095246‑M1/99‑60 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      <  1:1905844‑M1/99‑60 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      <  2:1134571‑M1/99‑60 (MQ=255)
                                                              gTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaa                                      <  1:801318‑M1/99‑60 (MQ=255)
                                                                 aaTGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCAtt                                   >  1:629673/1‑99 (MQ=255)
                                                                   tgatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattcca                                 <  1:1337074‑M1/99‑65 (MQ=255)
                                                                    gatgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccac                                >  1:674487‑M1/1‑34 (MQ=255)
                                                                     atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                               >  2:901971‑M1/1‑33 (MQ=255)
                                                                     atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                               >  1:426956‑M1/1‑33 (MQ=255)
                                                                     atgTAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccaca                               <  2:2054024‑M1/99‑67 (MQ=255)
                                                                        tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaa                            >  2:751171‑M1/1‑30 (MQ=255)
                                                                        tAGCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaa                            >  2:1085559‑M1/1‑30 (MQ=255)
                                                                          gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                          <  2:1546418‑M1/99‑72 (MQ=255)
                                                                          gCCGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaaca                          >  1:1609527‑M1/1‑28 (MQ=255)
                                                                            cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                        <  2:1901597‑M1/99‑74 (MQ=255)
                                                                            cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                        <  1:1956539‑M1/99‑74 (MQ=255)
                                                                            cGTCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacata                        <  1:1099710‑M1/99‑74 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                      >  1:1148355‑M1/1‑24 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                      >  2:2623567‑M1/1‑24 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                      >  1:2760193‑M1/1‑24 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                      >  1:1486630‑M1/1‑24 (MQ=255)
                                                                              tCAAGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacg                      >  1:2197340‑M1/1‑24 (MQ=255)
                                                                               cAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTa                     <  2:1449884/99‑1 (MQ=255)
                                                                                aaGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgag                    <  2:2150001‑M1/99‑78 (MQ=255)
                                                                                aaGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgag                    >  2:54928‑M1/1‑22 (MQ=255)
                                                                                 aGTTGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagc                   <  1:1709073‑M1/99‑79 (MQ=255)
                                                                                   ttGTCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccg                 <  1:1652642‑M1/99‑81 (MQ=255)
                                                                                      tCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaa              <  2:342802‑M1/99‑84 (MQ=255)
                                                                                      tCATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaa              >  2:80465‑M1/1‑16 (MQ=255)
                                                                                       cATAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaag             >  1:1000958‑M1/1‑15 (MQ=255)
                                                                                        aTAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagc            >  2:1292596‑M1/1‑14 (MQ=255)
                                                                                         tAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagca           <  2:50339‑M1/99‑87 (MQ=255)
                                                                                         tAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagca           >  1:2605346‑M1/1‑13 (MQ=255)
                                                                                         tAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagca           >  2:1239212‑M1/1‑13 (MQ=255)
                                                                                         tAATTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagca           >  2:1748168‑M1/1‑13 (MQ=255)
                                                                                         tAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCAtt           >  2:443673/1‑99 (MQ=255)
                                                                                          aaTTGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcat          >  2:952156‑M1/1‑12 (MQ=255)
                                                                                             tGGTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaa       <  2:1403306‑M1/99‑91 (MQ=255)
                                                                                              ggTACCGAgctcgaattcgtaatcatgtcatagctgtttcctgtgtgaaattgttatccgctcacaattccacacaacatacgagccggaagcataaag      <  2:2522077‑M1/99‑92 (MQ=255)
                                                                                                  |                                                                                                  
GAGTATTTAAAAAATGCACCGGGGCCAGCCCGAGCGAGTTTCGTGCCGGTTGTGAAGAAAAAGTGAATGATGTAGCCGTCAAGTTGTCATAATTGGTAACGAATCAGACAATTGACGGCTTGACGGAGTAGCATAGGGTTTGCAGAATCCCTGCTTCGTCCATTTGACAGGCACATTATGCAAGCATTGCTGGAACA  >  NC_000913/71174‑71370

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: