Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,470,410 (T)7→6 100% coding (1280/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,470,4040T.97.0% 439.2 / 5.8 101coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base . (50/48);  total (51/50)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 114 total aligned reads displayed.

TTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTAT  >  NC_000913/2470308‑2470500
                                                                                                |                                                                                                
tttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttt                                                                                               <  1:1306527/100‑1 (MQ=255)
  tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                             <  1:377886/100‑1 (MQ=255)
  tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                             >  1:229847/1‑100 (MQ=255)
    tgtgTCATTTTTCGGGGAAGGATATAGATATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                           <  2:1677127/100‑3 (MQ=255)
     gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACATTCGttttttaat                                                                                          >  1:1545530/1‑97 (MQ=255)
     gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                          >  1:179403/1‑97 (MQ=255)
     gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                          >  2:1033194/1‑97 (MQ=255)
     gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                          >  1:179434/1‑97 (MQ=255)
      tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                         <  2:289111/100‑5 (MQ=255)
      tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                         >  1:1188166/1‑96 (MQ=255)
       gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                       >  2:1638517/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                     >  1:149551/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                     >  2:1733745/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                     <  1:1787976/100‑1 (MQ=255)
          aTTTTTGGGGGAAGTATATAGATATTTAAGCAGGCTTTTATTTGTAATGTATTTTTCGTTCGACTTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                     <  2:2134356/99‑1 (MQ=255)
           tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaa                                                                                   <  1:407330/100‑1 (MQ=255)
           tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaa                                                                                   <  2:737693/100‑1 (MQ=255)
             ttttGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                                 <  1:807604/100‑1 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAATCAATTATTTGGCATGTATTTTCCGCTCGACCTTTGCTTATACCTAACAGTCG‑TTTCTTAATTTATaaaa                                                                                 >  2:1476525/1‑100 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAATCAATTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTCATTTATaaaa                                                                                 >  2:1476514/1‑100 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATaaa                                                                                  <  1:106860/100‑1 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAAGCAGTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTGATTTATaaaa                                                                                 >  2:1476524/1‑100 (MQ=255)
              ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTCGAATGTATTTTTCGTTCGACCTTTGCTTATACATATCAGTCG‑TTTTTTAATTTATAAaat                                                                                <  2:226943/100‑1 (MQ=255)
                cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGTAATGTATTTTTCGTTCGACTTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaataa                                                                              <  2:1222428/100‑1 (MQ=255)
                 ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaataat                                                                             <  1:1563487/100‑1 (MQ=255)
                 gaggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaataat                                                                             <  2:1357730/98‑1 (MQ=255)
                    gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCa                                                                          <  2:1682007/100‑1 (MQ=255)
                      aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaaa                                                                        >  2:1177599/1‑100 (MQ=255)
                       ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                       <  2:2099123/100‑1 (MQ=255)
                       ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTATTCGTGCTACCTTTTCTTATACATAAAAATAG‑TTTTTTAATTTATAAAATAATTTAAAg                                                                       >  2:1554601/1‑100 (MQ=255)
                        gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGa                                                                      <  2:391096/100‑1 (MQ=255)
                        gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGa                                                                      >  1:1989012/1‑100 (MQ=255)
                          tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGaaa                                                                    <  1:2064257/100‑1 (MQ=255)
                           ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   <  1:2069875/100‑1 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGNAATGNATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAANGAAATc                                                                  >  1:1715818/1‑100 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTATTTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   >  2:755819/1‑99 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCTTTCGATCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                  >  2:24585/1‑100 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   >  1:1141681/1‑99 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                  <  1:329260/100‑1 (MQ=255)
                               agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  1:1995353/100‑1 (MQ=255)
                                gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                              >  1:1458541/1‑100 (MQ=255)
                                gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                              >  2:2002288/1‑100 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                             <  1:709329/100‑1 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                             <  2:1259120/100‑1 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGCAATGTATTTTTCGTTTGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  2:1302798/98‑1 (MQ=255)
                                     aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                         >  1:1763848/1‑100 (MQ=255)
                                       gcaagcaTTTATTTGGAATGTATTTTTCGTTCGAACTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                       <  2:1134270/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  1:1859852/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  2:452248/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  2:2072061/100‑1 (MQ=255)
                                           gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                   <  1:1370976/100‑1 (MQ=255)
                                           gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                   <  1:20754/100‑1 (MQ=255)
                                              tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                <  1:636643/100‑1 (MQ=255)
                                                tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                               >  1:15355/1‑99 (MQ=255)
                                                tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                              >  2:162029/1‑100 (MQ=255)
                                                 atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  2:1097433/100‑1 (MQ=255)
                                                 atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  1:2230015/100‑1 (MQ=255)
                                                  tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                            >  2:863525/1‑100 (MQ=255)
                                                   ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                            >  1:803906/1‑99 (MQ=255)
                                                   ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                           <  1:1977715/100‑1 (MQ=255)
                                                    tGGAATGTATTTTTCGTTCGACCTATGCTTATACATAACAGTCG‑TTTTTTACTTTCAAAAATCATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCa                                            >  2:1582565/1‑98 (MQ=255)
                                                        aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                      >  1:476404/1‑100 (MQ=255)
                                                          gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                    >  2:926432/1‑100 (MQ=255)
                                                            aTTTTTCGTTCGACCTTTGCTTATACATAACAGACG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAAACAATAGCGATGTAAAGCACTAAGTTTaa                                  >  2:1342636/1‑100 (MQ=255)
                                                              ttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAAt                                <  2:760707/100‑1 (MQ=255)
                                                               tttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                               >  1:2176514/1‑100 (MQ=255)
                                                               tttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                               >  1:1414896/1‑100 (MQ=255)
                                                                   gttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           >  2:1067550/1‑100 (MQ=255)
                                                                   gttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAAGTTAGAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           <  2:892792/100‑1 (MQ=255)
                                                                    ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           >  2:869847/1‑100 (MQ=255)
                                                                     tcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                         <  1:270730/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                     <  2:190870/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                     <  1:708674/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGACG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCAATGTAAAGCACTAATTTTAAATTTCGCGCCAAt                     >  2:880722/1‑100 (MQ=255)
                                                                           tttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                   >  2:341404/1‑100 (MQ=255)
                                                                              gCTTATACATAATAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                <  1:1395791/100‑1 (MQ=255)
                                                                              gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                <  1:2187935/100‑1 (MQ=255)
                                                                              gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTACGCGCCAATCATgg                >  2:2198875/1‑100 (MQ=255)
                                                                               cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc               >  2:2025296/1‑100 (MQ=255)
                                                                               cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGAGCCAATCATgg                >  2:606664/1‑99 (MQ=255)
                                                                                  ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg            <  1:267854/100‑1 (MQ=255)
                                                                                   taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc           >  1:1397103/1‑100 (MQ=255)
                                                                                   taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc           >  2:562369/1‑100 (MQ=255)
                                                                                    aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca          <  1:500011/100‑1 (MQ=255)
                                                                                    aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca          <  1:153847/100‑1 (MQ=255)
                                                                                     cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:556669/1‑100 (MQ=255)
                                                                                     cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:825374/1‑100 (MQ=255)
                                                                                     cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCAATGTAAAACACAACGTAAAAATTGCACGCCAAACATGgcg              >  2:540937/1‑95 (MQ=255)
                                                                                      aTAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                  >  2:1910009/1‑90 (MQ=255)
                                                                                       tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa       >  1:1633030/1‑100 (MQ=255)
                                                                                       tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa       >  2:2198844/1‑100 (MQ=255)
                                                                                         aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     >  2:495276/1‑100 (MQ=255)
                                                                                         aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     <  1:1738050/100‑1 (MQ=255)
                                                                                         aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     <  1:1991565/100‑1 (MQ=255)
                                                                                          cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaa    <  2:1139720/100‑7 (MQ=37)
                                                                                          cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    <  2:738436/100‑1 (MQ=255)
                                                                                          cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    >  2:1201114/1‑100 (MQ=255)
                                                                                          cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    >  1:39438/1‑100 (MQ=255)
                                                                                           aGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta   >  1:2031895/1‑100 (MQ=255)
                                                                                             tCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtat  <  1:649379/100‑1 (MQ=255)
                                                                                                |                                                                                                
TTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTAT  >  NC_000913/2470308‑2470500

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