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)
      tGCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGa                                                                                    <  2:968635/100‑1 (MQ=255)
       gCTTTTGCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATAttttttttata                                                                                       >  1:2038358/1‑94 (MQ=255)
        cTTTTGCTACTTTTTTCTTTTCTAATAACTGCTATTGCTAAAAATCAATTTATTTTATTCCCACTAATGTTGTATTTATATTATAtttttttt                                                                                          <  1:1631536/93‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAAc                                                                             <  1:613921/100‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAAc                                                                             <  2:1326427/100‑1 (MQ=255)
             gCTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAAc                                                                             <  2:525352/100‑1 (MQ=255)
              cTATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACt                                                                            <  1:2185053/100‑1 (MQ=255)
               tATTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTa                                                                           <  2:119002/100‑1 (MQ=255)
                aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTaa                                                                          >  1:188980/1‑100 (MQ=255)
                aTTTTTTTCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTaa                                                                          >  1:1921596/1‑100 (MQ=255)
                  ttttttCGTTGCTAATAATATCTATTGGTAAAAATCAATTTATATTAACCCCACTAATAGTATATTCATATTATATTTTTTTGATAGACACAAACTAatt                                                                        <  1:2234913/100‑1 (MQ=255)
                   tttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAatta                                                                       <  2:2108456/100‑1 (MQ=255)
                    ttttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTATAGACACAAACTAattat                                                                      >  1:980918/1‑100 (MQ=255)
                    ttttCGTTGCTAATAATATCTACTGCTAAAAATCAATCTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAattat                                                                      <  1:1746893/100‑1 (MQ=255)
                     ttttgTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAattatt                                                                     <  2:914462/96‑1 (MQ=255)
                     tttCGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAattatt                                                                     <  1:499844/100‑1 (MQ=255)
                        cGTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTaaa                                                                  <  2:1304517/100‑1 (MQ=255)
                        agTTGCTAATAATATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGGGTAGTCATATTATATGATTTTGATAGACACAAATTAatgatt                                                                     >  1:213693/2‑97 (MQ=255)
                            gCTAATAATATCTACTGCTAAAAATCAATTTATATTAGCCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCt                                                               >  1:1389125/1‑99 (MQ=255)
                              taataaTATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTTTTAAAACTAACTAATTATTAAATCTGt                                                             >  1:819404/1‑99 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           >  1:1214878/1‑100 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           >  2:1521236/1‑100 (MQ=255)
                               aataatATCTACTGCTAAAAATCAATTTATATTAACCCCACGAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTaa                                                           <  1:2314776/100‑1 (MQ=255)
                                ataataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACATAAACTAATTATTAAATCTGTaa                                                           <  1:1500671/99‑1 (MQ=255)
                                 taataTCTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAata                                                         <  2:450055/100‑1 (MQ=255)
                                       cTACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   <  2:1583667/100‑1 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtg                                                   >  1:258823/1‑100 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTATATATACAAACTAATTATTAAATCTGTAATATGCgtg                                                   >  1:1892257/1‑99 (MQ=255)
                                        tACTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtgg                                                  <  2:2015273/100‑1 (MQ=255)
                                          cTGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCgtggtg                                                <  1:1088714/100‑1 (MQ=255)
                                           tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               >  2:2219125/1‑100 (MQ=255)
                                           tGCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGgtgt                                               <  1:2219110/100‑1 (MQ=255)
                                            gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  1:373346/1‑100 (MQ=255)
                                            gCTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGtgtg                                              >  2:827842/1‑100 (MQ=255)
                                             cTAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgc                                             <  2:694104/100‑1 (MQ=255)
                                              tAAAAATCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgct                                            >  2:2493837/1‑100 (MQ=255)
                                               aaaaaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgctt                                           <  1:2542554/100‑1 (MQ=255)
                                                  aaTCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtgcttgct                                        <  2:1521422/100‑1 (MQ=255)
                                                    tCAATTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGcttgcttg                                      <  1:2364225/100‑1 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTTTTTTTTCAAACTAATTATTAAATTTGTAATATGCGTGGTGTGCTTGCTTGCg                                    >  1:633712/1‑100 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATATACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTgcttgctt                                       >  1:1316279/1‑97 (MQ=255)
                                                      aaTTTATATTAACCCCACTAATAGAGTATTCATATTATATTTTTTTGATATACACAAACTAATTATTAAATCTGTAATATGCGTGGTGtactttct                                        >  1:2312138/1‑96 (MQ=255)
                                                       aTTTATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGt                                   <  1:975423/100‑1 (MQ=255)
                                                        tttATATTAACTCTACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTAGTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  1:643568/100‑1 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGTTTGCGTc                                  >  2:1549731/1‑100 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  2:272889/100‑1 (MQ=255)
                                                        tttATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  2:1208421/100‑1 (MQ=255)
                                                         ttATATTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTc                                  <  1:602219/99‑1 (MQ=255)
                                                          tataTTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCaa                                >  2:2212498/1‑100 (MQ=255)
                                                           atatTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAat                               >  2:1395098/1‑100 (MQ=255)
                                                            tatTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAACCTAATTCTTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAtta                              >  2:40563/1‑100 (MQ=255)
                                                            tatTAACCCCACTAATAGTGTATTCATATTATATTTGTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAata                              >  2:1822405/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATATTTTTTTTTTTTTTTTAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTAAatat                             >  1:2074518/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                             >  2:381792/1‑100 (MQ=255)
                                                             atTAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAatat                             >  1:909680/1‑100 (MQ=255)
                                                              ttAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATAtt                            >  1:1046105/1‑100 (MQ=255)
                                                               tAACCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTTCTTGCGTCAATAttt                           >  1:1667346/1‑100 (MQ=255)
                                                                aaCCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  2:2413012/100‑1 (MQ=255)
                                                                aaCCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          >  2:1088196/1‑100 (MQ=255)
                                                                 aCCCCACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         >  1:984790/1‑100 (MQ=255)
                                                                  ccccACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGc                         <  1:110046/99‑1 (MQ=255)
                                                                  ccccACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCa                        >  2:682246/1‑100 (MQ=255)
                                                                   cccACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCaa                       >  2:1715094/1‑100 (MQ=255)
                                                                    ccACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGATGTGCTTGCTTGCGTCAATATTTGCAat                      >  1:649930/1‑100 (MQ=255)
                                                                    ccACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTCGTGTGCTTGCTTGCGTCAATATTTg                          >  1:2331301/1‑96 (MQ=255)
                                                                     cACTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAata                     >  2:1635652/1‑100 (MQ=255)
                                                                       cTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                   >  1:1535185/1‑100 (MQ=255)
                                                                       cTAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATc                   <  1:20768/100‑1 (MQ=255)
                                                                        tAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  2:622894/100‑1 (MQ=255)
                                                                        tAATAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCt                  <  1:1519559/100‑1 (MQ=255)
                                                                         aaTAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTg                          <  1:2050717/91‑1 (MQ=255)
                                                                          aTAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                >  2:1469690/1‑100 (MQ=255)
                                                                          aTAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                <  1:1040983/100‑1 (MQ=255)
                                                                          aTAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTa                >  1:1323442/1‑100 (MQ=255)
                                                                           tAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:188981/100‑1 (MQ=255)
                                                                           tAGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAt               <  2:2438807/100‑1 (MQ=255)
                                                                            aGTGTATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTAtt              >  2:1556644/1‑100 (MQ=255)
                                                                              tgtATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCa            <  2:633715/100‑1 (MQ=255)
                                                                              tgtATTCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCa            <  2:1892267/100‑1 (MQ=255)
                                                                                 aTTCATATTATATTTTTTTGATAGACACAAACTAATTGTTAAATTTGTAATATTCGTGGTGTGCTTGCTTGAGACAATATTTGCAATATTTTTTTaaaa          >  1:1770162/1‑99 (MQ=255)
                                                                                   tCATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAgg       <  2:2141707/100‑1 (MQ=255)
                                                                                    cATATTATATTTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGgt      >  1:908516/1‑100 (MQ=255)
                                                                                        ttatattTTTTTGATAGACACAAACTAATTATTAAATCTGTAATATGCGTGGTGTGCTTGCTTGCGTCAATATTTGCAATATCTTATTCAAAAGGtgttg  >  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: