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.6% 564.3 / 4.8 125coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (3/0);  new base . (62/60);  total (65/60)
Fisher's exact test for biased strand distribution p-value = 2.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 136 total aligned reads displayed.

ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470306‑2470498
                                                                                                  |                                                                                              
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                           <  2:1029448/100‑1 (MQ=255)
     aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  1:2406219/100‑2 (MQ=255)
     aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                          <  2:472409/100‑2 (MQ=255)
      tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaa                                                                                         <  1:2412675/100‑3 (MQ=255)
       gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                        <  2:1917076/100‑4 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                       <  2:1181241/100‑5 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     <  1:2713105/100‑1 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     >  2:1545219/1‑100 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAAttt                                                                                     <  2:1984710/100‑1 (MQ=255)
         gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGCTTTTTAAttt                                                                                     >  1:2268926/1‑100 (MQ=255)
            aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTTTTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                  >  1:3078029/1‑100 (MQ=255)
             tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtat                                                                                   <  1:2530320/98‑1 (MQ=255)
             tttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaa                                                                                 >  2:2739656/1‑100 (MQ=255)
              ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                <  2:862452/100‑1 (MQ=255)
              ttttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaa                                                                                >  1:2545133/1‑100 (MQ=255)
                  cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                            >  1:2617193/1‑100 (MQ=255)
                        aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTTTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                      >  1:1184215/1‑100 (MQ=255)
                        aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                      <  2:897943/100‑1 (MQ=255)
                        aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaaa                                                                      <  1:551871/100‑1 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                     <  2:2599846/100‑1 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTAATTTATAAAATAATTCAAAGa                                                                    >  1:1429049/1‑100 (MQ=255)
                            tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  <  2:306154/100‑1 (MQ=255)
                            tataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  <  2:3215182/100‑1 (MQ=255)
                             ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCTTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                 >  1:3194139/1‑100 (MQ=255)
                              taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaaa                                                                  <  1:348305/98‑1 (MQ=255)
                              taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                >  2:239539/1‑100 (MQ=255)
                                 agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                             <  2:377723/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                              >  1:2693515/1‑98 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            <  2:2545152/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                            >  2:467152/1‑100 (MQ=255)
                                   aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  2:2736219/100‑1 (MQ=255)
                                     ttAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                           <  1:2097381/98‑1 (MQ=255)
                                      tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCGAAGAAATCAAGACAAt                                                        >  1:2022997/1‑100 (MQ=255)
                                      tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                        <  2:2023014/100‑1 (MQ=255)
                                        agcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                      >  1:265216/1‑100 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                      >  1:2388942/1‑99 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                     <  2:2686230/100‑1 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTTTTTTTTTTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACACTatc                                                     >  1:850475/1‑97 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    >  2:62562/1‑100 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                    <  2:2301459/100‑1 (MQ=255)
                                             gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                 <  1:1726904/100‑1 (MQ=255)
                                              caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                >  2:3438152/1‑100 (MQ=255)
                                              caTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGt                                                >  2:2262472/1‑100 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               <  1:2722689/100‑1 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                               <  2:2070409/100‑1 (MQ=255)
                                                tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                              >  2:93380/1‑100 (MQ=255)
                                                 ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                             <  2:2961287/100‑1 (MQ=255)
                                                  tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                            <  2:618313/100‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           >  2:3456622/1‑100 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                           >  1:2753541/1‑100 (MQ=255)
                                                     ttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCAc                                         <  2:3306772/100‑1 (MQ=255)
                                                      tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                        >  2:3310275/1‑100 (MQ=255)
                                                       ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCagat                                       >  2:2485541/1‑97 (MQ=255)
                                                       ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       <  2:1711676/100‑1 (MQ=255)
                                                       ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                       >  2:724173/1‑100 (MQ=255)
                                                         aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                     <  2:2523088/100‑1 (MQ=255)
                                                           tGTATTTTTTTTTTTACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGtt                                   >  1:2852135/1‑100 (MQ=255)
                                                             tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 <  1:2260119/100‑1 (MQ=255)
                                                             tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 <  2:2334190/100‑1 (MQ=255)
                                                             tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                 <  2:1548102/100‑1 (MQ=255)
                                                              aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAGTAGCGATGTAAAGCACTAAGTTTa                                 >  1:2855810/1‑99 (MQ=255)
                                                              aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                >  2:726184/1‑100 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  2:2136464/100‑1 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  2:2221931/100‑1 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:2796102/100‑1 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  2:596734/100‑1 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              >  1:1495584/1‑100 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAt                              <  1:3071982/100‑1 (MQ=255)
                                                                   tcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                           <  2:2954338/100‑1 (MQ=255)
                                                                    cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                          >  2:763706/1‑100 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  1:2225222/100‑1 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTATAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                         <  2:1206313/100‑1 (MQ=255)
                                                                      ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        <  2:850482/100‑1 (MQ=255)
                                                                      ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        >  2:2850875/1‑100 (MQ=255)
                                                                      ttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                        <  1:88402/100‑1 (MQ=255)
                                                                         gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                     >  2:2884705/1‑100 (MQ=255)
                                                                         gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                     >  2:1775827/1‑100 (MQ=255)
                                                                         gACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCa                     <  1:2825155/100‑1 (MQ=255)
                                                                          aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                    >  2:1235963/1‑100 (MQ=255)
                                                                            cTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                  >  1:1870943/1‑100 (MQ=255)
                                                                             tttGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                 <  1:1013684/100‑1 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg              >  2:2093484/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgac             <  1:1288350/100‑3 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:2812389/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:2721742/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc             >  2:1104276/1‑100 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc         >  1:1138943/1‑100 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTTTTTTTTTTTATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATCGCGCGca        >  1:2556645/1‑100 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGtcgcgc         >  2:750480/1‑99 (MQ=255)
                                                                                      aCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca        <  2:1460416/100‑1 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGAACTAAGTTTAAATTAAGAGCCAATCATgaagcgaa        >  1:805882/1‑92 (MQ=255)
                                                                                        aTAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcaca      <  1:837409/100‑1 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa     <  1:1235952/100‑1 (MQ=255)
                                                                                          aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    <  2:1262279/100‑1 (MQ=255)
                                                                                          aaCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAg    <  1:2769863/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc   <  2:1765803/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc   <  2:1289171/100‑1 (MQ=255)
                                                                                           aCAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc   <  1:1345407/100‑1 (MQ=255)
                                                                                            cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATTAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  >  1:1848188/1‑100 (MQ=255)
                                                                                            cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATGAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt  >  1:1849460/1‑100 (MQ=255)
                                                                                                  |                                                                                              
ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCT  >  NC_000913/2470306‑2470498

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