Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 366,301 Δ1 bp coding (5/3075 nt) lacZ ← beta‑galactosidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913366,3000G.97.7% 21869.6 / inf 5060coding (6/3075 nt)lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (59/58);  new base . (2503/2439);  total (2563/2497)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 5340 total aligned reads displayed.

GCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366202‑366398
                                                                                                  |                                                                                                  
ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                                    <  1:2291665‑M2/15‑1 (MQ=255)
ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                                    >  2:1613276‑M2/85‑99 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   <  1:2108571‑M2/16‑2 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   >  1:1097621‑M2/84‑98 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   <  1:44575‑M2/16‑2 (MQ=255)
 gccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAtgt                                                                                                   <  1:144061‑M2/16‑2 (MQ=255)
     gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                              >  1:193076‑M2/80‑99 (MQ=255)
     gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                              >  1:1482342‑M2/80‑99 (MQ=255)
       ccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGc                                                                                            <  1:2231238‑M2/22‑1 (MQ=255)
        cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                           >  2:792449‑M2/77‑99 (MQ=255)
         gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                                          <  1:2475785‑M2/24‑1 (MQ=255)
          agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         >  2:1282760‑M2/75‑99 (MQ=255)
          agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         <  1:616541‑M2/25‑1 (MQ=255)
          agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         >  2:2243593‑M2/75‑99 (MQ=255)
           gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                                        <  1:1050257‑M2/26‑1 (MQ=255)
            cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGttt                                                                                       >  2:1740636‑M2/73‑99 (MQ=255)
             gagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTc                                                                                      <  2:1649689‑M2/28‑1 (MQ=255)
              agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTcc                                                                                     <  1:706919‑M2/29‑1 (MQ=255)
               gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCt                                                                                    >  1:992111‑M2/70‑99 (MQ=255)
               gttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCt                                                                                    >  2:1697420‑M2/70‑99 (MQ=255)
                ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                                   >  2:1792818‑M2/69‑99 (MQ=255)
                 tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                                  <  2:1132491‑M2/32‑1 (MQ=255)
                   cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                                >  1:2387751‑M2/66‑99 (MQ=255)
                      gccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaa                                                                             >  2:321867‑M2/63‑99 (MQ=255)
                       ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                            >  1:1001957‑M2/62‑99 (MQ=255)
                       ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                            <  2:1034751‑M2/38‑1 (MQ=255)
                        cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                           >  1:398290‑M2/61‑99 (MQ=255)
                        cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                           <  2:847176‑M2/39‑1 (MQ=255)
                         ggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAtt                                                                          <  2:2142923‑M2/40‑1 (MQ=255)
                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  1:893876‑M2/59‑99 (MQ=255)
                          gttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTg                                                                         >  2:757617‑M2/59‑99 (MQ=255)
                           ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  1:2431985‑M2/42‑1 (MQ=255)
                             gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                                      <  1:1359244‑M2/44‑1 (MQ=255)
                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     <  1:266033‑M2/45‑1 (MQ=255)
                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     <  1:571736‑M2/45‑1 (MQ=255)
                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     >  2:1085617‑M2/55‑99 (MQ=255)
                              tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     >  1:2150519‑M2/55‑99 (MQ=255)
                               gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                                    <  1:2566766‑M2/46‑1 (MQ=255)
                               gaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATc                                                                    >  2:1941661‑M2/54‑99 (MQ=255)
                                aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   >  1:1180539‑M2/53‑99 (MQ=255)
                                aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   >  2:1997196‑M2/53‑99 (MQ=255)
                                aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   <  2:1326768‑M2/47‑1 (MQ=255)
                                    aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               >  2:2519367‑M2/49‑99 (MQ=255)
                                     aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                              >  2:2562399‑M2/48‑99 (MQ=255)
                                       agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            >  2:370882‑M2/46‑99 (MQ=255)
                                       agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            >  2:2391591‑M2/46‑99 (MQ=255)
                                        gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                           >  2:1518790‑M2/45‑99 (MQ=255)
                                        gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                           <  2:1418222‑M2/55‑1 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:1984124‑M2/44‑99 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          <  2:2090144‑M2/56‑1 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  2:1569568‑M2/44‑99 (MQ=255)
                                         tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          <  1:2184652‑M2/56‑1 (MQ=255)
                                          gaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                                         >  2:659567‑M2/43‑99 (MQ=255)
                                           aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        <  1:345219‑M2/58‑1 (MQ=255)
                                           aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        <  1:1137556‑M2/58‑1 (MQ=255)
                                             tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      >  2:420689‑M2/40‑99 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    >  1:2018912‑M2/38‑99 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    <  2:565547‑M2/62‑1 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    >  1:306285‑M2/38‑99 (MQ=255)
                                               atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                                    >  2:234145‑M2/38‑99 (MQ=255)
                                                tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                                   >  2:2293149‑M2/37‑99 (MQ=255)
                                                    gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:1582709‑M2/67‑1 (MQ=255)
                                                     ccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACAt                                              >  2:1367277‑M2/32‑99 (MQ=255)
                                                      cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             <  1:849388‑M2/69‑1 (MQ=255)
                                                      cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             <  2:1897265‑M2/69‑1 (MQ=255)
                                                         caagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGa                                          >  2:2188020‑M2/28‑99 (MQ=255)
                                                           agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                                        >  1:525245‑M2/26‑99 (MQ=255)
                                                            gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       <  2:2042043‑M2/75‑1 (MQ=255)
                                                            gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       <  1:946103‑M2/75‑1 (MQ=255)
                                                            gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       >  1:1685260‑M2/25‑99 (MQ=255)
                                                             ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                                      >  1:233823‑M2/24‑99 (MQ=255)
                                                             ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                                      <  1:1837730‑M2/76‑1 (MQ=255)
                                                                   taattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCa                                >  1:661537‑M2/18‑99 (MQ=255)
                                                                     attggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                              <  2:2471471‑M2/84‑1 (MQ=255)
                                                                       tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                            <  2:948108‑M2/86‑1 (MQ=255)
                                                                       tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                            <  1:1883861‑M2/86‑1 (MQ=255)
                                                                         gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt                          >  1:139510‑M2/12‑99 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  1:1304260‑M2/89‑1 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:515585‑M2/89‑1 (MQ=255)
                                                                          taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  1:1746719‑M2/89‑1 (MQ=255)
                                                                                 ctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc                  >  1:2515314‑M2/4‑99 (MQ=255)
                                                                                  tcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt                 >  1:1535436‑M2/3‑99 (MQ=255)
                                                                                  tcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt                 <  1:749265‑M2/97‑1 (MQ=255)
                                                                                     aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  1:1783203/1‑99 (MQ=255)
                                                                                     aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  1:483932/1‑99 (MQ=255)
                                                                                     aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              <  2:1469428/99‑1 (MQ=255)
                                                                                      attcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             <  2:464247/96‑1 (MQ=255)
                                                                                      attcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             <  2:602159/96‑1 (MQ=255)
                                                                                      attcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             <  2:1181225/96‑1 (MQ=255)
                                                                                      attcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgggg             >  1:1221663/4‑99 (MQ=255)
                                                                                       ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  2:1846759/97‑1 (MQ=255)
                                                                                          gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         <  2:904230/99‑1 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:274386/1‑99 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:70205/1‑99 (MQ=255)
                                                                                           tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:1235211/1‑99 (MQ=255)
                                                                                           tAATCATGACATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:2432782/1‑99 (MQ=255)
                                                                                            aaTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTa       >  1:1633661/1‑99 (MQ=255)
                                                                                            aaTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTa       <  1:1261582/99‑1 (MQ=255)
                                                                                             aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      <  1:206905/99‑1 (MQ=255)
                                                                                             aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  1:794750/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: