Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 16.9% 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.T16.9% 395.8 / 87.0 124intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (60/43);  new base T (8/13);  total (68/56)
Fisher's exact test for biased strand distribution p-value = 9.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 126 total aligned reads displayed.

ATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGA  >  NC_000913/4520405‑4520602
                                                                                                    |                                                                                                  
aTGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATg                                                                                                     >  2:1918563/1‑100 (MQ=255)
   aTTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACAATTTGAGCGATAGCAAGCTGTTACCCCATACTGTAAAGCCGGAGACATg                                                                                                     >  2:2447217/1‑97 (MQ=255)
      ttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCgg                                                                                              <  2:2911167/100‑1 (MQ=255)
       tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                             >  1:3026007/1‑100 (MQ=255)
       tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGc                                                                                             <  2:1856768/100‑1 (MQ=255)
       tataCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGc                                                                                             <  2:794034/100‑1 (MQ=39)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                           <  2:104223/100‑1 (MQ=255)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCtt                                                                                           <  2:795995/100‑1 (MQ=255)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                           <  1:2215314/100‑1 (MQ=38)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                           <  2:796265/100‑1 (MQ=38)
          aCATATATTTTATTCTTTTACATGAAATAAACCATCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                          <  2:929328/100‑1 (MQ=37)
          aCAGATATTTTATTCTTTTACATGATATAAACCAGCATACCATTTGATCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCtt                                                                                           <  2:1695867/99‑1 (MQ=37)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCttt                                                                                          <  1:1704532/100‑1 (MQ=255)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCttt                                                                                          <  2:325829/100‑1 (MQ=38)
             tcttccgatctCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                       <  1:1950119/90‑1 (MQ=25)
             gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACa                                                                                       <  1:2717882/100‑1 (MQ=255)
             gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                       <  1:1737504/100‑1 (MQ=38)
             gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACa                                                                                       <  1:567151/100‑1 (MQ=38)
              atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                      <  1:910657/100‑1 (MQ=38)
              atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAg                                                                                      <  1:1457707/100‑1 (MQ=38)
               tatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACAAA‑CTCCGGCTTTACAGt                                                                                     >  1:2924010/1‑100 (MQ=255)
                atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTa                                                                                    >  1:908986/1‑100 (MQ=255)
                actttATTCTTTTACATGAAATAAACCAGCATACCATTTGGTCGATCGCCCGCTTTTCCCCCATACTGTAAAGCCTGATACATG‑CTCCGGCTTTACAGt                                                                                     <  2:1481318/97‑1 (MQ=255)
                 ttttATTCTTTTACATGAAATAAACCAGCATGCCATTTGAGCGATCGCCCGCTTTTCCCCCATACTGTAAAGTCGTAGACATG‑CTCCGGCTTTACAGTa                                                                                    <  2:2682777/99‑1 (MQ=255)
                      ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACt                                                                               >  2:2184237/1‑99 (MQ=255)
                      ttCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAACCGGAGACATG‑CTCCGGCTTTACAGTATTACTg                                                                              >  1:2600720/1‑100 (MQ=255)
                       tCTTTTACATGAAATCAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTga                                                                             <  1:3000453/100‑1 (MQ=255)
                        cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                            >  1:122035/1‑100 (MQ=255)
                        cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgag                                                                            >  2:472127/1‑100 (MQ=255)
                        cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTAtgggggaa                                                                            <  2:2930499/100‑8 (MQ=21)
                         ttttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                           <  1:2840008/100‑1 (MQ=255)
                         ttttACATGAAATAAACCAGAATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTgaga                                                                           >  2:161558/1‑100 (MQ=255)
                          tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAAACATG‑CTCCGGCTTTACAGTATTACTGAGaa                                                                          >  1:502247/1‑100 (MQ=255)
                            taactGAACTAAACCAGCATACCATTTGTGCGATCGCCCGCTGTTGCCCCATACTGTAAAGCGTTAGACATG‑CTCCGGTTTTACAGTATGACTgag                                                                            <  2:2904867/92‑1 (MQ=255)
                              cATGAAATAAACCAGTATACCATTTGAGCGATCGCTCACTGTTCCCCCATACTGTAAATCCTGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:2660536/1‑100 (MQ=255)
                              cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACGGTATTACTGAGAATaca                                                                      >  1:2170524/1‑100 (MQ=255)
                              cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:1474810/1‑100 (MQ=255)
                              cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATaca                                                                      >  2:959972/1‑100 (MQ=255)
                              cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCACCATACTGTAAAGCCGGAGACCTG‑CTACGGCTTTACAGTATTATTGAGAATcca                                                                      >  2:245667/1‑100 (MQ=255)
                                tGAAATAAACCAGCATACCATTTGTGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACAcc                                                                    <  1:1290244/100‑1 (MQ=255)
                                tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAATCCGGAGACATG‑CTCCAGCTTTACAGTATTACTGAGAATACAcc                                                                    <  2:1136347/100‑1 (MQ=255)
                                tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATacac                                                                     <  2:407344/99‑1 (MQ=255)
                                tGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CGCCGGCTTTACAGTATTACTGAGAATACAcc                                                                    <  2:1996617/100‑1 (MQ=255)
                                  aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                                  >  1:2225646/1‑100 (MQ=255)
                                  aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                                  >  2:2254686/1‑100 (MQ=255)
                                  aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAATATTACTGAGAATACACCCa                                                                  >  1:632569/1‑100 (MQ=255)
                                   aaTACACCAGCATACCATTTGAGCGATCGCCCGCTGTTCACCCATACTGTAAAGCAGGAGACATA‑CTCCGCCTTTACAGTATTACAGAAAAAACACCCAt                                                                 >  2:2154147/1‑100 (MQ=255)
                                       aaCCAGCATACCATTTGAGCGATCGCCCGGTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                             <  1:2219203/100‑1 (MQ=255)
                                       aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCCTACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                             >  1:3033402/1‑100 (MQ=255)
                                       aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                             >  1:2331508/1‑100 (MQ=255)
                                        aCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCAGACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGc                                                            <  2:2331481/100‑1 (MQ=255)
                                              aTACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATc                                                      >  1:1496197/1‑100 (MQ=255)
                                               tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                     >  1:1082136/1‑100 (MQ=255)
                                               tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                     <  1:2195176/100‑1 (MQ=255)
                                               tACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                     >  2:1176147/1‑100 (MQ=255)
                                                aCCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCg                                                     <  2:2225623/99‑1 (MQ=255)
                                                aCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGa                                                    >  2:673321/1‑100 (MQ=255)
                                                 ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                   >  1:1671159/1‑100 (MQ=255)
                                                 ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                   >  1:2789301/1‑100 (MQ=255)
                                                   aTTTGAGCGATCGCCCGCTGTTCCTCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAgg                                                 <  2:533105/100‑1 (MQ=255)
                                                    tttGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                   >  2:1306386/1‑97 (MQ=255)
                                                      tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                              <  1:2706918/100‑1 (MQ=255)
                                                       gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                             >  1:2385567/1‑100 (MQ=255)
                                                        aGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCaa                                            <  2:230406/100‑1 (MQ=255)
                                                         gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                           >  1:2221446/1‑100 (MQ=255)
                                                              cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                      <  2:2113555/100‑1 (MQ=255)
                                                               gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAGTACACTCATCAAGCATAATCGAAGGTTCAATATCGc                                      <  2:2114975/99‑1 (MQ=255)
                                                               gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                     >  2:943851/1‑100 (MQ=255)
                                                               gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                     >  2:2905147/1‑100 (MQ=255)
                                                                cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATT‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                    <  1:902725/100‑1 (MQ=255)
                                                                cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                    <  2:2614236/100‑1 (MQ=255)
                                                                     ttctcagtaaTACTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                               <  1:1573612/91‑1 (MQ=21)
                                                                       ttCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                             <  2:2700587/100‑1 (MQ=255)
                                                                           cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                         >  1:3019787/1‑100 (MQ=255)
                                                                           cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                         <  1:2689461/100‑1 (MQ=255)
                                                                           cccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                         >  2:1705731/1‑100 (MQ=255)
                                                                            ccATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGACCCTTTTTGATCGc                        >  2:2067240/1‑100 (MQ=255)
                                                                             cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                       >  1:81511/1‑100 (MQ=255)
                                                                             cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                       <  2:2265585/100‑1 (MQ=255)
                                                                             cATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                       >  1:2907344/1‑100 (MQ=255)
                                                                               tACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                     >  1:2169233/1‑100 (MQ=255)
                                                                                aCTGTAAAGCCGGAG‑CGTGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGAATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                     >  1:586548/1‑99 (MQ=37)
                                                                                 cTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   >  2:768364/1‑100 (MQ=255)
                                                                                 cTGTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACACATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                   >  1:23834/1‑100 (MQ=37)
                                                                                  tGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGa                  >  2:170106/1‑100 (MQ=255)
                                                                                   gTAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCACCAGCGa                  >  1:1197632/1‑99 (MQ=37)
                                                                                    tAAAGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                >  2:1474302/1‑100 (MQ=38)
                                                                                     aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  2:3043967/1‑100 (MQ=255)
                                                                                     aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGATTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  1:2397734/1‑100 (MQ=255)
                                                                                     aaaGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCAAAGGTTCAATATCGCAACTTTTTGATCGCCATCGAtt                >  1:1362652/1‑99 (MQ=255)
                                                                                     aaaGCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc               >  2:2977919/1‑100 (MQ=38)
                                                                                        gCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCATTTTTGATCGCCAGCGATTCGTc            >  1:64006/1‑100 (MQ=255)
                                                                                        gCCGGAG‑CATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTc            >  1:2568817/1‑100 (MQ=39)
                                                                                          cGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATAGCCCCTTTTTGATCGCCAGCGATTCGTCAc          >  2:1165096/1‑100 (MQ=255)
                                                                                             agaCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATAGCCAGCGATTCGTCACCg        >  1:1481021/1‑99 (MQ=255)
                                                                                               aCATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc     <  1:803074/100‑1 (MQ=255)
                                                                                                cATG‑CTCCGTCTTTACAGTATTACTGAGAATACACTCATTAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa    <  1:287950/100‑1 (MQ=255)
                                                                                                cATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa    >  1:492485/1‑100 (MQ=255)
                                                                                                  tG‑CTCCGGCTTTACAGTATTACTGATAATACTCCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAg   <  1:2533447/99‑1 (MQ=255)
                                                                                                  tG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTTAATATCGCCCCTTTTTGATCCCCAGCGATTCGTCACCGTGCAGa  >  2:3032629/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
ATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATG‑CTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGA  >  NC_000913/4520405‑4520602

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