Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,238,103 T→C K396R (AAA→AGA)  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,1030TC99.2% 461.3 / NA 129K396R (AAA→AGA) mglAmethyl‑galactoside ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (66/62);  minor base A (1/0);  total (67/62)
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 129 total aligned reads displayed.

AATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATATTGGA  >  NC_000913/2238007‑2238197
                                                                                                |                                                                                              
aaTATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTa                                                                                             >  1:323351/1‑100 (MQ=255)
 atatTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTAc                                                                                            <  1:1314937/100‑1 (MQ=255)
  tatTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTAcc                                                                                           <  2:194912/100‑1 (MQ=255)
   atTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTAccc                                                                                          >  2:2278041/1‑100 (MQ=255)
        tGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATc                                                                                     <  2:2521816/100‑1 (MQ=255)
         ggTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCg                                                                                    <  1:767316/100‑1 (MQ=255)
           ttGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGCCCGATGACCCGGCGTTCTTACCCGCATCGaa                                                                                  >  1:1638926/1‑100 (MQ=255)
            tGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaat                                                                                 >  1:2093925/1‑100 (MQ=255)
            tGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaat                                                                                 >  1:1980644/1‑100 (MQ=255)
            tGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaat                                                                                 <  2:1974226/100‑1 (MQ=255)
            tGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaat                                                                                 <  2:1351623/100‑1 (MQ=255)
             gatctagtagCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatc                                                                                <  2:1764088/96‑1 (MQ=255)
             gCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatc                                                                                <  2:1091326/100‑1 (MQ=255)
              cGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatca                                                                               >  2:1056002/1‑100 (MQ=255)
              cGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCCaatc                                                                                >  1:509454/1‑99 (MQ=255)
              cGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCACATCGaatc                                                                                >  1:134502/1‑99 (MQ=255)
               gTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatcaa                                                                              >  1:1891828/1‑100 (MQ=255)
               gTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatcaa                                                                              >  1:1036438/1‑100 (MQ=255)
               gTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatcaa                                                                              >  1:2259997/1‑100 (MQ=255)
                 tagtagCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatcaatc                                                                            >  1:759072/1‑100 (MQ=255)
                  agtagCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGaatcagatc                                                                           >  1:1764062/1‑96 (MQ=255)
                   gtagCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCAc                                                                          >  2:726885/1‑100 (MQ=255)
                     agCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCAccc                                                                        <  1:924897/100‑1 (MQ=255)
                      gCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCa                                                                       <  1:1007393/100‑1 (MQ=255)
                        naGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACt                                                                     >  2:832548/2‑100 (MQ=255)
                         aGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTg                                                                    >  2:2647881/1‑100 (MQ=255)
                          gCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTgg                                                                   <  1:1451720/100‑1 (MQ=255)
                          gCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTgg                                                                   <  2:1766391/100‑1 (MQ=255)
                           cGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTggg                                                                  <  2:470831/100‑1 (MQ=255)
                            gACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGt                                                                 <  2:21328/100‑1 (MQ=255)
                             aCCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTa                                                                <  2:764381/100‑1 (MQ=255)
                             aCCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCAACCACTg                                                                    >  1:2836782/1‑96 (MQ=255)
                              ccAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTAt                                                               >  2:1643290/1‑100 (MQ=255)
                              ccAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTCCGTCCGATTACCCGGCGTTCTTACCCGCATCGCATCAATCACCCACTGGGTAt                                                               >  1:1573033/1‑100 (MQ=255)
                               cAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATc                                                              <  2:2826201/100‑1 (MQ=255)
                                aataatCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTAttg                                                             >  1:1290160/1‑98 (MQ=255)
                                aataatCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCg                                                             >  2:2449720/1‑100 (MQ=255)
                                  taatCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCt                                                           <  2:367935/100‑1 (MQ=255)
                                   aatCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCtt                                                          <  2:1166830/100‑1 (MQ=255)
                                     tCACCTTTTGCTGATTACTACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCtttt                                                        <  2:2768785/100‑1 (MQ=255)
                                     tCACCTTTTGCTGATTACCCCCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCtttt                                                        <  2:2906591/100‑1 (MQ=255)
                                     tCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCtttt                                                        <  2:132420/100‑1 (MQ=255)
                                      cACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTGCCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTc                                                       <  1:2727646/100‑1 (MQ=255)
                                      cACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTc                                                       <  1:2094799/100‑1 (MQ=255)
                                       aCCTTTTGCTGATTACCACCGGAGCGCGAACCAATCTGCGTCCGATGACCCGGCGTTCTTCCTCGCATCGAATCAATCACCCACTGGGTAACGCTTTTc                                                       >  1:2551294/1‑99 (MQ=255)
                                        ccTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCCGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCAt                                                     >  1:916949/1‑100 (MQ=255)
                                         cTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCa                                                      <  1:541013/98‑1 (MQ=255)
                                             tGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCcgcg                                                <  2:1329907/100‑1 (MQ=255)
                                                 gATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCAAATGGGTATCGCTTTTAATAAGAGAGtt                                            >  2:1985317/1‑100 (MQ=255)
                                                     accaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATcc                                        <  1:2870882/100‑1 (MQ=255)
                                                     accaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATcc                                        <  2:1290178/100‑1 (MQ=255)
                                                     accaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATcc                                        <  1:289791/100‑1 (MQ=255)
                                                     accaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATc                                         >  2:1504459/1‑99 (MQ=255)
                                                      ccaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCa                                       >  1:1643842/1‑100 (MQ=255)
                                                      ccaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCa                                       >  2:300286/1‑100 (MQ=255)
                                                      ccaccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCa                                       >  1:72022/1‑100 (MQ=255)
                                                        accGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGt                                     >  2:2892029/1‑100 (MQ=255)
                                                        accGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGt                                     <  1:2018405/100‑1 (MQ=255)
                                                        accGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGt                                     <  1:1169767/100‑1 (MQ=255)
                                                        accGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGt                                     <  2:1036451/100‑1 (MQ=255)
                                                        accGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGt                                     >  1:1040934/1‑100 (MQ=255)
                                                         ccGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTa                                    >  2:1359275/1‑100 (MQ=255)
                                                          cGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTaa                                   >  2:2691336/1‑100 (MQ=255)
                                                          cGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTaa                                   >  2:727358/1‑100 (MQ=255)
                                                            gagaGCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAAc                                 <  1:2251747/100‑1 (MQ=255)
                                                              gagCGAACCAATTTGCGTCCGATGGCCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCa                               <  2:1835679/100‑1 (MQ=255)
                                                              gagCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCa                               <  1:1475979/100‑1 (MQ=255)
                                                              gagCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCa                               <  2:1602311/100‑1 (MQ=255)
                                                                gCGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAAc                             >  2:81527/1‑100 (MQ=255)
                                                                 cGAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACt                            <  1:548012/100‑1 (MQ=255)
                                                                  gAACCAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACtt                           >  1:1973808/1‑100 (MQ=255)
                                                                     ccAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAt                        <  1:2392153/100‑1 (MQ=255)
                                                                     ccAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAt                        >  2:1003076/1‑100 (MQ=255)
                                                                     ccAATTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAt                        >  1:2157608/1‑100 (MQ=255)
                                                                       aaTTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAttt                      <  1:81525/100‑1 (MQ=255)
                                                                       aaTTTGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAttt                      <  2:2208138/100‑1 (MQ=255)
                                                                         tttGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTAttttt                    >  2:623265/1‑100 (MQ=255)
                                                                          ttGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTg                   <  1:1019500/100‑1 (MQ=255)
                                                                          ttGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTg                   <  1:2849652/100‑1 (MQ=255)
                                                                          ttGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTg                   <  2:2191148/100‑1 (MQ=255)
                                                                          ttGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTg                   <  1:1343864/100‑1 (MQ=255)
                                                                           tGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGt                  <  2:2378595/100‑1 (MQ=255)
                                                                           tGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGt                  >  2:2352896/1‑100 (MQ=255)
                                                                           tGCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGt                  <  2:2888829/100‑1 (MQ=255)
                                                                            gCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTa                 >  2:1558470/1‑100 (MQ=255)
                                                                            gCGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTa                 >  1:1856726/1‑100 (MQ=255)
                                                                             cGTCCGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAg                <  1:2427939/100‑1 (MQ=255)
                                                                                ccGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGtt              >  2:2198753/1‑99 (MQ=255)
                                                                                ccGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTg             >  1:2073648/1‑100 (MQ=255)
                                                                                 cGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGc            >  1:137131/1‑100 (MQ=255)
                                                                                 cGATGACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGc            >  1:1992088/1‑100 (MQ=255)
                                                                                     gACCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAat        <  1:2285194/100‑1 (MQ=255)
                                                                                      aCCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAata       >  2:2609529/1‑100 (MQ=255)
                                                                                      aCCCGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAata       >  1:1747789/1‑100 (MQ=255)
                                                                                       cccGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAatat      >  1:1968135/1‑100 (MQ=255)
                                                                                       cccGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAatat      >  1:1170845/1‑100 (MQ=255)
                                                                                        ccGGCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATAtt     <  2:411454/100‑1 (MQ=255)
                                                                                          ggCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATATTgg   <  2:2260024/100‑1 (MQ=255)
                                                                                          ggCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATATTgg   <  1:339676/100‑1 (MQ=255)
                                                                                           gCGTTCTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATATTGGa  <  2:1625212/100‑1 (MQ=255)
                                                                                                |                                                                                              
AATATTTCTGGTTGCGTTAGTAGCCAGCGACCAATAATCACCTTTTGCTGATTACCACCGGAGAGCGAACCAATTTGCGTCCGATGACCCGGCGTTTTTACCCGCATCGAATCAATCACCCACTGGGTATCGCTTTTCATCCGCGAGTTATCCAGTAAACCAACTTTATTTTTGTAGTTGCGAATATTGGA  >  NC_000913/2238007‑2238197

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