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.

GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  NC_000913/4520403‑4520598
                                                                                                 |                                                                                                   
 cATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAt                                                                                                  <  2:847379/100‑1 (MQ=255)
         tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                         <  2:3444975/100‑1 (MQ=255)
          ataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCt                                                                                        <  2:1698542/100‑1 (MQ=38)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                       >  1:3079780/1‑100 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                       <  2:3339849/100‑1 (MQ=255)
           taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                       <  1:2027020/100‑1 (MQ=255)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                      <  2:1330498/100‑1 (MQ=255)
            aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                      <  1:3215406/100‑1 (MQ=38)
             cAGTTATTTTATTCTGTTACATGACATAAACCAGCCGACCATTTGCGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                         <  2:3321553/96‑1 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTTAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                     >  2:2906664/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                     >  2:107047/1‑100 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                     <  1:2746184/100‑1 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATGTGAGCGATCGCCCGCTGTTCCCCCATACAGTAAAGCCGGAG‑CATGTCTCCggaatta                                                                                     <  2:1092708/100‑6 (MQ=255)
             cAGATATTTTATTCTTTTACATGAAATAAACCAACATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTATAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                     >  2:1017216/1‑100 (MQ=255)
              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                    <  1:243310/100‑1 (MQ=255)
              aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTAc                                                                                    <  1:1731486/100‑1 (MQ=38)
               gATATTTTATTCTTTTACATGACATAAACCACCATACGATTTGAGCGGTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                   <  2:2461398/100‑1 (MQ=25)
               gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                   <  1:1892695/100‑1 (MQ=38)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                  <  1:2946939/100‑1 (MQ=38)
                atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                  <  2:1098705/100‑1 (MQ=38)
                 tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                   <  1:1245595/98‑1 (MQ=38)
                  atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                >  1:3102914/1‑100 (MQ=255)
                  atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                >  1:731538/1‑100 (MQ=255)
                   ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                               >  1:3469785/1‑100 (MQ=255)
                      tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                            >  1:1112156/1‑100 (MQ=255)
                       aTTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAtgggg                                                                           <  1:3080790/100‑5 (MQ=25)
                        ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                          <  1:302076/100‑1 (MQ=255)
                         tCTTTTACCTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                         >  2:3243788/1‑100 (MQ=255)
                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                         <  1:3341945/100‑1 (MQ=255)
                         tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                          >  1:3361999/1‑99 (MQ=255)
                          cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                        <  1:3059356/100‑1 (MQ=255)
                          cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                          >  2:346789/1‑98 (MQ=255)
                            tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                      >  1:2437231/1‑100 (MQ=255)
                               aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                   >  1:49185/1‑100 (MQ=255)
                                cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                  >  1:2202898/1‑100 (MQ=255)
                                 aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATacac                                                                 >  2:1752573/1‑100 (MQ=255)
                                 aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAGG‑CTACGGCCTTACACTATTACTGAGAATaca                                                                  >  2:393325/1‑99 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGTGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                >  2:3024862/1‑100 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  1:1734784/100‑1 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  2:40964/100‑1 (MQ=255)
                                  tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                <  2:1178406/100‑1 (MQ=255)
                                     aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                             >  1:2385953/1‑100 (MQ=255)
                                       tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGGGAATACACCCATCa                                                           <  1:2277764/100‑1 (MQ=255)
                                        aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                          <  1:3244911/100‑1 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                       <  2:2095246/100‑1 (MQ=255)
                                           ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                       <  2:3622547/100‑1 (MQ=255)
                                               cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAt                                                   >  1:263646/1‑100 (MQ=255)
                                                aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                  <  2:2103714/100‑1 (MQ=255)
                                                   ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                               >  1:1622078/1‑100 (MQ=255)
                                                     aTTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATT‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                             >  1:1523919/1‑100 (MQ=255)
                                                        tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGAAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACCCCCATCAAGCATAATCGAAGGTTc                                          <  2:3064338/100‑1 (MQ=255)
                                                           gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                        >  1:3075243/1‑99 (MQ=255)
                                                           gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                       >  2:3094302/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  1:2806483/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  1:194483/1‑100 (MQ=255)
                                                             gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                     >  2:762049/1‑100 (MQ=255)
                                                              aTCGCCCGCTGTTCCCCCATACTGTACAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                    >  1:2212440/1‑100 (MQ=255)
                                                               tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                   <  2:276407/100‑1 (MQ=255)
                                                                cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                  <  2:3404768/100‑1 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACATTATTACTGAGAATACACCCATCAAGCATAATCTAAGTTTCAATATCGcc                                 <  1:2836733/100‑1 (MQ=255)
                                                                 gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                 >  1:1038678/1‑100 (MQ=255)
                                                                   ccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcccc                               >  1:440341/1‑100 (MQ=255)
                                                                    cGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcccc                               >  1:2777304/1‑99 (MQ=255)
                                                                    cGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCt                              >  1:2938202/1‑100 (MQ=255)
                                                                       tGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGATAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                           <  2:1617857/100‑1 (MQ=255)
                                                                        gtccccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAACAAGCATAATCGAAGGTTCAATATCGCCCCttttt                          <  1:1339634/97‑1 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  1:656403/1‑100 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  2:2642745/1‑100 (MQ=255)
                                                                          tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        <  2:728609/100‑1 (MQ=255)
                                                                           cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                       >  1:2869308/1‑100 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                        >  2:1650004/1‑98 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                       <  2:1621562/99‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      <  2:2097622/100‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      <  1:407255/100‑1 (MQ=255)
                                                                            ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                      >  1:805526/1‑100 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     <  1:1374043/100‑1 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  1:2623079/1‑100 (MQ=255)
                                                                             cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                     >  1:3333574/1‑100 (MQ=255)
                                                                                  aCTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGc                <  1:3446465/100‑1 (MQ=255)
                                                                                   cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:3524380/1‑100 (MQ=255)
                                                                                   cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:656959/1‑100 (MQ=38)
                                                                                   cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg               >  1:444636/1‑100 (MQ=38)
                                                                                    tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa              >  2:68737/1‑100 (MQ=38)
                                                                                     gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  2:1559326/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:2582759/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:3051447/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt             >  1:65565/1‑100 (MQ=255)
                                                                                     gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATAGCCAGCGAt             >  1:2718260/1‑100 (MQ=38)
                                                                                     gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAATACTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGACCGCCAGCGAt             >  2:2394062/1‑100 (MQ=37)
                                                                                      tAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt            <  1:22859/100‑1 (MQ=255)
                                                                                       aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  2:939502/1‑100 (MQ=255)
                                                                                       aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:3450254/1‑100 (MQ=38)
                                                                                       aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:1710705/1‑100 (MQ=38)
                                                                                       aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc           >  1:1583491/1‑100 (MQ=38)
                                                                                        aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg          >  1:3312889/1‑100 (MQ=255)
                                                                                        aaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg          >  1:1703006/1‑100 (MQ=255)
                                                                                         aGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGt         <  1:497342/100‑1 (MQ=255)
                                                                                          gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  1:990878/1‑100 (MQ=255)
                                                                                          gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc        >  2:18737/1‑100 (MQ=39)
                                                                                            cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc      >  2:2566566/1‑100 (MQ=255)
                                                                                                 |                                                                                                   
GCATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTG  >  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: ‑