Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,469,733 (T)8→7 100% 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.92.5% 300.2 / 18.4 80coding (596/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (6/0);  new base . (34/40);  total (40/40)
Fisher's exact test for biased strand distribution p-value = 2.55e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.10e-01

AATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTG  >  NC_000913/2469633‑2469821
                                                                                             |                                                                                               
aaTCTATGCTTTTGCTATTTTTTTCGTTTCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           <  1:2240873/100‑1 (MQ=255)
aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           >  1:806074/1‑100 (MQ=255)
aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           <  1:780883/100‑1 (MQ=255)
aaTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           <  1:339950/100‑1 (MQ=255)
 aTCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttt                                                                                           <  1:634603/99‑1 (MQ=255)
  tCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttga                                                                                         <  2:2432416/100‑3 (MQ=255)
   cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgat                                                                                        >  2:2519546/1‑97 (MQ=255)
   cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgat                                                                                        >  2:1420300/1‑97 (MQ=255)
   cTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgat                                                                                        <  1:1913700/100‑4 (MQ=255)
    tATGCTTTTGCTATTTTTTTTGTTGCTAATAATATCTACTGCTTAAAATCAATTTATATTAACTCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                       >  1:1305879/1‑96 (MQ=255)
    tATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                       <  1:638852/100‑5 (MQ=255)
    tATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAtttttttgata                                                                                       <  1:254615/100‑5 (MQ=255)
      tGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGa                                                                                    <  2:968635/100‑1 (MQ=255)
       gCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttttata                                                                                       >  1:2038358/1‑94 (MQ=255)
        cTTTTGCTACTTTTTTCTTTTCTAATAACTGCTATTGCTAAAAATCAATTTATTTTATTCCCACTAATGTTGTATTTATATTATAtttttttt                                                                                          <  1:1631536/93‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAAc                                                                             <  1:613921/100‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAAc                                                                             <  2:1326427/100‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAAc                                                                             <  2:525352/100‑1 (MQ=255)
              cTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACt                                                                            <  1:2185053/100‑1 (MQ=255)
               tATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTa                                                                           <  2:119002/100‑1 (MQ=255)
                aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTaa                                                                          >  1:188980/1‑100 (MQ=255)
                aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTaa                                                                          >  1:1921596/1‑100 (MQ=255)
                  ttttttCGTTGCTAATAATATCTATTGGTAAAAATCAATTTATATTAACCCCACTAATAGTATATTCATATTATA‑TTTTTTTGATAGACACAAACTAatt                                                                        <  1:2234913/100‑1 (MQ=255)
                   tttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAatta                                                                       <  2:2108456/100‑1 (MQ=255)
                    ttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑ATAGACACAAACTAattat                                                                      >  1:980918/1‑100 (MQ=255)
                    ttttCGTTGCTAATAATATCTACTGCTAAAAATCAATCTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAattat                                                                      <  1:1746893/100‑1 (MQ=255)
                     ttttgTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAattatt                                                                     <  2:914462/96‑1 (MQ=255)
                     tttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAattatt                                                                     <  1:499844/100‑1 (MQ=255)
                        cGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTaaa                                                                  <  2:1304517/100‑1 (MQ=255)
                        agTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGGGTAGTCATATTATA‑TGATTTTGATAGACACAAATTAatgatt                                                                     >  1:213693/2‑97 (MQ=255)
                            gCTAATAATATCTACTGCTAAAAATCAATTTATATTAGCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCt                                                               >  1:1389125/1‑99 (MQ=255)
                              taataaTATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTTTTA‑AAACTAACTAATTATTAAATCTGt                                                             >  1:819404/1‑99 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           >  1:1214878/1‑100 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           >  2:1521236/1‑100 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACGAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           <  1:2314776/100‑1 (MQ=255)
                                ataataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACATAAACTAATTATTAAATCTGTaa                                                           <  1:1500671/99‑1 (MQ=255)
                                 taataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAata                                                         <  2:450055/100‑1 (MQ=255)
                                       cTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   <  2:1583667/100‑1 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   >  1:258823/1‑100 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑ATATATACAAACTAATTATTAAATCTGTAATATGCgtg                                                   >  1:1892257/1‑99 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtgg                                                  <  2:2015273/100‑1 (MQ=255)
                                          cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                <  1:1088714/100‑1 (MQ=255)
                                           tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               >  2:2219125/1‑100 (MQ=255)
                                           tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               <  1:2219110/100‑1 (MQ=255)
                                            gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  1:373346/1‑100 (MQ=255)
                                            gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  2:827842/1‑100 (MQ=255)
                                             cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                             <  2:694104/100‑1 (MQ=255)
                                              tAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgct                                            >  2:2493837/1‑100 (MQ=255)
                                               aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                           <  1:2542554/100‑1 (MQ=255)
                                                  aaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgcttgct                                        <  2:1521422/100‑1 (MQ=255)
                                                    tCAATTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGcttgcttg                                      <  1:2364225/100‑1 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑TTTTTTTCAAACTAATTATTAAATTTGTAATATGCGTGGTGTGCTTGCTTGCg                                    >  1:633712/1‑100 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATATACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTgcttgctt                                       >  1:1316279/1‑97 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGAGTATTCATATTATA‑TTTTTTTGATATACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtactttct                                        >  1:2312138/1‑96 (MQ=255)
                                                       aTTTATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                   <  1:975423/100‑1 (MQ=255)
                                                        tttATATTAACTCTACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTAGTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  1:643568/100‑1 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGTTTGCGTc                                  >  2:1549731/1‑100 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  2:272889/100‑1 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  2:1208421/100‑1 (MQ=255)
                                                         ttATATTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  1:602219/99‑1 (MQ=255)
                                                          tataTTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCaa                                >  2:2212498/1‑100 (MQ=255)
                                                           atatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAat                               >  2:1395098/1‑100 (MQ=255)
                                                            tatTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAACCTAATTCTTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAtta                              >  2:40563/1‑100 (MQ=255)
                                                            tatTAACCCCACTAATAGTGTATTCATATTATA‑TTTGTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAata                              >  2:1822405/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATATTTTTTTT‑TTTTTTTTAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTAAatat                             >  1:2074518/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                             >  2:381792/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                             >  1:909680/1‑100 (MQ=255)
                                                              ttAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATAtt                            >  1:1046105/1‑100 (MQ=255)
                                                               tAACCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTTCTTGCGTCAATAttt                           >  1:1667346/1‑100 (MQ=255)
                                                                aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  2:2413012/100‑1 (MQ=255)
                                                                aaCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          >  2:1088196/1‑100 (MQ=255)
                                                                 aCCCCACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         >  1:984790/1‑100 (MQ=255)
                                                                  ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         <  1:110046/99‑1 (MQ=255)
                                                                  ccccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                        >  2:682246/1‑100 (MQ=255)
                                                                   cccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCaa                       >  2:1715094/1‑100 (MQ=255)
                                                                    ccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGATGTGCTTGCTTGCGTCAATATTTGCAat                      >  1:649930/1‑100 (MQ=255)
                                                                    ccACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTCGTGTGCTTGCTTGCGTCAATATTTg                          >  1:2331301/1‑96 (MQ=255)
                                                                     cACTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAata                     >  2:1635652/1‑100 (MQ=255)
                                                                       cTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                   >  1:1535185/1‑100 (MQ=255)
                                                                       cTAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                   <  1:20768/100‑1 (MQ=255)
                                                                        tAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  2:622894/100‑1 (MQ=255)
                                                                        tAATAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  1:1519559/100‑1 (MQ=255)
                                                                         aaTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  1:2050717/91‑1 (MQ=255)
                                                                          aTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                >  2:1469690/1‑100 (MQ=255)
                                                                          aTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                <  1:1040983/100‑1 (MQ=255)
                                                                          aTAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                >  1:1323442/1‑100 (MQ=255)
                                                                           tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:188981/100‑1 (MQ=255)
                                                                           tAGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:2438807/100‑1 (MQ=255)
                                                                            aGTGTATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAtt              >  2:1556644/1‑100 (MQ=255)
                                                                              tgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCa            <  2:633715/100‑1 (MQ=255)
                                                                              tgtATTCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCa            <  2:1892267/100‑1 (MQ=255)
                                                                                 aTTCATATTATA‑TTTTTTTGATAGACACAAACTAATTGTTAAATTTGTAATATTCGTGGTGTGCTTGCTTGAGACAATATTTGCAATATTTTTTTaaaa          >  1:1770162/1‑99 (MQ=255)
                                                                                   tCATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAgg       <  2:2141707/100‑1 (MQ=255)
                                                                                    cATATTATA‑TTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGgt      >  1:908516/1‑100 (MQ=255)
                                                                                        ttata‑ttTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgttg  >  1:1015228/1‑100 (MQ=255)
                                                                                             |                                                                                               
AATCTATGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGTGTTG  >  NC_000913/2469633‑2469821

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