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.99.7% 194478.7 / inf 56659coding (6/3075 nt)lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base G (67/68);  new base . (28751/27749);  total (28827/27832)
Fisher's exact test for biased strand distribution p-value = 7.96e-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 59781 total aligned reads displayed.

CGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366203‑366398
                                                                                                 |                                                                                                  
  cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                                 <  2:2512890‑M2/18‑4 (MQ=255)
    gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                              <  1:371550‑M2/20‑1 (MQ=255)
    gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                              >  1:1929498‑M2/80‑99 (MQ=255)
     cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                             >  2:2792734‑M2/79‑99 (MQ=255)
      ccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGc                                                                                            <  1:1216460‑M2/22‑1 (MQ=255)
      ccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGc                                                                                            <  1:1268938‑M2/22‑1 (MQ=255)
       cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                           <  2:2848997‑M2/23‑1 (MQ=255)
        gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                                          >  1:2739231‑M2/76‑99 (MQ=255)
         agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         <  2:2904613‑M2/25‑1 (MQ=255)
         agcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGt                                                                                         >  1:2042358‑M2/75‑99 (MQ=255)
          gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                                        <  2:1251441‑M2/26‑1 (MQ=255)
           cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGttt                                                                                       <  2:83429‑M2/27‑1 (MQ=255)
             agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTcc                                                                                     <  2:1635372‑M2/29‑1 (MQ=255)
                 ccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtg                                                                                 <  2:1693100‑M2/33‑1 (MQ=255)
                  cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                                >  2:1855212‑M2/66‑99 (MQ=255)
                  cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                                <  1:319373‑M2/34‑1 (MQ=255)
                     gccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaa                                                                             <  1:1051927‑M2/37‑1 (MQ=255)
                     gccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaa                                                                             >  1:538378‑M2/63‑99 (MQ=255)
                      ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                            <  2:1798610‑M2/38‑1 (MQ=255)
                      ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                            <  1:877036‑M2/38‑1 (MQ=255)
                       cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                           >  1:2152148‑M2/61‑99 (MQ=255)
                          ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  1:93102‑M2/42‑1 (MQ=255)
                          ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        >  2:1419922‑M2/58‑99 (MQ=255)
                          ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  2:449892‑M2/42‑1 (MQ=255)
                          ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        >  1:1325238‑M2/58‑99 (MQ=255)
                          ttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGt                                                                        <  2:1910715‑M2/42‑1 (MQ=255)
                           tgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGtt                                                                       <  1:653919‑M2/43‑1 (MQ=255)
                            gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                                      <  1:2466895‑M2/44‑1 (MQ=255)
                            gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                                      <  2:1144939‑M2/44‑1 (MQ=255)
                             tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     <  2:393167‑M2/45‑1 (MQ=255)
                             tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                                     >  2:2626946‑M2/55‑99 (MQ=255)
                               aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                                   <  1:2099534‑M2/47‑1 (MQ=255)
                                agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                                  <  2:1197844‑M2/48‑1 (MQ=255)
                                 gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 <  1:1599768‑M2/49‑1 (MQ=255)
                                 gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 <  1:1108816‑M2/49‑1 (MQ=255)
                                 gaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGc                                                                 <  2:1527706‑M2/49‑1 (MQ=255)
                                  aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                >  1:2257303‑M2/50‑99 (MQ=255)
                                  aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                <  1:1161544‑M2/50‑1 (MQ=255)
                                  aaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCt                                                                >  2:1475492‑M2/50‑99 (MQ=255)
                                   aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                               <  2:2679384‑M2/51‑1 (MQ=255)
                                     aagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                             >  2:1042590‑M2/47‑99 (MQ=255)
                                     aagtgaatgatgtaaccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcac                                                             >  2:2957511‑M2/47‑99 (MQ=255)
                                      agtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTcaca                                                            >  1:2791750‑M2/46‑99 (MQ=255)
                                       gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                           <  2:2350638‑M2/55‑1 (MQ=255)
                                        tgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAt                                                          >  1:2897521‑M2/44‑99 (MQ=255)
                                         gaatgatgtagccgtctagttctcataattggtaccgaggtcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAAtt                                                         <  2:1307978‑M2/57‑1 (MQ=255)
                                          aatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTc                                                        <  1:1761402‑M2/58‑1 (MQ=255)
                                           atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                                       >  1:1374233‑M2/41‑99 (MQ=255)
                                           atgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTcc                                                       >  1:2629080‑M2/41‑99 (MQ=255)
                                            tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  2:2294188‑M2/60‑1 (MQ=255)
                                            tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  1:762414‑M2/60‑1 (MQ=255)
                                            tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                                      <  2:338581‑M2/60‑1 (MQ=255)
                                                  agccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaac                                                >  2:724057‑M2/34‑99 (MQ=255)
                                                   gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:484758‑M2/67‑1 (MQ=255)
                                                   gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                               <  1:1000081‑M2/67‑1 (MQ=255)
                                                     cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             >  2:136783‑M2/31‑99 (MQ=255)
                                                     cgtcaagttgtcataatgcgtaccgagctcGAATTCGTAATCATGTCATAGCCGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                             <  1:1430756‑M2/69‑1 (MQ=255)
                                                      gtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            <  2:564611‑M2/70‑1 (MQ=255)
                                                      gtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            <  1:1487298‑M2/70‑1 (MQ=255)
                                                      gtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            >  1:182033‑M2/30‑99 (MQ=255)
                                                      gtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATAc                                            <  1:2684314‑M2/70‑1 (MQ=255)
                                                       tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  2:66885‑M2/29‑99 (MQ=255)
                                                       tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  1:819029‑M2/29‑99 (MQ=255)
                                                       tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                           >  1:2843856‑M2/29‑99 (MQ=255)
                                                          agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                                        >  2:2404040‑M2/26‑99 (MQ=255)
                                                          agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                                        >  1:1705513‑M2/26‑99 (MQ=255)
                                                          agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                                        <  2:2735340‑M2/74‑1 (MQ=255)
                                                           gttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGcc                                       <  1:427007‑M2/75‑1 (MQ=255)
                                                            ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                                      >  2:944208‑M2/24‑99 (MQ=255)
                                                               tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                                   >  2:2568860‑M2/21‑99 (MQ=255)
                                                                 ataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGc                                 >  1:300009‑M2/19‑99 (MQ=255)
                                                                   aattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCAt                               <  2:1584318‑M2/83‑1 (MQ=255)
                                                                     ttggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                             >  1:2202820‑M2/15‑99 (MQ=255)
                                                                         taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg                         <  2:3018104‑M2/89‑1 (MQ=255)
                                                                          accgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                        >  1:2361758‑M2/10‑99 (MQ=255)
                                                                          accgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtgt                        >  2:2456120‑M2/10‑99 (MQ=255)
                                                                           ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                       <  1:1987553‑M2/91‑1 (MQ=255)
                                                                           ccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTa                       >  2:1092832‑M2/9‑99 (MQ=255)
                                                                            cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa                      <  2:2240042‑M2/92‑1 (MQ=255)
                                                                             gagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa                     <  2:2023485‑M2/93‑1 (MQ=255)
                                                                             gagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaaa                     >  2:2130797‑M2/7‑99 (MQ=255)
                                                                              agctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    <  2:776028‑M2/94‑1 (MQ=255)
                                                                              agctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg                    <  2:2186916‑M2/94‑1 (MQ=255)
                                                                               gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   <  1:988155‑M2/95‑1 (MQ=255)
                                                                               gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   >  2:889201‑M2/5‑99 (MQ=255)
                                                                               gctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGc                   >  1:2521457‑M2/5‑99 (MQ=255)
                                                                                  cGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg                <  1:592457‑M2/98‑1 (MQ=255)
                                                                                    aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  1:1653455/1‑99 (MQ=255)
                                                                                    aaTTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTggg              >  2:1965588/1‑99 (MQ=255)
                                                                                      ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            <  1:1544236/97‑1 (MQ=255)
                                                                                      ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  1:2309203/3‑99 (MQ=255)
                                                                                      ttcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGt            >  2:616451/3‑99 (MQ=255)
                                                                                       tcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTg           <  2:2077601/98‑1 (MQ=255)
                                                                                       tcGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTg           <  1:709036/98‑1 (MQ=255)
                                                                                        cGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGc          <  1:2574338/99‑1 (MQ=255)
                                                                                         gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         >  1:2949024/1‑99 (MQ=255)
                                                                                         gTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGcc         <  2:1746806/99‑1 (MQ=255)
                                                                                          tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  2:507738/1‑99 (MQ=255)
                                                                                          tAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCt        >  1:1877415/1‑99 (MQ=255)
                                                                                            aTCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTaa      >  1:2414077/1‑99 (MQ=255)
                                                                                                 |                                                                                                  
CGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGGTGCCTAATGAG  >  NC_000913/366203‑366398

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