Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,469,733 (T)8→7 coding (603/1332 nt) gtrS → serotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,469,7260T.96.1% 524.6 / 13.5 128coding (596/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (52/71);  total (56/72)
Fisher's exact test for biased strand distribution p-value = 1.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 151 total aligned reads displayed.

CAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTG  >  NC_000913/2469628‑2469821
                                                                                                  |                                                                                               
cAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCatattatatt                                                                                                >  1:2199083/1‑100 (MQ=255)
cAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCatattatatt                                                                                                >  2:31084/1‑100 (MQ=255)
    aaaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttt                                                                                            <  1:2636777/100‑1 (MQ=255)
    aaaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAAAATATAAAAATAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttt                                                                                            <  2:1029859/100‑1 (MQ=255)
     tctCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           <  1:2087713/98‑1 (MQ=255)
      aTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttg                                                                                          <  1:1129082/100‑2 (MQ=255)
      aTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttg                                                                                          <  2:2728710/100‑2 (MQ=255)
       tCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttga                                                                                         <  1:2650817/100‑3 (MQ=255)
       tCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttga                                                                                         >  1:1602304/1‑98 (MQ=255)
        cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttaga                                                                                        >  2:2087700/1‑97 (MQ=255)
         tATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                       <  2:986315/100‑5 (MQ=255)
          aTGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAg                                                                                     <  2:2217663/100‑1 (MQ=255)
           tGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGa                                                                                    <  1:2430579/100‑1 (MQ=255)
            gCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGac                                                                                   >  1:2723053/1‑100 (MQ=255)
             cTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGaca                                                                                  >  1:1569662/1‑100 (MQ=255)
                 tGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACaaa                                                                              <  2:1025840/100‑1 (MQ=255)
                  gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAAc                                                                             <  1:305486/100‑1 (MQ=255)
                  gCTATTTTTTTCGATGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACaaa                                                                              >  2:1791511/1‑99 (MQ=255)
                   tattttttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAAc                                                                             >  1:1562858/4‑100 (MQ=255)
                   cTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑ATAGACACAAAc                                                                             >  2:1731991/1‑99 (MQ=255)
                    tATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑ATTGACACTAACt                                                                            >  2:1310853/1‑99 (MQ=255)
                    tATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTa                                                                           >  1:358751/1‑100 (MQ=255)
                    tATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTa                                                                           >  2:2637120/1‑100 (MQ=255)
                     aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTaa                                                                          <  1:1617983/100‑1 (MQ=255)
                     aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTaa                                                                          >  1:1195209/1‑100 (MQ=255)
                     aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTaa                                                                          <  2:943941/100‑1 (MQ=255)
                       ttttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAatt                                                                        <  2:1433731/100‑1 (MQ=255)
                        tttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATCTTATA‑TTTTTTTGATAGACACAAACTAatt                                                                        >  2:110928/1‑99 (MQ=255)
                        tttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAatta                                                                       >  1:1373368/1‑100 (MQ=255)
                           ttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAAttatta                                                                    >  2:343302/1‑100 (MQ=255)
                            tCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTaa                                                                   <  1:1835850/100‑1 (MQ=255)
                             cGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTaaa                                                                  >  1:2288452/1‑100 (MQ=255)
                                tGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCt                                                               >  1:149739/1‑100 (MQ=255)
                                tGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCt                                                               >  1:2132554/1‑100 (MQ=255)
                                ngCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCt                                                               >  2:934872/2‑100 (MQ=255)
                                 gCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTg                                                              <  2:360790/100‑1 (MQ=255)
                                  cTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGt                                                             <  1:1231961/100‑1 (MQ=255)
                                       aataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAatat                                                        >  2:419855/1‑100 (MQ=255)
                                       aataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAatat                                                        >  2:651800/1‑100 (MQ=255)
                                       aataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAatat                                                        >  2:2225494/1‑100 (MQ=255)
                                          atCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCg                                                     <  2:2198620/100‑1 (MQ=255)
                                          atCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCg                                                     <  2:1602294/100‑1 (MQ=255)
                                           tCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgt                                                    <  1:1728581/100‑1 (MQ=255)
                                           cgatcTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑‑TTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   <  1:825529/96‑1 (MQ=255)
                                            cTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   <  1:651802/100‑1 (MQ=255)
                                             tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtgg                                                  <  2:2024128/100‑1 (MQ=255)
                                              aCTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggt                                                 <  1:2394656/100‑1 (MQ=255)
                                               cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                <  2:2392587/100‑1 (MQ=255)
                                               cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                <  1:440333/100‑1 (MQ=255)
                                               cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                <  1:1934442/100‑1 (MQ=255)
                                                tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               <  1:31084/100‑1 (MQ=255)
                                                 gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               >  2:866472/1‑99 (MQ=255)
                                                 gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  2:2494404/1‑100 (MQ=255)
                                                 gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  2:903964/1‑100 (MQ=255)
                                                  cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                             >  1:2716201/1‑100 (MQ=255)
                                                  cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                             >  2:2311600/1‑100 (MQ=255)
                                                  cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                             <  2:2591471/100‑1 (MQ=255)
                                                      aaaTCAATTTATCTTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGACAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgcttgc                                         >  1:1903387/1‑100 (MQ=255)
                                                      aaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgcttgc                                         >  2:444436/1‑100 (MQ=255)
                                                         tCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGcttgcttg                                      <  1:2554370/100‑1 (MQ=255)
                                                          cAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCttgcttgc                                     <  1:536527/100‑1 (MQ=255)
                                                           aaTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCg                                    >  1:1851883/1‑100 (MQ=255)
                                                            aTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                   <  2:761438/100‑1 (MQ=255)
                                                            aTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                   >  1:972796/1‑100 (MQ=255)
                                                             tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  >  1:1707459/1‑100 (MQ=255)
                                                                 tatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAata                              <  2:358750/100‑1 (MQ=255)
                                                                  atTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                             <  2:1952836/100‑1 (MQ=255)
                                                                   ttAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATAtt                            >  1:2576595/1‑100 (MQ=255)
                                                                     taCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  1:1290117/99‑1 (MQ=255)
                                                                     aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  1:2637139/100‑1 (MQ=255)
                                                                     aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  2:1918433/100‑1 (MQ=255)
                                                                      aCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACATAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         <  2:1856340/100‑1 (MQ=255)
                                                                      aCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTga                         >  2:1290112/1‑99 (MQ=255)
                                                                       ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                        <  1:324207/100‑1 (MQ=255)
                                                                        cccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         <  1:1877902/98‑1 (MQ=255)
                                                                        cccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCaa                       >  1:1968621/1‑100 (MQ=255)
                                                                          cACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAata                     >  1:1081900/1‑100 (MQ=255)
                                                                          cACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAat                      <  1:279464/99‑1 (MQ=255)
                                                                            cTAATAGTGTATTCATATTATATTTTTTTT‑TTTTTTTTTAACTAATTATTAAATCTGTAATATTCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                   >  2:432894/1‑100 (MQ=255)
                                                                             tAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  2:526440/100‑1 (MQ=255)
                                                                             tAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  1:1055873/100‑1 (MQ=255)
                                                                             tAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  1:648227/100‑1 (MQ=255)
                                                                               aTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                >  2:2686597/1‑100 (MQ=255)
                                                                                tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:2282099/100‑1 (MQ=255)
                                                                                tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               >  2:613940/1‑100 (MQ=255)
                                                                                tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  1:1733498/100‑1 (MQ=255)
                                                                                tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:2735536/100‑1 (MQ=255)
                                                                                tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:293434/100‑1 (MQ=255)
                                                                                 aGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTTGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:1233852/99‑1 (MQ=255)
                                                                                  gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc             >  1:2262846/1‑100 (MQ=255)
                                                                                  gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc             >  1:1775474/1‑100 (MQ=255)
                                                                                  gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc             <  1:2691538/100‑1 (MQ=255)
                                                                                  gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc             <  1:908346/100‑1 (MQ=255)
                                                                                  gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc             <  1:521742/100‑1 (MQ=255)
                                                                                      aTTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaaaa         <  2:2132541/100‑1 (MQ=255)
                                                                                      aTTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaaaa         <  1:545633/100‑1 (MQ=255)
                                                                                         cATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGgt      >  1:2433273/1‑100 (MQ=255)
                                                                                            attatattTTTTTT‑ATATACTCAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt   >  2:2392102/1‑100 (MQ=255)
                                                                                            attata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt   >  1:2859324/1‑100 (MQ=255)
                                                                                             ttata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgttg  >  2:1322593/1‑100 (MQ=255)
                                                                                                  |                                                                                               
CAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTG  >  NC_000913/2469628‑2469821

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: ‑