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.1% 31093.5 / 74.4 8831D6N (GAT→AAT) lacZbeta‑galactosidase
Reads supporting (aligned to +/- strand):  ref base C (43/38);  new base T (4157/4589);  total (4202/4629)
Fisher's exact test for biased strand distribution p-value = 3.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 8831 total aligned reads displayed.

GCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGG  >  NC_000913/366194‑366387
                                                                                                |                                                                                                 
tgcaccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCg                                                                                                 <  2:1043832‑M2/7‑1 (MQ=255)
  caccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTa                                                                                               >  2:1496127‑M2/91‑99 (MQ=255)
   accggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTaa                                                                                              <  1:1110660‑M2/10‑1 (MQ=255)
    ccggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAAt                                                                                             <  2:856876‑M2/11‑1 (MQ=255)
      ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                           <  1:93805‑M2/13‑1 (MQ=255)
      ggggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCa                                                                                           <  1:1148061‑M2/13‑1 (MQ=255)
       gggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCAt                                                                                          <  2:123466‑M2/14‑1 (MQ=255)
        ggccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATg                                                                                         <  2:415953‑M2/15‑1 (MQ=255)
          ccagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCatgtc                                                                                       >  2:302263‑M2/83‑97 (MQ=255)
          ccagcccgagcgagttccgcgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCatgtc                                                                                       >  2:1633975‑M2/83‑97 (MQ=255)
           cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                      <  1:1485313‑M2/18‑4 (MQ=255)
           cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                      <  2:166973‑M2/18‑4 (MQ=255)
           cagcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagcagtcaagttgtcataattggtaccgagctcGAATTCGTAATcatgtca                                                                                      <  1:1807961‑M2/18‑4 (MQ=255)
             gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                   <  2:1855369‑M2/20‑1 (MQ=255)
             gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                   >  1:1228434‑M2/80‑99 (MQ=255)
             gcccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATa                                                                                   >  1:601264‑M2/80‑99 (MQ=255)
              cccgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAg                                                                                  <  2:936785‑M2/21‑1 (MQ=255)
               ccgagcgagttacgtgccggttgcgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGc                                                                                 >  2:2035807‑M2/78‑99 (MQ=255)
                cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                >  2:636926‑M2/77‑99 (MQ=255)
                cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                <  1:410248‑M2/23‑1 (MQ=255)
                cgagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCt                                                                                >  1:904350‑M2/77‑99 (MQ=255)
                 gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                               <  2:1777498‑M2/24‑1 (MQ=255)
                 gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                               <  1:786415‑M2/24‑1 (MQ=255)
                 gagcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTg                                                                               >  2:1306536‑M2/76‑99 (MQ=255)
                   gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                             <  2:1996970‑M2/26‑1 (MQ=255)
                   gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                             <  1:1881305‑M2/26‑1 (MQ=255)
                   gcgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGtt                                                                             <  1:1376424‑M2/26‑1 (MQ=255)
                    cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGttt                                                                            <  2:1818542‑M2/27‑1 (MQ=255)
                    cgagttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGttt                                                                            <  1:1705156‑M2/27‑1 (MQ=255)
                      agttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTcc                                                                          <  1:1448042‑M2/29‑1 (MQ=255)
                        ttccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtg                                                                        >  1:1068658‑M2/69‑99 (MQ=255)
                         tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                       <  2:1422919‑M2/32‑1 (MQ=255)
                         tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                       <  1:335708‑M2/32‑1 (MQ=255)
                         tccgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgt                                                                       <  1:827581‑M2/32‑1 (MQ=255)
                           cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                     <  2:601349‑M2/34‑1 (MQ=255)
                           cgtgccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCtgtgt                                                                     >  2:226497‑M2/66‑99 (MQ=255)
                               ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                 <  2:1892011‑M2/38‑1 (MQ=255)
                               ccggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGaaa                                                                 <  2:674525‑M2/38‑1 (MQ=255)
                                cggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAt                                                                >  2:1458622‑M2/61‑99 (MQ=255)
                                 ggttgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAAtt                                                               >  2:1739157‑M2/60‑99 (MQ=255)
                                    tgtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGtt                                                            >  1:373369‑M2/57‑99 (MQ=255)
                                    tAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGt                                                             <  2:108117/99‑1 (MQ=255)
                                     gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                           <  2:196349‑M2/44‑1 (MQ=255)
                                     gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                           <  2:780419‑M2/44‑1 (MQ=255)
                                     gtgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTa                                                           >  2:454927‑M2/56‑99 (MQ=255)
                                      tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                          >  1:1000540‑M2/55‑99 (MQ=255)
                                      tgaagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTAt                                                          >  1:1635950‑M2/55‑99 (MQ=255)
                                        aagaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATcc                                                        <  1:1265724‑M2/47‑1 (MQ=255)
                                         agaaaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCg                                                       <  1:1412045‑M2/48‑1 (MQ=255)
                                            aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                    >  1:1987239‑M2/49‑99 (MQ=255)
                                            aaaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTc                                                    >  2:74681‑M2/49‑99 (MQ=255)
                                             aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                   >  1:560495‑M2/48‑99 (MQ=255)
                                             aaagtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTca                                                   >  1:679692‑M2/48‑99 (MQ=255)
                                                gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                <  1:133970‑M2/55‑1 (MQ=255)
                                                gtgaatgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACaa                                                <  1:1519783‑M2/55‑1 (MQ=255)
                                                     tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                           <  1:752609‑M2/60‑1 (MQ=255)
                                                     tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                           >  2:341613‑M2/40‑99 (MQ=255)
                                                     tgatgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCca                                           <  1:258612‑M2/60‑1 (MQ=255)
                                                       atgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcaca                                         <  2:1125970‑M2/62‑1 (MQ=255)
                                                        tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                        >  1:1596929‑M2/37‑99 (MQ=255)
                                                        tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                        >  2:1344612‑M2/37‑99 (MQ=255)
                                                        tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                        >  1:296176‑M2/37‑99 (MQ=255)
                                                        tgtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacac                                        >  2:378527‑M2/37‑99 (MQ=255)
                                                         gtagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCcacaca                                       >  1:19471‑M2/36‑99 (MQ=255)
                                                          tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                      <  1:22263‑M2/65‑1 (MQ=255)
                                                          tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                      <  2:1169876‑M2/65‑1 (MQ=255)
                                                          tagccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaa                                      <  1:977925‑M2/65‑1 (MQ=255)
                                                            gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                    <  2:1929041‑M2/67‑1 (MQ=255)
                                                            gccgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACacaaca                                    <  1:1672002‑M2/67‑1 (MQ=255)
                                                              cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                  >  2:1704409‑M2/31‑99 (MQ=255)
                                                              cgtcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATa                                  >  1:1034424‑M2/31‑99 (MQ=255)
                                                                tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                >  2:816662‑M2/29‑99 (MQ=255)
                                                                tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                >  2:1010447‑M2/29‑99 (MQ=255)
                                                                tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                >  1:2029145‑M2/29‑99 (MQ=255)
                                                                tcaagttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACg                                <  1:445118‑M2/71‑1 (MQ=255)
                                                                   agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                             <  2:1205320‑M2/74‑1 (MQ=255)
                                                                   agttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGc                             >  1:1560328‑M2/26‑99 (MQ=255)
                                                                     ttgtcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCg                           >  2:1530544‑M2/24‑99 (MQ=255)
                                                                        tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                        <  2:1569106‑M2/79‑1 (MQ=255)
                                                                        tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                        >  2:1383900‑M2/21‑99 (MQ=255)
                                                                        tcataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGaa                        <  1:1304921‑M2/79‑1 (MQ=255)
                                                                         cataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                       >  2:710540‑M2/20‑99 (MQ=255)
                                                                         cataattggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAg                       <  1:717234‑M2/80‑1 (MQ=255)
                                                                             attggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATa                   >  1:1743458‑M2/16‑99 (MQ=255)
                                                                              ttggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                  >  1:520181‑M2/15‑99 (MQ=255)
                                                                              ttggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaa                  <  1:480581‑M2/85‑1 (MQ=255)
                                                                               tggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATaaa                 <  2:896337‑M2/86‑1 (MQ=255)
                                                                                ggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                <  1:118893‑M2/87‑1 (MQ=255)
                                                                                ggtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAg                <  1:868453‑M2/87‑1 (MQ=255)
                                                                                 gtaccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgt               <  2:568162‑M2/88‑1 (MQ=255)
                                                                                  taccgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAgtg              >  2:1601921‑M2/11‑99 (MQ=255)
                                                                                     cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa           >  2:967323‑M2/8‑99 (MQ=255)
                                                                                     cgagctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTaa           <  2:493655‑M2/92‑1 (MQ=255)
                                                                                       agctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg         <  1:1780646‑M2/94‑1 (MQ=255)
                                                                                       agctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAg         <  2:260050‑M2/94‑1 (MQ=255)
                                                                                         ctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc       <  1:1916930‑M2/96‑1 (MQ=255)
                                                                                         ctcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGcc       >  1:63976‑M2/4‑99 (MQ=255)
                                                                                          tcGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCt      <  2:529417‑M2/97‑1 (MQ=255)
                                                                                           cGAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTg     <  2:1267678‑M2/98‑1 (MQ=255)
                                                                                            gAATTCGTAATCATGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTgg    >  1:1846955/1‑99 (MQ=255)
                                                                                                |                                                                                                 
GCCAGCTGGCGAAAGGGGGATGTGCTGCAAGGCGATTAAGTTGGGTAACGCCAGGGTTTTCCCAGTCACGACGTTGTAAAACGACGGCCAGTGAATCCGTAATCATGGTCATAGCTGTTTCCTGTGTGAAATTGTTATCCGCTCACAATTCCACACAACATACGAGCCGGAAGCATAAAGTGTAAAGCCTGGGG  >  NC_000913/366194‑366387

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: