Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 17.5% intergenic (‑180/+167) 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,5041.T17.5% 445.4 / 107.4 143intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (69/49);  new base T (12/13);  total (81/62)
Fisher's exact test for biased strand distribution p-value = 3.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.32e-01

Only 100 of 146 total aligned reads displayed.

TTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520409‑4520603
                                                                                                |                                                                                                   
ttttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGCGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTcc                                                                                                 <  1:3455145/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)
     taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                            <  1:4798/100‑1 (MQ=38)
      aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                           <  2:1330498/100‑1 (MQ=255)
      aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttn                                                                                           <  1:3240244/100‑2 (MQ=255)
       cAGTTATTTTATTCTGTTACATGACATAAACCAGCCGACCATTTGCGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                              <  2:3321553/96‑1 (MQ=255)
       cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTTAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                          >  2:2906664/1‑100 (MQ=255)
       cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                          >  1:2191827/1‑100 (MQ=255)
       cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                          >  2:107047/1‑100 (MQ=255)
       cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTa                                                                                          <  1:3402931/100‑1 (MQ=38)
       cAGATATTTTATTCTTTTACATGAAATAAACCAACATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTATAAAGCCGGAGACATG‑CTCCGGCTTTa                                                                                          >  2:1017216/1‑100 (MQ=255)
        aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTAc                                                                                         <  1:243310/100‑1 (MQ=255)
         gATATTTTATTCTTTTACATGACATAAACCACCATACGATTTGAGCGGTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                        <  2:2461398/100‑1 (MQ=25)
         gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                        >  1:1726526/1‑100 (MQ=255)
         gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                        <  1:1892695/100‑1 (MQ=38)
          atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                       <  1:3093287/100‑1 (MQ=38)
          atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                       <  1:2946939/100‑1 (MQ=38)
          atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                       <  2:1098705/100‑1 (MQ=38)
            atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                     >  1:3102914/1‑100 (MQ=255)
             ttttATTCTTTTACCTGAAATAAACCAGCATACCATTTTTGCTATCGCCCGCTGTTGCTCCATACTGTCAAGCCGGAGACATG‑CTCCTGCTTTACAGTAt                                                                                    >  2:777378/1‑100 (MQ=255)
             ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTAt                                                                                    >  1:262920/1‑100 (MQ=255)
                tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTAc                                                                                 >  1:1112156/1‑100 (MQ=255)
                  ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                               <  1:302076/100‑1 (MQ=255)
                   tCTTTTACCTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                              >  2:3243788/1‑100 (MQ=255)
                   tCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                               >  1:3361999/1‑99 (MQ=255)
                    cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                               >  2:346789/1‑98 (MQ=255)
                      tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                           >  1:2437231/1‑100 (MQ=255)
                      tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                           <  1:1377526/100‑1 (MQ=255)
                       ttACATAAAATAAACCAGCATACCATTTGAGCGAGCGACCGCTGTTCCCCCATACTGTAAAGCCGGAGACATT‑CTCCGGCTTTACAGTATTACTGAGAAt                                                                          <  2:1065204/100‑1 (MQ=39)
                         aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                        >  1:276409/1‑100 (MQ=255)
                         aCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATac                                                                        >  1:49185/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:3068486/100‑1 (MQ=255)
                            tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     <  1:1734784/100‑1 (MQ=255)
                            tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                     <  2:40964/100‑1 (MQ=255)
                               aaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAt                                                                  >  1:2385953/1‑100 (MQ=255)
                                aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATc                                                                 <  1:2259924/100‑1 (MQ=255)
                                  aaaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCaa                                                               <  1:3244911/100‑1 (MQ=255)
                                     ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                            <  2:3622547/100‑1 (MQ=255)
                                     ccAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCa                                                            <  2:2095246/100‑1 (MQ=255)
                                        gCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                         >  2:1369976/1‑100 (MQ=255)
                                         cATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAAt                                                        >  1:263646/1‑100 (MQ=255)
                                          aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                       <  2:2103714/100‑1 (MQ=255)
                                          aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                       <  2:1123526/100‑1 (MQ=255)
                                             ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                    >  1:3328597/1‑100 (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)
                                                        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)
                                                                 tGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGATAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                                <  2:1617857/100‑1 (MQ=255)
                                                                  gtccccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCAACAAGCATAATCGAAGGTTCAATATCGCCCCttttt                               <  1:1339634/97‑1 (MQ=255)
                                                                    tCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  2:728609/100‑1 (MQ=255)
                                                                     cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            >  1:909482/1‑100 (MQ=255)
                                                                     cccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            >  1:2809441/1‑100 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             >  2:1650004/1‑98 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            <  2:1621562/99‑1 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           <  1:407255/100‑1 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           <  1:2900702/100‑1 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           >  1:805526/1‑100 (MQ=255)
                                                                      ccccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           <  2:2097622/100‑1 (MQ=255)
                                                                       cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          >  1:3333574/1‑100 (MQ=255)
                                                                       cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          <  1:186447/100‑1 (MQ=255)
                                                                             cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  1:3524380/1‑100 (MQ=255)
                                                                             cTGTAAAGCCGGAG‑TATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGCTCAATATCCCCCCTTTTTGACCACCAGc                     >  2:716526/1‑99 (MQ=25)
                                                                             cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  1:656959/1‑100 (MQ=38)
                                                                              tGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                   >  1:694447/1‑100 (MQ=38)
                                                                               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)
                                                                                    gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc             >  2:18737/1‑100 (MQ=39)
                                                                                      cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAc           >  2:2566566/1‑100 (MQ=255)
                                                                                        gagaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCg         <  1:1622454/100‑1 (MQ=255)
                                                                                        ccgatctG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCg         <  1:2552325/94‑1 (MQ=255)
                                                                                          gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg       >  2:2732911/1‑100 (MQ=255)
                                                                                          gaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTg       >  1:1278354/1‑100 (MQ=255)
                                                                                            cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  1:2909360/1‑100 (MQ=255)
                                                                                              tG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAg    >  1:3031895/1‑99 (MQ=255)
                                                                                              tG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGa   >  1:1061957/1‑100 (MQ=255)
                                                                                               g‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATGGCCCCTGTTTGATCGCCAGCGATTCGTCACCgagatcgg  >  2:2552309/1‑93 (MQ=255)
                                                                                               g‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  1:624129/1‑100 (MQ=255)
                                                                                                |                                                                                                   
TTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520409‑4520603

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: ‑