Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,500 Δ1 bp 18.1% intergenic (‑176/+171) yjhU ← / ← yjhF putative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,520,5000A.18.1% 394.0 / 85.7 138intergenic (‑176/+171)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base A (65/48);  new base . (12/13);  total (77/61)
Fisher's exact test for biased strand distribution p-value = 5.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

Only 100 of 142 total aligned reads displayed.

GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  NC_000913/4520403‑4520598
                                                                                                 |                                                                                                   
 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                  <  2:847379/100‑1 (MQ=255)
         tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                         <  2:3444975/100‑1 (MQ=255)
          ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCt                                                                                        <  2:1698542/100‑1 (MQ=38)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                       >  1:3079780/1‑100 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                       <  2:3339849/100‑1 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                       <  1:2027020/100‑1 (MQ=255)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                      <  2:1330498/100‑1 (MQ=255)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                      <  1:3215406/100‑1 (MQ=38)
             cAGTTATTTTATTCTGTTACATGACATAAACCAGCCGACCATTTGCGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGc                                                                                         <  2:3321553/96‑1 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTTAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     >  2:2906664/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     >  2:107047/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     <  1:2746184/100‑1 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATGTGAGCGATCGCCCGCTGTTCCCCCATACAGTAAAGCCGGAGCATGTCTCCggaatta                                                                                     <  2:1092708/100‑6 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAACATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTATAAAGCCGGAGACATGCTCCGGCTTTa                                                                                     >  2:1017216/1‑100 (MQ=255)
              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                    <  1:243310/100‑1 (MQ=255)
              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTAc                                                                                    <  1:1731486/100‑1 (MQ=38)
               gATATTTTATTCTTTTACATGACATAAACCACCATACGATTTGAGCGGTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                   <  2:2461398/100‑1 (MQ=25)
               gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                   <  1:1892695/100‑1 (MQ=38)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                  <  1:2946939/100‑1 (MQ=38)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                  <  2:1098705/100‑1 (MQ=38)
                 tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                   <  1:1245595/98‑1 (MQ=38)
                  atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                >  1:3102914/1‑100 (MQ=255)
                  atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                >  1:731538/1‑100 (MQ=255)
                   ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAt                                                                               >  1:3469785/1‑100 (MQ=255)
                      tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                            >  1:1112156/1‑100 (MQ=255)
                       aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTAtgggg                                                                           <  1:3080790/100‑5 (MQ=25)
                        ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                          <  1:302076/100‑1 (MQ=255)
                         tCTTTTACCTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTga                                                                         >  2:3243788/1‑100 (MQ=255)
                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTga                                                                         <  1:3341945/100‑1 (MQ=255)
                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                          >  1:3361999/1‑99 (MQ=255)
                          cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                        <  1:3059356/100‑1 (MQ=255)
                          cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTg                                                                          >  2:346789/1‑98 (MQ=255)
                            tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                      >  1:2437231/1‑100 (MQ=255)
                               aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATac                                                                   >  1:49185/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  1:2202898/1‑100 (MQ=255)
                                 aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                 >  2:1752573/1‑100 (MQ=255)
                                 aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAGGCTACGGCCTTACACTATTACTGAGAATaca                                                                  >  2:393325/1‑99 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGTGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                >  2:3024862/1‑100 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  1:1734784/100‑1 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  2:40964/100‑1 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  2:1178406/100‑1 (MQ=255)
                                     aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                             >  1:2385953/1‑100 (MQ=255)
                                       tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGGGAATACACCCATCa                                                           <  1:2277764/100‑1 (MQ=255)
                                        aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                          <  1:3244911/100‑1 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                       <  2:2095246/100‑1 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                       <  2:3622547/100‑1 (MQ=255)
                                               cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAt                                                   >  1:263646/1‑100 (MQ=255)
                                                aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                  <  2:2103714/100‑1 (MQ=255)
                                                   ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                               >  1:1622078/1‑100 (MQ=255)
                                                     aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                             >  1:1523919/1‑100 (MQ=255)
                                                        tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGAAGACATGCTCCGGCTTTACAGTATTACTGAGAATACCCCCATCAAGCATAATCGAAGGTTc                                          <  2:3064338/100‑1 (MQ=255)
                                                           gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                        >  1:3075243/1‑99 (MQ=255)
                                                           gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                       >  2:3094302/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  1:2806483/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  1:194483/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  2:762049/1‑100 (MQ=255)
                                                              aTCGCCCGCTGTTCCCCCATACTGTACAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                    >  1:2212440/1‑100 (MQ=255)
                                                               tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                   <  2:276407/100‑1 (MQ=255)
                                                                cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                  <  2:3404768/100‑1 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACATTATTACTGAGAATACACCCATCAAGCATAATCTAAGTTTCAATATCGcc                                 <  1:2836733/100‑1 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                 >  1:1038678/1‑100 (MQ=255)
                                                                   ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcccc                               >  1:440341/1‑100 (MQ=255)
                                                                    cGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcccc                               >  1:2777304/1‑99 (MQ=255)
                                                                    cGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                              >  1:2938202/1‑100 (MQ=255)
                                                                       tGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGATAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                           <  2:1617857/100‑1 (MQ=255)
                                                                        gtccccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCAACAAGCATAATCGAAGGTTCAATATCGCCCCttttt                          <  1:1339634/97‑1 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  1:656403/1‑100 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  2:2642745/1‑100 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        <  2:728609/100‑1 (MQ=255)
                                                                           cccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                       >  1:2869308/1‑100 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  2:1650004/1‑98 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                       <  2:1621562/99‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      <  2:2097622/100‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      <  1:407255/100‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      >  1:805526/1‑100 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     <  1:1374043/100‑1 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  1:2623079/1‑100 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  1:3333574/1‑100 (MQ=255)
                                                                                  aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                <  1:3446465/100‑1 (MQ=255)
                                                                                   cTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:3524380/1‑100 (MQ=255)
                                                                                   cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:656959/1‑100 (MQ=38)
                                                                                   cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:444636/1‑100 (MQ=38)
                                                                                    tGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa              >  2:68737/1‑100 (MQ=38)
                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  2:1559326/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:2582759/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:3051447/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:65565/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATAGCCAGCGAt             >  1:2718260/1‑100 (MQ=38)
                                                                                     gTAAAGCCGGAGCATGTCTCCGGCTTTACAATACTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGACCGCCAGCGAt             >  2:2394062/1‑100 (MQ=37)
                                                                                      tAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt            <  1:22859/100‑1 (MQ=255)
                                                                                       aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  2:939502/1‑100 (MQ=255)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:3450254/1‑100 (MQ=38)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:1710705/1‑100 (MQ=38)
                                                                                       aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:1583491/1‑100 (MQ=38)
                                                                                        aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg          >  1:3312889/1‑100 (MQ=255)
                                                                                        aaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg          >  1:1703006/1‑100 (MQ=255)
                                                                                         aGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt         <  1:497342/100‑1 (MQ=255)
                                                                                          gCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  1:990878/1‑100 (MQ=255)
                                                                                          gCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  2:18737/1‑100 (MQ=39)
                                                                                            cGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc      >  2:2566566/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  NC_000913/4520403‑4520598

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: