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)
       gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                       >  2:1638517/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                     >  1:149551/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                     >  2:1733745/1‑100 (MQ=255)
         cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                     <  1:1787976/100‑1 (MQ=255)
          aTTTTTGGGGGAAGTATATAGATATTTAAGCAGGCTTTTATTTGTAATGTATTTTTCGTTCGACTTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                     <  2:2134356/99‑1 (MQ=255)
           tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                   <  1:407330/100‑1 (MQ=255)
           tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                   <  2:737693/100‑1 (MQ=255)
             ttttGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaaa                                                                                 <  1:807604/100‑1 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAATCAATTATTTGGCATGTATTTTCCGCTCGACCTTTGCTTATACCTAACAGTCGTTTCTTAATTTATaaaa                                                                                 >  2:1476525/1‑100 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAATCAATTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTCATTTATaaaa                                                                                 >  2:1476514/1‑100 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATaaa                                                                                  <  1:106860/100‑1 (MQ=255)
             tttCGGGGAAGGATATAGAGATTTAAGCAAGCAGTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTGATTTATaaaa                                                                                 >  2:1476524/1‑100 (MQ=255)
              ttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTCGAATGTATTTTTCGTTCGACCTTTGCTTATACATATCAGTCGTTTTTTAATTTATAAaat                                                                                <  2:226943/100‑1 (MQ=255)
                cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGTAATGTATTTTTCGTTCGACTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                              <  2:1222428/100‑1 (MQ=255)
                 ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                             <  1:1563487/100‑1 (MQ=255)
                 gaggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                             <  2:1357730/98‑1 (MQ=255)
                    gAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCa                                                                          <  2:1682007/100‑1 (MQ=255)
                      aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                        >  2:1177599/1‑100 (MQ=255)
                       ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                       <  2:2099123/100‑1 (MQ=255)
                       ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTATTCGTGCTACCTTTTCTTATACATAAAAATAGTTTTTTAATTTATAAAATAATTTAAAg                                                                       >  2:1554601/1‑100 (MQ=255)
                        gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                      <  2:391096/100‑1 (MQ=255)
                        gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                      >  1:1989012/1‑100 (MQ=255)
                          tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                    <  1:2064257/100‑1 (MQ=255)
                           ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   <  1:2069875/100‑1 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGNAATGNATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAANGAAATc                                                                  >  1:1715818/1‑100 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTATTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   >  2:755819/1‑99 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCTTTCGATCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                  >  2:24585/1‑100 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                   >  1:1141681/1‑99 (MQ=255)
                            taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                  <  1:329260/100‑1 (MQ=255)
                               agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  1:1995353/100‑1 (MQ=255)
                                gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                              >  1:1458541/1‑100 (MQ=255)
                                gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                              >  2:2002288/1‑100 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                             <  1:709329/100‑1 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                             <  2:1259120/100‑1 (MQ=255)
                                 aTTTAAGCAAGCATTTATTTGCAATGTATTTTTCGTTTGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                               <  2:1302798/98‑1 (MQ=255)
                                     aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                         >  1:1763848/1‑100 (MQ=255)
                                       gcaagcaTTTATTTGGAATGTATTTTTCGTTCGAACTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                       <  2:1134270/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  1:1859852/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  2:452248/100‑1 (MQ=255)
                                         aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                     <  2:2072061/100‑1 (MQ=255)
                                           gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                   <  1:1370976/100‑1 (MQ=255)
                                           gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                   <  1:20754/100‑1 (MQ=255)
                                              tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                <  1:636643/100‑1 (MQ=255)
                                                tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                               >  1:15355/1‑99 (MQ=255)
                                                tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                              >  2:162029/1‑100 (MQ=255)
                                                 atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  2:1097433/100‑1 (MQ=255)
                                                 atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                             <  1:2230015/100‑1 (MQ=255)
                                                  tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                            >  2:863525/1‑100 (MQ=255)
                                                   ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                            >  1:803906/1‑99 (MQ=255)
                                                   ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                           <  1:1977715/100‑1 (MQ=255)
                                                    tGGAATGTATTTTTCGTTCGACCTATGCTTATACATAACAGTCGTTTTTTACTTTCAAAAATCATTCAAAGAAATCAAGACAAAAGCGATGTAAAGCa                                            >  2:1582565/1‑98 (MQ=255)
                                                        aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                      >  1:476404/1‑100 (MQ=255)
                                                          gTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGttt                                    >  2:926432/1‑100 (MQ=255)
                                                            aTTTTTCGTTCGACCTTTGCTTATACATAACAGACGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAAACAATAGCGATGTAAAGCACTAAGTTTaa                                  >  2:1342636/1‑100 (MQ=255)
                                                              ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCTATGTAAAGCACTAAGTTTAAAt                                <  2:760707/100‑1 (MQ=255)
                                                               tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                               >  1:2176514/1‑100 (MQ=255)
                                                               tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                               >  1:1414896/1‑100 (MQ=255)
                                                                   gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           >  2:1067550/1‑100 (MQ=255)
                                                                   gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAAGTTAGAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           <  2:892792/100‑1 (MQ=255)
                                                                    ttcgACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                           >  2:869847/1‑100 (MQ=255)
                                                                     tcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                         <  1:270730/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                     <  2:190870/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAAt                     <  1:708674/100‑1 (MQ=255)
                                                                         ccTTTGCTTATACATAACAGACGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCAATGTAAAGCACTAATTTTAAATTTCGCGCCAAt                     >  2:880722/1‑100 (MQ=255)
                                                                           tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                   >  2:341404/1‑100 (MQ=255)
                                                                              gCTTATACATAATAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                <  1:1395791/100‑1 (MQ=255)
                                                                              gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg                <  1:2187935/100‑1 (MQ=255)
                                                                              gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTACGCGCCAATCATgg                >  2:2198875/1‑100 (MQ=255)
                                                                               cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc               >  2:2025296/1‑100 (MQ=255)
                                                                               cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGAGCCAATCATgg                >  2:606664/1‑99 (MQ=255)
                                                                                  ataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg            <  1:267854/100‑1 (MQ=255)
                                                                                   taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc           >  1:1397103/1‑100 (MQ=255)
                                                                                   taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc           >  2:562369/1‑100 (MQ=255)
                                                                                    aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca          <  1:500011/100‑1 (MQ=255)
                                                                                    aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca          <  1:153847/100‑1 (MQ=255)
                                                                                     cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:556669/1‑100 (MQ=255)
                                                                                     cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac         >  1:825374/1‑100 (MQ=255)
                                                                                     cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCAATGTAAAACACAACGTAAAAATTGCACGCCAAACATGgcg              >  2:540937/1‑95 (MQ=255)
                                                                                      aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCAt                  >  2:1910009/1‑90 (MQ=255)
                                                                                       tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa       >  1:1633030/1‑100 (MQ=255)
                                                                                       tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa       >  2:2198844/1‑100 (MQ=255)
                                                                                         aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATATCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     >  2:495276/1‑100 (MQ=255)
                                                                                         aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     <  1:1738050/100‑1 (MQ=255)
                                                                                         aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc     <  1:1991565/100‑1 (MQ=255)
                                                                                          cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatttaa    <  2:1139720/100‑7 (MQ=37)
                                                                                          cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    <  2:738436/100‑1 (MQ=255)
                                                                                          cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    >  2:1201114/1‑100 (MQ=255)
                                                                                          cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt    >  1:39438/1‑100 (MQ=255)
                                                                                           aGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCta   >  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: