Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 366,300 Δ1 bp 100% coding (6/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,3000G.95.1% 13675.5 / inf 3300coding (6/3075 nt)lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (85/75);  new base . (1593/1546);  total (1678/1622)
Fisher's exact test for biased strand distribution p-value = 5.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 3470 total aligned reads displayed.

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398
                                                                                                  |                                                                                                  
ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                                    >  1:1964823‑M2/85‑99 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTTGTAATCAtgt                                                                                                   <  2:1665187‑M2/16‑2 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   <  1:1683308‑M2/16‑2 (MQ=255)
    agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                                >  2:91953‑M2/81‑95 (MQ=255)
    agcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAtcatgtcat                                                                                                >  1:609914‑M2/81‑95 (MQ=255)
         gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                                          >  2:1188419‑M2/76‑99 (MQ=255)
          agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         >  2:2153541‑M2/75‑99 (MQ=255)
          agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         >  2:1860580‑M2/75‑99 (MQ=255)
           gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                                        <  1:1000384‑M2/26‑1 (MQ=255)
            tgtgCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGtt                                                                                        >  2:1260564/1‑99 (MQ=255)
             gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTc                                                                                      >  2:530752‑M2/72‑99 (MQ=255)
                ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                                   >  1:132715‑M2/69‑99 (MQ=255)
                ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                                   <  1:807123‑M2/31‑1 (MQ=255)
                 tgcAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCtg                                                                                   <  1:1471782/99‑1 (MQ=255)
                       ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                            >  2:1221572‑M2/62‑99 (MQ=255)
                        cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                           <  2:368419‑M2/39‑1 (MQ=255)
                        cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                           >  1:2002743‑M2/61‑99 (MQ=255)
                         ggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAtt                                                                          <  2:1786399‑M2/40‑1 (MQ=255)
                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTg                                                                         <  1:1568218‑M2/41‑1 (MQ=255)
                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  1:364983‑M2/59‑99 (MQ=255)
                           ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  1:2170010‑M2/42‑1 (MQ=255)
                           ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  2:804349‑M2/42‑1 (MQ=255)
                             gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                                      <  1:778117‑M2/44‑1 (MQ=255)
                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     <  1:1281738‑M2/45‑1 (MQ=255)
                               gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                                    >  1:651259‑M2/54‑99 (MQ=255)
                               gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                                    >  1:1518044‑M2/54‑99 (MQ=255)
                                aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   >  1:1729078‑M2/53‑99 (MQ=255)
                                  gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 >  1:1838106‑M2/51‑99 (MQ=255)
                                  gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 <  1:2097597‑M2/49‑1 (MQ=255)
                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                <  2:1143502‑M2/50‑1 (MQ=255)
                                   aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                <  1:477368‑M2/50‑1 (MQ=255)
                                    aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               >  2:599818‑M2/49‑99 (MQ=255)
                                      aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                             <  2:2122153‑M2/53‑1 (MQ=255)
                                      aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                             >  2:1448591‑M2/47‑99 (MQ=255)
                                       agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            >  1:842909‑M2/46‑99 (MQ=255)
                                       agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            >  2:488342‑M2/46‑99 (MQ=255)
                                        gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                           <  2:1937254‑M2/55‑1 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          <  1:1360739‑M2/56‑1 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          <  2:1338882‑M2/56‑1 (MQ=255)
                                           cAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                                         <  1:106706/99‑1 (MQ=255)
                                           aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        >  1:2061522‑M2/42‑99 (MQ=255)
                                           aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        >  2:105173‑M2/42‑99 (MQ=255)
                                            atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                                       <  2:674752‑M2/59‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  2:1634415‑M2/60‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  1:1870646‑M2/60‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  1:1196378‑M2/60‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  1:1775535‑M2/60‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      >  2:870856‑M2/40‑99 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  1:2205881‑M2/60‑1 (MQ=255)
                                              gatgtagccgtcaagttgtcataatttgtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcac                                                     >  1:1104502‑M2/39‑99 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    <  1:926469‑M2/62‑1 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    >  2:1985058‑M2/38‑99 (MQ=255)
                                                tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   <  2:764521‑M2/63‑1 (MQ=255)
                                                tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   >  2:1530146‑M2/37‑99 (MQ=255)
                                                tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   >  2:249071‑M2/37‑99 (MQ=255)
                                                  tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                                 >  1:287127‑M2/35‑99 (MQ=255)
                                                    gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               >  1:1035109‑M2/33‑99 (MQ=255)
                                                      cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             >  2:452960‑M2/31‑99 (MQ=255)
                                                      cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             <  1:2140241‑M2/69‑1 (MQ=255)
                                                        tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  2:2199468‑M2/29‑99 (MQ=255)
                                                        tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  1:1242440‑M2/29‑99 (MQ=255)
                                                          aagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                                         <  2:1105211‑M2/73‑1 (MQ=255)
                                                          aagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAg                                         <  1:1438257‑M2/73‑1 (MQ=255)
                                                            gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       >  2:1401993‑M2/25‑99 (MQ=255)
                                                             ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                                      <  2:1823048‑M2/76‑1 (MQ=255)
                                                               gtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGa                                    >  1:1644010‑M2/22‑99 (MQ=255)
                                                                tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                                   >  2:294904‑M2/21‑99 (MQ=255)
                                                                      ttggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                             >  1:881564‑M2/15‑99 (MQ=255)
                                                                       tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                            <  1:730392‑M2/86‑1 (MQ=255)
                                                                        ggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                           <  2:3377‑M2/87‑1 (MQ=255)
                                                                         gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                          >  1:139573‑M2/12‑99 (MQ=255)
                                                                         gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                          <  1:1601729‑M2/88‑1 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:2048932‑M2/89‑1 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:634264‑M2/89‑1 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:1297710‑M2/89‑1 (MQ=255)
                                                                           accgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                        >  1:556607‑M2/10‑99 (MQ=255)
                                                                           accgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                        <  2:1372152‑M2/90‑1 (MQ=255)
                                                                            ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                       <  2:2075602‑M2/91‑1 (MQ=255)
                                                                               agctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    >  2:718230‑M2/6‑99 (MQ=255)
                                                                                gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  2:965458‑M2/95‑1 (MQ=255)
                                                                                gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  2:1488282‑M2/95‑1 (MQ=255)
                                                                                gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  1:43755‑M2/95‑1 (MQ=255)
                                                                                   cGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg                <  2:1013780‑M2/98‑1 (MQ=255)
                                                                                     aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  2:414695/1‑99 (MQ=255)
                                                                                      attcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             <  2:2022811/96‑1 (MQ=255)
                                                                                       ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  1:250381/3‑99 (MQ=255)
                                                                                       ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  2:1901859/3‑99 (MQ=255)
                                                                                       ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  2:561775/3‑99 (MQ=255)
                                                                                       ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  2:169918/97‑1 (MQ=255)
                                                                                        ccGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  1:1147103/99‑1 (MQ=255)
                                                                                          gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         >  1:692724/1‑99 (MQ=255)
                                                                                          gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         <  2:156876/99‑1 (MQ=255)
                                                                                          gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         <  2:1737266/99‑1 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  1:398090/1‑99 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:840011/1‑99 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  1:329339/1‑99 (MQ=255)
                                                                                             aTCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTa       <  1:513131/99‑1 (MQ=255)
                                                                                             aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      <  1:975993/99‑1 (MQ=255)
                                                                                             aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      <  2:1599229/99‑1 (MQ=255)
                                                                                             aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  1:1918332/1‑99 (MQ=255)
                                                                                                  |                                                                                                  
GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398

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: