Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 366,290 C→T 100% D6N (GAT→AAT)  lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,2900CT99.0% 41691.7 / 124.3 11689D6N (GAT→AAT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (52/68);  new base T (5799/5765);  total (5853/5836)
Fisher's exact test for biased strand distribution p-value = 1.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 11691 total aligned reads displayed.

CCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366195‑366388
                                                                                               |                                                                                                  
  accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                               >  2:2974154‑M2/90‑99 (MQ=255)
   ccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAt                                                                                              <  1:2607006‑M2/11‑1 (MQ=255)
     ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  1:2772582‑M2/13‑1 (MQ=255)
     ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                            <  1:2266632‑M2/13‑1 (MQ=255)
        gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         >  2:2520475‑M2/84‑98 (MQ=255)
        gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         <  1:1482707‑M2/16‑2 (MQ=255)
        gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                         <  2:433177‑M2/16‑2 (MQ=255)
          cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                       <  1:736227‑M2/18‑4 (MQ=255)
          cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                       <  2:722742‑M2/18‑4 (MQ=255)
           agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                      >  1:580600‑M2/81‑95 (MQ=255)
            gtccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                    >  2:1280588‑M2/80‑99 (MQ=255)
            gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                    >  2:10019‑M2/80‑99 (MQ=255)
            gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagcccGAATTCGTAATCATGTCATa                                                                                    >  2:1834958‑M2/80‑99 (MQ=255)
             cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                   >  2:1225807‑M2/79‑99 (MQ=255)
             cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                   <  1:2000465‑M2/21‑1 (MQ=255)
             cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                   <  1:1161722‑M2/21‑1 (MQ=255)
              ccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGc                                                                                  >  1:1716080‑M2/78‑99 (MQ=255)
                 agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                               <  1:1768724‑M2/25‑1 (MQ=255)
                  gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                              <  1:100000‑M2/26‑1 (MQ=255)
                  gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                              <  1:2715351‑M2/26‑1 (MQ=255)
                    gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTc                                                                            <  1:838481‑M2/28‑1 (MQ=255)
                     agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTcc                                                                           <  1:1938861‑M2/29‑1 (MQ=255)
                      gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCt                                                                          >  1:789820‑M2/70‑99 (MQ=255)
                       ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                         <  2:1776935‑M2/31‑1 (MQ=255)
                       ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                         >  2:1443094‑M2/69‑99 (MQ=255)
                       ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                         >  2:2466819‑M2/69‑99 (MQ=255)
                       ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                         >  2:845057‑M2/69‑99 (MQ=255)
                        tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                        <  2:2570150‑M2/32‑1 (MQ=255)
                         gcAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCtgt                                                                        >  1:685114/1‑99 (MQ=255)
                          cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                      >  1:1378729‑M2/66‑99 (MQ=255)
                          cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                      >  2:2302240‑M2/66‑99 (MQ=255)
                          cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                      <  2:1336369‑M2/34‑1 (MQ=255)
                            tgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGa                                                                    <  2:1552445‑M2/36‑1 (MQ=255)
                            tgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGa                                                                    >  1:1831264‑M2/64‑99 (MQ=255)
                              ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                  <  1:892409‑M2/38‑1 (MQ=255)
                                  ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                              <  1:1327157‑M2/42‑1 (MQ=255)
                                    gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                            <  2:3026418‑M2/44‑1 (MQ=255)
                                     tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                           >  1:1888829‑M2/55‑99 (MQ=255)
                                     tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                           >  2:1893419‑M2/55‑99 (MQ=255)
                                     tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                           <  2:325398‑M2/45‑1 (MQ=255)
                                        agaaaaagtgaatgatgtagccgtcatgttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                        <  2:942297‑M2/48‑1 (MQ=255)
                                        agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                        >  2:898316‑M2/52‑99 (MQ=255)
                                        agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                        >  2:2848439‑M2/52‑99 (MQ=255)
                                         gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                       >  1:1220128‑M2/51‑99 (MQ=255)
                                          aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                      <  2:1724211‑M2/50‑1 (MQ=255)
                                           aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                     <  1:428721‑M2/51‑1 (MQ=255)
                                            aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                    >  1:528025‑M2/48‑99 (MQ=255)
                                             aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                   >  2:2683006‑M2/47‑99 (MQ=255)
                                             aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                   >  1:2044738‑M2/47‑99 (MQ=255)
                                               gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                 <  2:2407188‑M2/55‑1 (MQ=255)
                                               gtgaatgatgtagccgacaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                 <  1:2821243‑M2/55‑1 (MQ=255)
                                                  aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                              >  2:782086‑M2/42‑99 (MQ=255)
                                                  aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAACCGCTCACAATTc                                              <  2:2176213‑M2/58‑1 (MQ=255)
                                                   atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                             >  1:1045114‑M2/41‑99 (MQ=255)
                                                    tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                            >  1:166446‑M2/40‑99 (MQ=255)
                                                    tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCAAAATTCca                                            >  2:2901980‑M2/40‑99 (MQ=255)
                                                     gatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcac                                           <  2:2008186‑M2/61‑1 (MQ=255)
                                                       tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                         >  1:2370473‑M2/37‑99 (MQ=255)
                                                        gtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacaca                                        <  1:1425030‑M2/64‑1 (MQ=255)
                                                         tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                       >  2:1105889‑M2/35‑99 (MQ=255)
                                                         tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                       >  1:2495324‑M2/35‑99 (MQ=255)
                                                         tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                       <  1:2981947‑M2/65‑1 (MQ=255)
                                                         tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                       <  1:2163730‑M2/65‑1 (MQ=255)
                                                          agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                      <  1:1544844‑M2/66‑1 (MQ=255)
                                                          agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                      >  2:211539‑M2/34‑99 (MQ=255)
                                                          agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                      <  2:2797777‑M2/66‑1 (MQ=255)
                                                             cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                   <  1:3024921‑M2/69‑1 (MQ=255)
                                                             cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                   <  1:1273729‑M2/69‑1 (MQ=255)
                                                               tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:479104‑M2/29‑99 (MQ=255)
                                                               tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:2934885‑M2/29‑99 (MQ=255)
                                                               tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:2877191‑M2/29‑99 (MQ=255)
                                                               tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  1:2663165‑M2/29‑99 (MQ=255)
                                                               tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                 >  2:2625832‑M2/29‑99 (MQ=255)
                                                                caagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGa                                <  2:498074‑M2/72‑1 (MQ=255)
                                                                 aagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                               >  2:1948251‑M2/27‑99 (MQ=255)
                                                                  agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                              <  1:2319786‑M2/74‑1 (MQ=255)
                                                                   gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                             >  2:160597‑M2/25‑99 (MQ=255)
                                                                      gtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGa                          <  2:554948‑M2/78‑1 (MQ=255)
                                                                      gtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGa                          >  1:1099590‑M2/22‑99 (MQ=255)
                                                                        cataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                        <  2:1493973‑M2/80‑1 (MQ=255)
                                                                          taattggtaccgagctcGAATTCGTAGTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTAACAATTCCACACAACATACGAGCCGGAAgct                      >  2:1051603‑M2/18‑98 (MQ=255)
                                                                          taattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                      <  2:2356896‑M2/82‑1 (MQ=255)
                                                                           aattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                     >  1:2428977‑M2/17‑99 (MQ=255)
                                                                           aattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                     >  2:1661655‑M2/17‑99 (MQ=255)
                                                                            attggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    <  2:1163631‑M2/84‑1 (MQ=255)
                                                                            attggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                    <  1:2217711‑M2/84‑1 (MQ=255)
                                                                              tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                  <  1:628587‑M2/86‑1 (MQ=255)
                                                                               ggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                 <  2:612344‑M2/87‑1 (MQ=255)
                                                                                 taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg               >  2:2742740‑M2/11‑99 (MQ=255)
                                                                                 gacgGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                <  2:1385895/99‑1 (MQ=255)
                                                                                   ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa             >  2:382214‑M2/9‑99 (MQ=255)
                                                                                    cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa            <  1:210873‑M2/92‑1 (MQ=255)
                                                                                     gagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa           <  1:376470‑M2/93‑1 (MQ=255)
                                                                                       gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc         >  2:2064535‑M2/5‑99 (MQ=255)
                                                                                          cGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg      <  1:1611742‑M2/98‑1 (MQ=255)
                                                                                           gAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgg     >  2:2236864/1‑99 (MQ=255)
                                                                                            aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg    >  2:671137/1‑99 (MQ=255)
                                                                                            aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg    >  1:948799/1‑99 (MQ=255)
                                                                                            aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg    >  1:2550524/1‑99 (MQ=255)
                                                                                            aaTCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgg     >  1:325823/1‑99 (MQ=255)
                                                                                               |                                                                                                  
CCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGT  >  NC_000913/366195‑366388

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