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.98.6% 640.1 / 1.6 138coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base . (71/65);  total (72/66)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

Only 100 of 154 total aligned reads displayed.

TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470307‑2470499
                                                                                                 |                                                                                               
ttttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttt                                                                                               <  1:2278033/100‑1 (MQ=255)
  ttATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttt                                                                                             <  1:640368/100‑1 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                           <  2:662730/100‑2 (MQ=255)
    aTGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttta                                                                                           <  1:2143961/100‑2 (MQ=255)
     tgtgTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTa                                                                                          <  1:3129604/100‑1 (MQ=255)
      gtgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         >  1:2932633/1‑97 (MQ=255)
       tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                        <  2:165237/100‑5 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaat                                                                                         <  2:3218332/98‑4 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                      <  2:2302978/100‑1 (MQ=255)
        gtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAAttt                                                                                      <  1:221050/100‑1 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTta                                                                                     <  2:1380368/100‑1 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTta                                                                                     >  2:1066861/1‑100 (MQ=255)
         tCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTta                                                                                     <  1:1912343/100‑1 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                    >  1:2810987/1‑100 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                    >  1:2375236/1‑100 (MQ=255)
          cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                    >  1:1936934/1‑100 (MQ=255)
           aTTTTTCTGGTAAGGATATAGATATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTTGACCTTTGCTTATACATAACAGTCGttttttatttttt                                                                                     >  2:416301/1‑92 (MQ=255)
           aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtata                                                                                   <  2:84560/100‑1 (MQ=255)
           aTTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTtat                                                                                    >  2:2928274/1‑99 (MQ=255)
              tttCGGGGAAGGATATAGAGATTTAAGCATGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                                <  2:1968399/100‑1 (MQ=255)
              tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                                >  1:927586/1‑100 (MQ=255)
              tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATaaaa                                                                                >  1:2795842/1‑100 (MQ=255)
               ttCGGTGAAAGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAaat                                                                               <  2:2551138/100‑1 (MQ=255)
                   gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                           <  1:1841817/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTTTTCTTTTGCTTATACATAACCGTCG‑TTTTTTCATTTATCAAATAATTc                                                                          >  2:1891711/1‑100 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAAtt                                                                           <  2:1704149/99‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  1:1028234/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  1:1240703/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  1:373860/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  1:2980622/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          <  1:2961990/100‑1 (MQ=255)
                    ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTc                                                                          >  1:1689014/1‑100 (MQ=255)
                       aGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCaaa                                                                       >  1:1713083/1‑100 (MQ=255)
                        ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                      <  1:1822783/100‑1 (MQ=255)
                        ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAg                                                                      >  1:3081235/1‑100 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 >  2:2993185/1‑100 (MQ=255)
                             taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                 >  2:2716454/1‑100 (MQ=255)
                              agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                >  2:728477/1‑100 (MQ=255)
                                agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                               <  2:1881698/99‑1 (MQ=255)
                                agaTTTAAGCAAGCATTTATTTGAAATGTATTTTTCGTTAGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                              <  2:2797497/100‑1 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             >  1:659558/1‑100 (MQ=255)
                                 gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                             <  2:631725/100‑1 (MQ=255)
                                   tttgcgcaagcaTTTATTTGGCATGCATTTTTCGTTCGAACTTTGCATATACGTACCAGTCG‑TTGTTTAATTTATACAATAATTAAAAGAAAt                                                                  >  2:3236658/6‑93 (MQ=255)
                                     tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                         >  2:2202155/1‑100 (MQ=255)
                                      aagcaagcATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                        >  1:548589/1‑100 (MQ=255)
                                        gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAg                                                       <  2:1518716/99‑1 (MQ=255)
                                        gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                      >  1:177064/1‑100 (MQ=255)
                                         caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                      <  2:2933108/99‑1 (MQ=255)
                                         caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                     >  2:3152036/1‑100 (MQ=255)
                                          aagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                    <  2:2365862/100‑1 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                               >  1:3129115/1‑100 (MQ=255)
                                               tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                               >  1:548035/1‑100 (MQ=255)
                                                ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                              >  2:2798776/1‑100 (MQ=255)
                                                 tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                             <  1:1364726/100‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            >  2:308831/1‑100 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGCCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  1:912301/100‑1 (MQ=255)
                                                  atttGGAATGTATTTTTCGTTCGACATTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                            <  2:2398599/100‑1 (MQ=255)
                                                   tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                           <  2:537124/100‑1 (MQ=255)
                                                     tGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACt                                         >  2:2549312/1‑100 (MQ=255)
                                                       gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                       >  1:1481826/1‑100 (MQ=255)
                                                       gAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTaa                                       >  1:3228275/1‑100 (MQ=255)
                                                        aaTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAg                                      >  1:2330329/1‑100 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                  >  2:1587193/1‑100 (MQ=255)
                                                            tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                  >  1:20225/1‑100 (MQ=255)
                                                             aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                 <  1:1458118/100‑1 (MQ=255)
                                                              tttttCGTTCGACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                >  1:2542084/1‑100 (MQ=255)
                                                                   cgttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                           >  1:132056/1‑100 (MQ=255)
                                                                    gttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                          <  2:596551/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcg                         <  2:175576/100‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                          <  2:2212691/99‑1 (MQ=255)
                                                                     ttcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                           >  2:1284645/1‑98 (MQ=255)
                                                                      tcgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                        >  2:2269239/1‑100 (MQ=255)
                                                                       cgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                       >  1:142446/1‑100 (MQ=255)
                                                                       cgACCTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                       >  1:2581908/1‑100 (MQ=255)
                                                                           cTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                     <  2:2449215/98‑1 (MQ=255)
                                                                           cTTTGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                   <  1:2918811/100‑1 (MQ=255)
                                                                            tttGCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCa                  <  1:2695805/100‑1 (MQ=255)
                                                                               gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg               <  2:548588/100‑1 (MQ=255)
                                                                               gCTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATgg               >  1:855194/1‑100 (MQ=255)
                                                                                cTTATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc              >  1:1640960/1‑100 (MQ=255)
                                                                                 ttATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg             >  1:1064481/1‑100 (MQ=255)
                                                                                 ttATACATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg             >  1:805010/1‑100 (MQ=255)
                                                                                  tataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTTAAACAACTCAAGACAATAGCACTGTACAGCACTAAGCTTAAATTCCACGCCAATCAtctag             >  2:1747310/1‑95 (MQ=255)
                                                                                   ataCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg           <  1:2717531/100‑1 (MQ=255)
                                                                                    taCATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg           >  1:2743213/1‑100 (MQ=255)
                                                                                    taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc          >  1:569414/1‑100 (MQ=255)
                                                                                    taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc          <  2:1206276/100‑1 (MQ=255)
                                                                                    taCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc          >  2:2637813/1‑100 (MQ=255)
                                                                                     aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca         <  2:425217/100‑1 (MQ=255)
                                                                                     aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca         >  1:2689565/1‑100 (MQ=255)
                                                                                     aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca         >  1:1769148/1‑100 (MQ=255)
                                                                                     aCATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca         <  1:1597858/100‑1 (MQ=255)
                                                                                      cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac        >  2:1213405/1‑100 (MQ=255)
                                                                                      cATAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac        <  2:1481822/100‑1 (MQ=255)
                                                                                      cATAACAGTCG‑CTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac        >  1:111299/1‑100 (MQ=255)
                                                                                        tAACAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa      >  1:622592/1‑100 (MQ=255)
                                                                                          aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaattta    <  1:2047387/100‑6 (MQ=37)
                                                                                          aCAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGc    <  2:211803/100‑1 (MQ=255)
                                                                                           cAGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt   <  2:2351309/100‑1 (MQ=255)
                                                                                            aGTCG‑TTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCACGCCAATCATGGCGCGCACAAGCta  >  2:2528567/1‑100 (MQ=255)
                                                                                                 |                                                                                               
TTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTA  >  NC_000913/2470307‑2470499

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