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.95.3% 532.7 / 16.5 128coding (1274/1332 nt)gtrSserotype‑specific glucosyl transferase, CPS‑53 (KpLE1) prophage
Reads supporting (aligned to +/- strand):  ref base T (4/1);  new base . (55/67);  total (59/69)
Fisher's exact test for biased strand distribution p-value = 1.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 136 total aligned reads displayed.

ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATA  >  NC_000913/2470306‑2470501
                                                                                                  |                                                                                                 
    tATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGtttttt                                                                                              <  2:2072403/100‑1 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatt                                                                                          <  2:2252322/100‑5 (MQ=255)
        tgtCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGttttttaatn                                                                                          <  2:1807424/100‑5 (MQ=255)
           cATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAAAAGTCGTTTTTTAATTtat                                                                                      >  1:979630/1‑100 (MQ=255)
            tccgatctgggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTtata                                                                                     <  1:1257767/92‑1 (MQ=255)
               ttttGGGGAAGTATATAGAGATTTGAGCAAGCAGTTATTTGGATTGTATTTTTTGTTCGAACTTTGCTTATACATAACAGTCGttttttaatt                                                                                          <  1:2715831/93‑5 (MQ=255)
               tttCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaaa                                                                                  <  1:137654/100‑1 (MQ=255)
                  cGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                               >  1:550584/1‑100 (MQ=255)
                  cGGGGAAGGATATAGAGATTTAAGCAAGCATTAATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                               >  2:2765718/1‑100 (MQ=255)
                   ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTGTTTTTAGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataa                                                                               <  1:1247856/99‑1 (MQ=255)
                   ggggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAaataat                                                                              <  2:2038593/100‑1 (MQ=255)
                    gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaagatcgga                                                                             >  2:1257772/1‑93 (MQ=255)
                    gggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAAtt                                                                             >  1:1152115/1‑100 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATaaagatcgga                                                                            >  2:981511/1‑93 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                            >  2:1231589/1‑100 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                            <  2:1647381/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                            <  2:946925/100‑1 (MQ=255)
                     ggAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTc                                                                            <  2:2545013/100‑1 (MQ=255)
                      gAAGGATATAGTTGTTTAAGCAAGCATTTAGTTGGCATGTATTTTTCTTTCGAACTTTGTTTATGCATAGTAGTCGTCTTTTAATTTATAAAATAATTc                                                                            <  1:930157/99‑1 (MQ=255)
                       aaGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCaa                                                                          <  2:1665637/100‑1 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                        <  2:1824781/100‑1 (MQ=255)
                         ggATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAg                                                                        <  2:1699748/100‑1 (MQ=255)
                          gATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGa                                                                       <  1:2123867/100‑1 (MQ=255)
                           atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaa                                                                      >  2:2141123/1‑100 (MQ=255)
                           atatAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGaa                                                                      <  1:2868172/100‑1 (MQ=255)
                             ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTTTTTCTACTTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                    >  1:2543662/1‑100 (MQ=255)
                             ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTATTTTATAAAATAATTCAAAGAAAt                                                                    >  1:2543943/1‑100 (MQ=255)
                             ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                    <  2:2347652/100‑1 (MQ=255)
                             ataGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAAt                                                                    <  2:743540/100‑1 (MQ=255)
                              taGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATc                                                                   >  2:27869/1‑100 (MQ=255)
                               agagATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCa                                                                  >  2:1123526/1‑100 (MQ=255)
                                gagaTTTAAGCAAGCATTTATTTGGAATTTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCaa                                                                 <  2:2635653/100‑1 (MQ=255)
                                 agaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAg                                                                <  1:2246472/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                               >  2:840513/1‑100 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                               >  2:851754/1‑100 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                               <  2:1762148/100‑1 (MQ=255)
                                  gaTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGa                                                               <  2:1322444/100‑1 (MQ=255)
                                   aTTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                              <  2:501340/100‑1 (MQ=255)
                                     ttAAGCAACTATTTATTTGGAATGTATTTCTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGAc                                                              <  1:2722536/98‑1 (MQ=255)
                                      tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                           >  1:1950791/1‑100 (MQ=255)
                                      tAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAAt                                                           >  1:430706/1‑100 (MQ=255)
                                        agcaagcaTTTATTTGGAATGTATTTTTTTTTCTTTTTTTGATTATAAAAAACACTCGTATTTTATTTTGTGAATGAATTCAAGGAAATcgacaccgtag                                                         >  1:1548241/1‑90 (MQ=255)
                                        agcaagcaTTTATTTGGAATGTATTTTTTGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATa                                                          >  1:32705/1‑99 (MQ=255)
                                         gcaagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGc                                                        >  2:18271/1‑100 (MQ=255)
                                          caagcaTTTATTTTGAATGTATTTTTCGTTCGTCCTTAGCTTATACATAACAGTCGTTTTTTAATGTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                       <  1:1848633/100‑1 (MQ=255)
                                          caagcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCg                                                       >  2:2327868/1‑100 (MQ=255)
                                             gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                    >  2:145147/1‑100 (MQ=255)
                                             gcaTTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATg                                                    >  2:2729099/1‑100 (MQ=255)
                                              caTTTATTTGGAATGTATTTTTCGTTCGACCATTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGa                                                      <  1:287148/97‑1 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTGTAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                  <  1:1555279/100‑1 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                  >  2:2487347/1‑100 (MQ=255)
                                               atttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTa                                                  >  1:551216/1‑100 (MQ=255)
                                                tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                 <  1:940687/100‑1 (MQ=255)
                                                tttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaa                                                 <  2:658018/100‑1 (MQ=255)
                                                 ttatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTaaa                                                <  2:1152120/100‑1 (MQ=255)
                                                  tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                               <  1:665477/100‑1 (MQ=255)
                                                  tatttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                               <  2:1195566/100‑1 (MQ=255)
                                                  natttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAg                                               >  1:2785295/2‑100 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGc                                              <  2:2397159/100‑1 (MQ=255)
                                                   atttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTAGTTTTTTAATTTATAAAATAATTCAAAGAAATCGAGACAATAGCGATGTAAAGc                                              <  1:1777899/100‑1 (MQ=255)
                                                    tttGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCa                                             >  2:1934090/1‑100 (MQ=255)
                                                       ggAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTa                                          <  2:1986669/100‑1 (MQ=255)
                                                          aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                       <  1:2286975/100‑1 (MQ=255)
                                                          aTGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGt                                       <  1:795004/100‑1 (MQ=255)
                                                           tGTATTTTTCGTTTTACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACACTAGCGATGTAAAGCACTAAGtt                                      >  1:2427679/1‑100 (MQ=255)
                                                             tATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTa                                    >  2:1729539/1‑100 (MQ=255)
                                                              aTTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATACAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaa                                   >  2:2599917/1‑100 (MQ=255)
                                                                ttttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTaaa                                  <  2:378756/99‑1 (MQ=255)
                                                                 tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                >  2:751096/1‑100 (MQ=255)
                                                                 tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                >  1:135257/1‑100 (MQ=255)
                                                                 tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAAtt                                >  2:1740073/1‑100 (MQ=255)
                                                                 tttCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAACCACTAAATTTaaa                                  >  2:211232/1‑98 (MQ=255)
                                                                  ttcgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTg                               <  1:1729531/100‑1 (MQ=255)
                                                                   tcttccgatcTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCGAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgc                              <  2:2055857/91‑1 (MQ=255)
                                                                    cgttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcg                             <  2:1381516/100‑1 (MQ=255)
                                                                     gttcgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgc                            <  1:1631226/100‑1 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                          >  2:1636590/1‑99 (MQ=255)
                                                                        cgACCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGcc                         <  2:1726272/100‑1 (MQ=255)
                                                                         gACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTgcgcgc                          >  2:664007/1‑99 (MQ=255)
                                                                          aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                       >  1:2398683/1‑100 (MQ=255)
                                                                          aCCTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCaa                       <  2:618520/100‑1 (MQ=255)
                                                                            cTTTGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATc                     >  2:2353931/1‑100 (MQ=255)
                                                                               tGCTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                  <  1:2353916/100‑1 (MQ=255)
                                                                                gCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTAAAAGAAATCAAGACAATAGGGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATg                  >  1:1589877/1‑100 (MQ=255)
                                                                                 cTTATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgc                >  2:263961/1‑100 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCatgattggcg               <  1:161256/100‑10 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg               <  2:41525/100‑1 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg               <  1:693790/100‑1 (MQ=255)
                                                                                  ttATACATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcg               <  1:2729081/100‑1 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc              <  1:1406309/100‑1 (MQ=255)
                                                                                   tataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgc              >  2:2602976/1‑100 (MQ=255)
                                                                                    ataCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcg             <  2:470921/100‑1 (MQ=255)
                                                                                     taCATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGgcgcgc            >  2:2845484/1‑100 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGACATAAAGACAATATCGATATAAAGTACTAATTTTAAATTGCAAGCCAAATATGgcgtg             >  1:396788/1‑95 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGcac          >  2:1626798/1‑100 (MQ=255)
                                                                                       cATAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGca           >  1:1516528/1‑99 (MQ=255)
                                                                                         tAACAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACaa        >  1:2434580/1‑100 (MQ=255)
                                                                                          aaCAGTCGTTTTTTAATTTATAAAATAATCCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCATGATTGGCGCGcaatt        <  2:55712/99‑4 (MQ=37)
                                                                                            cAGTCGTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCt     <  2:327743/100‑1 (MQ=255)
                                                                                                cGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCtata  <  2:823876/100‑1 (MQ=255)
                                                                                                  |                                                                                                 
ATTTTATGTGTCATTTTTCGGGGAAGGATATAGAGATTTAAGCAAGCATTTATTTGGAATGTATTTTTCGTTCGACCTTTGCTTATACATAACAGTCGTTTTTTTAATTTATAAAATAATTCAAAGAAATCAAGACAATAGCGATGTAAAGCACTAAGTTTAAATTGCGCGCCAATCATGGCGCGCACAAGCTATA  >  NC_000913/2470306‑2470501

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