Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,238,656 T→A 5.8% I212F (ATC→TTC)  mglA ← methyl‑galactoside ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,238,6560TA 5.8% 362.1 / 10.8 103I212F (ATC→TTC) mglAmethyl‑galactoside ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (45/52);  new base A (3/3);  total (48/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

Only 100 of 108 total aligned reads displayed.

TTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGC  >  NC_000913/2238510‑2238803
                                                                                                                                                  |                                                                                                                                                   
tttgtttTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAttt                                                                                                                                                   >  4:1325242/1‑149 (MQ=255)
  tgtttTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttct                                                                                                                                                 <  2:622345/149‑1 (MQ=255)
   gtttTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                <  3:150427/149‑1 (MQ=255)
   gtttTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                <  1:278956/149‑1 (MQ=255)
      ttCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttc                                                                                                                                                  <  3:927689/144‑1 (MQ=255)
      ttCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCa                                                                                                                                             >  3:296286/1‑149 (MQ=255)
            gTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTAt                                                                                                                                       >  3:468137/1‑149 (MQ=255)
             tCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATg                                                                                                                                      >  1:1381627/1‑149 (MQ=255)
               aGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAttt                                                                                                                                                   <  4:816789/134‑1 (MQ=255)
               aGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAttt                                                                                                                                                   >  3:816789/1‑134 (MQ=255)
                gggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTTCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCga                                                                                                                                   <  3:1013365/149‑1 (MQ=255)
                 ggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCgag                                                                                                                                  <  1:872431/149‑1 (MQ=255)
                 ggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCgag                                                                                                                                  <  4:296286/149‑1 (MQ=255)
                 ggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCgag                                                                                                                                  >  2:1055589/1‑149 (MQ=255)
                 ggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAAATCTTCCa                                                                                                                                             >  1:960524/1‑138 (MQ=255)
                 ggAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAAATCTTCCa                                                                                                                                             <  2:960524/138‑1 (MQ=255)
                  gAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCgaga                                                                                                                                 <  1:477995/149‑1 (MQ=255)
                   aaaCGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTAt                                                                                                                                       >  1:134724/1‑142 (MQ=255)
                   aaaCGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTAt                                                                                                                                       <  2:134724/142‑1 (MQ=255)
                    aaCGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGata                                                                                                                               <  4:795033/149‑1 (MQ=255)
                     aCGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGatat                                                                                                                              <  2:1381627/149‑1 (MQ=255)
                      cGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGAtata                                                                                                                             >  1:731384/1‑149 (MQ=255)
                           gTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGGCCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCAt                                                                                                                                            <  3:69582/129‑1 (MQ=255)
                           gTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGGCCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCAt                                                                                                                                            >  4:69582/1‑129 (MQ=255)
                                gagcACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATAc                                                                                                                     >  3:150634/1‑147 (MQ=255)
                                gagaACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttc                                                                                                                                                  <  3:619344/118‑1 (MQ=255)
                                gagaACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttc                                                                                                                                                  >  4:619344/1‑118 (MQ=255)
                                gagaACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATAcc                                                                                                                    >  4:880592/1‑148 (MQ=255)
                                gagaACGCCCAACCATCATGCCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATAc                                                                                                                     <  4:150634/147‑1 (MQ=255)
                                  gaACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCa                                                                                                                 >  2:353187/1‑149 (MQ=255)
                                   aaCGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAg                                                                                                                <  4:493688/149‑1 (MQ=255)
                                      gCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCcg                                                                                                             >  2:871146/1‑149 (MQ=255)
                                       cccAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCcgc                                                                                                            <  1:1055589/149‑1 (MQ=255)
                                        ccAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAt                                                                                                                                                     >  1:1415308/1‑107 (MQ=255)
                                          aaCCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                <  1:295256/110‑1 (MQ=255)
                                          aaCCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                >  2:295256/1‑110 (MQ=255)
                                          aaCCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGct                                                                                                         >  3:1350132/1‑149 (MQ=255)
                                              atcatGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCga                                                                                                                                   >  2:1132604/1‑119 (MQ=255)
                                              atcatGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCga                                                                                                                                   <  1:1132604/119‑1 (MQ=255)
                                              atcatGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCttt                                                                                                     >  3:526393/1‑149 (MQ=255)
                                                  tGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAAt                                                                                                 <  1:353187/149‑1 (MQ=255)
                                                    gCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAAttt                                                                                               <  4:526393/149‑1 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                <  1:1089551/98‑1 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATttctt                                                                                                                                                >  2:1089551/1‑98 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATAc                                                                                                                     <  3:650908/125‑1 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATAc                                                                                                                     >  4:650908/1‑125 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTAc                                                                                             >  3:465863/1‑149 (MQ=255)
                                                      gatgatCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTAc                                                                                             >  4:855932/1‑149 (MQ=255)
                                                             ttGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGaataat                                                                                      >  1:904843/1‑149 (MQ=255)
                                                                tCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTaa                                                                                                  >  3:1298007/1‑134 (MQ=255)
                                                                tCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTaa                                                                                                  <  4:1298007/134‑1 (MQ=255)
                                                                tCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGt                                                                                   <  3:1325242/149‑1 (MQ=255)
                                                                   aTCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTa                                                                                                                                        <  4:264629/93‑1 (MQ=255)
                                                                   aTCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTa                                                                                                                                        >  3:264629/1‑93 (MQ=255)
                                                                   aTCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAAATCTTCCa                                                                                                                                             <  1:162984/88‑1 (MQ=255)
                                                                   aTCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGAAATCTTCCa                                                                                                                                             >  2:162984/1‑88 (MQ=255)
                                                                      gTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCa                                                                                                                                             >  4:724704/1‑85 (MQ=255)
                                                                      gTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCa                                                                                                                                             <  3:724704/85‑1 (MQ=255)
                                                                       tCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTaa                                                                                                  <  3:850244/127‑1 (MQ=255)
                                                                       tCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTaa                                                                                                  >  4:850244/1‑127 (MQ=255)
                                                                       tCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTg                                                                                  <  2:401368/143‑1 (MQ=255)
                                                                       tCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTg                                                                                  >  1:401368/1‑143 (MQ=255)
                                                                       tCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAga                                                                            <  4:465863/149‑1 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTATTTTACGaataat                                                                                      >  1:217169/1‑135 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAAt                                                                                                 <  1:1020708/124‑1 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAAt                                                                                                 >  2:1020708/1‑124 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTAc                                                                                             >  3:1394695/1‑128 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTAc                                                                                             <  4:1394695/128‑1 (MQ=255)
                                                                           tCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGaataat                                                                                      <  2:217169/135‑1 (MQ=255)
                                                                                       ttCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCttttt                                                            >  3:360145/1‑149 (MQ=255)
                                                                                         cGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCg                                                          >  4:1100108/1‑149 (MQ=255)
                                                                                              gCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTaa                                                     <  4:360145/149‑1 (MQ=255)
                                                                                                gATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACg                                                   >  1:1005785/1‑149 (MQ=255)
                                                                                                 aTCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCga                                                                                                                                   >  3:754348/1‑68 (MQ=255)
                                                                                                 aTCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCga                                                                                                                                   <  4:754348/68‑1 (MQ=255)
                                                                                                 aTCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACga                                                  >  1:851545/1‑149 (MQ=255)
                                                                                                      cTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCg                                                          <  1:102394/136‑1 (MQ=255)
                                                                                                      cTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCg                                                          >  2:102394/1‑136 (MQ=255)
                                                                                                       tGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCg                                                          >  1:422123/1‑135 (MQ=255)
                                                                                                       tGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCg                                                          <  2:422123/135‑1 (MQ=255)
                                                                                                          ccGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCgg                                         <  2:904843/149‑1 (MQ=255)
                                                                                                           cGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACctc                                                                 >  4:341233/1‑124 (MQ=255)
                                                                                                           cGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGt                                        <  4:1350132/149‑1 (MQ=255)
                                                                                                             tCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGt                                                                                   <  3:314571/104‑1 (MQ=255)
                                                                                                             tCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGt                                                                                   >  4:314571/1‑104 (MQ=255)
                                                                                                              cgcgCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCa                                     <  3:855932/149‑1 (MQ=255)
                                                                                                                 gcAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATcc                                  <  2:731384/149‑1 (MQ=255)
                                                                                                                                 atcaCATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTcg                  >  1:1030045/1‑149 (MQ=255)
                                                                                                                                  tcaCATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTcgc                 <  2:1030045/149‑1 (MQ=255)
                                                                                                                                      aTAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACTATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGAcc                                                                   <  1:308964/95‑1 (MQ=255)
                                                                                                                                      aTAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGAcc                                                                   >  2:308964/1‑95 (MQ=255)
                                                                                                                                      aTAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATaa                              >  2:1440196/1‑132 (MQ=255)
                                                                                                                                      aTAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATaa                              <  1:1440196/132‑1 (MQ=255)
                                                                                                                                         aatGGAAGAAATCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACgaagaa                                              >  1:842205/3‑113 (MQ=255)
                                                                                                                                         aCTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAg          <  2:194048/149‑1 (MQ=255)
                                                                                                                                           tGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCgcgc                                                                                                          >  2:101661/1‑51 (MQ=255)
                                                                                                                                           tGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCgcgc                                                                                                          <  1:101661/51‑1 (MQ=255)
                                                                                                                                            ggAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTcgc                 <  2:739735/139‑1 (MQ=255)
                                                                                                                                            ggAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTcgc                 >  1:739735/1‑139 (MQ=255)
                                                                                                                                                 aTTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGc  <  2:1005785/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
TTTGTTTTCTTTGTCAGGGAAACGCTGGTTAAGAGAACGCCCAACCATCATGGCGATGATCTTGTCCATCGTCAGTCCTGCCAGCGGTTCGGTGGCGATCCACTGACCGTCGCGCAATACGGTAACTTCATCACATAACTGGAAGATTTCTTCCATTTTATGCGAGATATAAACAATACCGCAGCCGCGCTCTTTTAATTTACGAATAATAGTGAACAGATGATTGACCTCTTTTTCGGTTAACGAAGAAGTCGGTTCATCCATAATCACAATTTTCGCGTTATAGGAAAACGC  >  NC_000913/2238510‑2238803

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