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.3% 588.4 / 14.6 135coding (596/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (5/0);  new base . (72/58);  total (77/58)
Fisher's exact test for biased strand distribution p-value = 7.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.55e-01

Only 100 of 167 total aligned reads displayed.

ACAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTGT  >  NC_000913/2469627‑2469822
                                                                                                   |                                                                                                
aCAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCatattatat                                                                                                  >  2:1449795/1‑100 (MQ=255)
 cAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCatattatatt                                                                                                 >  1:2514359/1‑100 (MQ=255)
  aaTAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttt                                                                                                >  1:1924068/1‑100 (MQ=255)
  aaTAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttt                                                                                                >  1:122404/1‑100 (MQ=255)
   aTAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttt                                                                                               >  2:3076016/1‑100 (MQ=255)
      attCTATTCTTTTGCTATTTTTTTCTTTGCTAATAATATCTACTGCTAAAAATCAATTTATTTTATCCCCACTAATAGTGTATTCATATTATAttttttt                                                                                            <  2:2441364/98‑1 (MQ=255)
      aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATGTTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                            <  2:2101125/100‑1 (MQ=255)
      aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                            <  1:1078073/100‑1 (MQ=255)
      aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                            >  1:453324/1‑100 (MQ=255)
       aTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttg                                                                                           <  1:2589842/100‑2 (MQ=255)
         cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgat                                                                                         <  2:2514351/100‑4 (MQ=255)
         cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgat                                                                                         >  1:2440082/1‑97 (MQ=255)
          tctGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                        <  2:2067060/98‑5 (MQ=255)
          tATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                        <  2:1248391/100‑5 (MQ=255)
          tATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                        <  2:117395/100‑5 (MQ=255)
           aTGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAg                                                                                      <  2:692192/100‑1 (MQ=255)
           aTGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAg                                                                                      <  1:310013/100‑1 (MQ=255)
             gCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGac                                                                                    >  1:133528/1‑100 (MQ=255)
              cTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGaca                                                                                   >  1:813975/1‑100 (MQ=255)
                       tttttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAat                                                                          <  1:307413/100‑1 (MQ=255)
                       tttttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAat                                                                          <  1:3088520/100‑1 (MQ=255)
                          ttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAattat                                                                       <  1:1609931/100‑1 (MQ=255)
                            ttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAAttatta                                                                     <  1:216657/100‑1 (MQ=255)
                              cGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTaaa                                                                   >  2:3216931/1‑100 (MQ=255)
                               gTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑TTATTCTTAAACTAATTATTAAAt                                                                  >  2:3168756/1‑100 (MQ=255)
                                  gCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTg                                                               <  2:3336775/100‑1 (MQ=255)
                                   cTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑ATAGACACAAACTAATTATTAAATTTGt                                                              >  2:1082855/1‑100 (MQ=255)
                                     aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                            >  1:3181698/1‑100 (MQ=255)
                                      ataataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAat                                                           >  1:282388/1‑100 (MQ=255)
                                      ataataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAat                                                           <  1:61426/100‑1 (MQ=255)
                                       taataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAata                                                          <  2:284501/100‑1 (MQ=255)
                                        aataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTTATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAatat                                                         <  2:765665/100‑1 (MQ=255)
                                        aataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAatat                                                         >  2:567079/1‑100 (MQ=255)
                                         atatCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATg                                                        <  1:1506905/100‑1 (MQ=255)
                                          tatCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATATACACAAACTAATTATTAAATCTGTAATATGc                                                       >  2:3050795/1‑100 (MQ=255)
                                            tCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATgaga                                                     >  2:2419516/1‑97 (MQ=255)
                                             cTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑TTTTTTTTAAACTAATTATTAAATCTGTAATATGCgtg                                                    >  2:1679715/1‑100 (MQ=255)
                                              tACTGCTCAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtgg                                                   >  1:1796930/1‑100 (MQ=255)
                                                cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                 >  1:1615447/1‑100 (MQ=255)
                                                cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                 >  2:592325/1‑100 (MQ=255)
                                                 tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                                <  2:1798197/100‑1 (MQ=255)
                                                 tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                                <  2:1498443/100‑1 (MQ=255)
                                                 tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                                <  1:2502563/100‑1 (MQ=255)
                                                  gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                               >  1:1631197/1‑100 (MQ=255)
                                                  gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                               >  2:944128/1‑100 (MQ=255)
                                                  gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                               >  1:401116/1‑100 (MQ=255)
                                                  gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                               >  1:3182111/1‑100 (MQ=255)
                                                  gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                               >  2:2945576/1‑100 (MQ=255)
                                                   cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                              >  1:717745/1‑100 (MQ=255)
                                                    tAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgct                                             >  1:320803/1‑100 (MQ=255)
                                                    tAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgct                                             >  1:113094/1‑100 (MQ=255)
                                                     aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑TTTTTTTCAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                            >  2:488824/1‑100 (MQ=255)
                                                     aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                            >  1:633541/1‑100 (MQ=255)
                                                     aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                            >  1:1192329/1‑100 (MQ=255)
                                                     aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                            <  1:3367406/100‑1 (MQ=255)
                                                        aaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgcttgct                                         <  2:320802/100‑1 (MQ=255)
                                                            aaTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCg                                     <  1:1893972/100‑1 (MQ=255)
                                                             aTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                    <  2:2508212/100‑1 (MQ=255)
                                                             aTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                    >  1:852405/1‑100 (MQ=255)
                                                              tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                   <  2:1982219/100‑1 (MQ=255)
                                                              tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                   >  1:678553/1‑100 (MQ=255)
                                                              tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                   <  1:782542/100‑1 (MQ=255)
                                                                 atatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAat                                <  2:404370/100‑1 (MQ=255)
                                                                  tatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAGATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAata                               >  2:2288011/1‑100 (MQ=255)
                                                                  tatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAata                               <  1:29484/100‑1 (MQ=255)
                                                                   atTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                              >  2:1471513/1‑100 (MQ=255)
                                                                     tAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATAttt                            >  1:2684573/1‑100 (MQ=255)
                                                                     tAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATAttt                            <  1:814760/100‑1 (MQ=255)
                                                                      aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                           >  1:1377721/1‑100 (MQ=255)
                                                                      aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                           <  2:1093799/100‑1 (MQ=255)
                                                                       aCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                          <  2:1515941/100‑1 (MQ=255)
                                                                        ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                         >  1:2006571/1‑100 (MQ=255)
                                                                        ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                         <  1:2189974/100‑1 (MQ=255)
                                                                        ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                         >  1:2006472/1‑100 (MQ=255)
                                                                        ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                         <  2:79718/100‑1 (MQ=255)
                                                                         cccTCTACTTGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                         <  1:2838507/99‑1 (MQ=255)
                                                                          ccACTAATAGTGTACTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAat                       >  1:31371/1‑100 (MQ=255)
                                                                             cTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGTGTTAATATTTGCAatat                     >  2:779954/1‑99 (MQ=255)
                                                                             cTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATGATTAAATCTGTAATATTCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                    >  2:2841627/1‑100 (MQ=255)
                                                                               aaTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCtn                  <  2:1160308/100‑2 (MQ=255)
                                                                                 tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt                <  2:2154107/100‑1 (MQ=255)
                                                                                 tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt                <  1:1623006/100‑1 (MQ=255)
                                                                                   gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc              <  2:2684564/100‑1 (MQ=255)
                                                                                   gtgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTc              <  2:401115/100‑1 (MQ=255)
                                                                                    tgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCa             >  1:2653099/1‑100 (MQ=255)
                                                                                     gtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaa            >  2:665473/1‑100 (MQ=255)
                                                                                      tATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaaa           <  2:739420/100‑1 (MQ=255)
                                                                                      tATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaaa           >  1:2093271/1‑100 (MQ=255)
                                                                                        ttCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAg         >  2:2475652/1‑100 (MQ=255)
                                                                                        ttCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTAGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCaaa           >  1:2949629/1‑98 (MQ=255)
                                                                                        ctCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAg         <  2:1047305/99‑1 (MQ=255)
                                                                                         tCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAgg        >  2:98177/1‑100 (MQ=255)
                                                                                          cATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAgag       >  1:1047306/1‑98 (MQ=255)
                                                                                          cATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGgt       <  1:2932385/100‑1 (MQ=255)
                                                                                            tattata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGgtg      >  2:1309652/1‑99 (MQ=255)
                                                                                             attata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt    >  1:2316258/1‑100 (MQ=255)
                                                                                             attata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt    >  1:77649/1‑100 (MQ=255)
                                                                                             attata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt    >  2:1703277/1‑100 (MQ=255)
                                                                                             attata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgtt    >  2:1500590/1‑100 (MQ=255)
                                                                                              ttata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgttg   <  2:1937992/100‑1 (MQ=255)
                                                                                                   |                                                                                                
ACAATAAATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTGT  >  NC_000913/2469627‑2469822

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