Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,520,504:1 +T 16.1% 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.1% 753.1 / 157.9 230intergenic (‑180/+167)yjhU/yjhFputative DNA‑binding transcriptional regulator; KpLE2 phage‑like element/putative transporter
Reads supporting (aligned to +/- strand):  ref base . (94/99);  new base T (21/16);  total (115/115)
Fisher's exact test for biased strand distribution p-value = 4.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.89e-02

Only 100 of 237 total aligned reads displayed.

ATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520405‑4520603
                                                                                                    |                                                                                                   
     tttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCg                                                                                                >  1:1171961/1‑100 (MQ=255)
     tttATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCg                                                                                                <  2:1141897/100‑1 (MQ=255)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                            <  1:2795120/100‑1 (MQ=255)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCtt                                                                                            <  2:1920551/100‑1 (MQ=255)
         taCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCtt                                                                                            <  1:3790756/100‑1 (MQ=38)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  1:4298643/100‑1 (MQ=255)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  1:4366636/100‑1 (MQ=255)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCttt                                                                                           <  1:1792081/100‑1 (MQ=255)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  2:1865170/100‑1 (MQ=38)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  1:1443811/100‑1 (MQ=38)
          aCAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCttt                                                                                           <  1:2288719/100‑1 (MQ=38)
           naGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  1:3139242/2‑100 (MQ=255)
           cAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTa                                                                                          >  1:3646712/1‑100 (MQ=255)
            aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  2:195055/100‑1 (MQ=255)
            aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  2:481475/100‑1 (MQ=255)
            aGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTAc                                                                                         <  1:13651/100‑1 (MQ=255)
             gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACa                                                                                        >  2:452859/1‑100 (MQ=255)
             gATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACa                                                                                        <  2:2206060/100‑1 (MQ=38)
              atatTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  2:2280689/100‑1 (MQ=38)
              atatTTTATTCTGTTACATGAAGTAAACCAGCATACCATTTGGGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAg                                                                                       <  1:1281781/100‑1 (MQ=25)
                atTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTa                                                                                     >  2:957686/1‑100 (MQ=255)
                 ttttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTAt                                                                                    >  2:3112009/1‑100 (MQ=255)
                   ttATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTa                                                                                  >  2:723491/1‑100 (MQ=255)
                    tATTCTTTTTCATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTCGTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                                 <  1:247941/100‑1 (MQ=255)
                    tATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTAc                                                                                 >  2:1783866/1‑100 (MQ=255)
                        cTTTTACATGAAATAAACCAGCATCCCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             <  2:3144359/100‑1 (MQ=255)
                        cTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTgag                                                                             >  2:3650660/1‑100 (MQ=255)
                          tttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGaa                                                                           <  2:2141775/100‑1 (MQ=255)
                           ttACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAAt                                                                          <  2:571738/100‑1 (MQ=255)
                            tACATGAAATAAACCAGCAGTCCAGTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGGATTTCTGAGAATa                                                                         <  1:185137/100‑1 (MQ=255)
                              cATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                       >  2:4163884/1‑100 (MQ=255)
                              cATGAAATAAAACAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCAGACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATaca                                                                       <  1:2172729/100‑1 (MQ=255)
                               aTGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATacac                                                                      >  1:885981/1‑100 (MQ=255)
                                tGAAATAAACCAGCATACCATTTGAGCGTTCTCCTGCTGTTCCCCCATACTGTATTGCCGGAGACATGGTCCGGCCTTACAGTATTACTGAGAAt                                                                          <  1:3898339/95‑1 (MQ=255)
                                  aaaTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCa                                                                   >  2:1353392/1‑100 (MQ=255)
                                    aTAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATc                                                                 <  2:4029325/100‑1 (MQ=255)
                                     tAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCa                                                                <  2:4369904/100‑1 (MQ=255)
                                       aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                              <  2:3694973/100‑1 (MQ=255)
                                       aaCCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAg                                                              >  1:425982/1‑100 (MQ=255)
                                          cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                           <  2:4048184/100‑1 (MQ=255)
                                          cAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCAt                                                           >  1:3448212/1‑100 (MQ=255)
                                            gCATACCATTTGAGCGATCGCCCTCTGTTCCCCCATCCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATa                                                          <  2:2075174/99‑1 (MQ=255)
                                            gCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATaa                                                         <  1:3349217/100‑1 (MQ=255)
                                                 ccATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGaa                                                    <  1:152071/100‑1 (MQ=255)
                                                  cATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAg                                                   <  1:2518831/100‑1 (MQ=255)
                                                      taaGCGATCGCCCGCTGTCCCCCCATACTGGAAAGCCGCACACATGCTCCGGCTTTAAAGTATTACTGATAATCCCCCCATCAAGCATAAGCGAAGGtt                                                <  1:1551041/97‑1 (MQ=255)
                                                      tGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTc                                               <  1:1808991/100‑1 (MQ=255)
                                                       gAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCa                                              >  2:3491683/1‑100 (MQ=255)
                                                         gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            >  1:4071711/1‑100 (MQ=255)
                                                         gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            >  2:2698334/1‑100 (MQ=255)
                                                         gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            >  2:3422004/1‑100 (MQ=255)
                                                         gCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAat                                            >  1:1411473/1‑100 (MQ=255)
                                                           gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAatat                                          >  1:3053569/1‑100 (MQ=255)
                                                           gATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACAACCATCAAGCATAATCGAAGGTTCAatat                                          >  1:3560100/1‑100 (MQ=255)
                                                            aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                         >  1:551586/1‑100 (MQ=255)
                                                            aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                         >  2:2858634/1‑100 (MQ=255)
                                                            aTCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATc                                         >  2:2516791/1‑100 (MQ=255)
                                                             tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                        <  2:2789475/100‑1 (MQ=255)
                                                             tCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCg                                        <  1:765567/100‑1 (MQ=255)
                                                              cGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGc                                       >  1:853175/1‑100 (MQ=255)
                                                               gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      >  2:20189/1‑100 (MQ=255)
                                                               gCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGcc                                      <  1:3185543/100‑1 (MQ=255)
                                                                cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGGAGGTTCAATATCGc                                       <  1:346182/98‑1 (MQ=255)
                                                                cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  2:1041086/1‑100 (MQ=255)
                                                                cccGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGccc                                     >  2:2335187/1‑100 (MQ=255)
                                                                     tGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATTCACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCtttt                                <  1:2393414/100‑1 (MQ=255)
                                                                      tctcagtaaTACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCttttt                               >  2:3283974/9‑100 (MQ=21)
                                                                       ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              <  2:144287/100‑1 (MQ=255)
                                                                       ttCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              <  1:1591248/100‑1 (MQ=255)
                                                                       ttCCCCCATACTGTAAAGCCGCATACATGCTCCGTGTTTACAGTATTACTGATAATATACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTg                              <  1:2645696/100‑1 (MQ=255)
                                                                        tCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGa                             <  2:1341325/100‑1 (MQ=255)
                                                                         cccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGAt                            <  2:1481114/100‑1 (MQ=255)
                                                                          ccccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGan                            <  2:1607254/99‑2 (MQ=255)
                                                                          ccccATACTGTAAAGACGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATc                           >  1:2061642/1‑100 (MQ=255)
                                                                            ccATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGc                         >  1:1264333/1‑100 (MQ=255)
                                                                             caggcggTAAAGCCGGATACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCg                          <  2:1902578/92‑1 (MQ=255)
                                                                             cATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGcc                        >  2:3545586/1‑100 (MQ=255)
                                                                               tACTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  2:1688737/1‑100 (MQ=37)
                                                                                aCTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAg                      >  1:937768/1‑99 (MQ=255)
                                                                                 cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  1:3978099/1‑100 (MQ=38)
                                                                                 cTGTAAAGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCg                    >  2:3627567/1‑100 (MQ=38)
                                                                                   gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  2:3605030/1‑100 (MQ=255)
                                                                                   gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  <  2:78941/100‑1 (MQ=255)
                                                                                   gTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAt                  >  2:4013744/1‑100 (MQ=255)
                                                                                    tAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGAtt                 <  1:620915/100‑1 (MQ=255)
                                                                                     aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:2334644/1‑100 (MQ=255)
                                                                                     aaaGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:3908637/1‑100 (MQ=255)
                                                                                     aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:3720126/1‑100 (MQ=38)
                                                                                     aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  1:2734013/1‑100 (MQ=38)
                                                                                     aaaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTc                >  2:2498852/1‑100 (MQ=38)
                                                                                      aaGCCGGAGCATGTCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCg               >  2:83870/1‑100 (MQ=38)
                                                                                           ggAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCAcc          >  2:2240318/1‑100 (MQ=255)
                                                                                             agaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGt        <  1:1349485/100‑1 (MQ=255)
                                                                                              gaCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGt        >  1:4370054/1‑99 (MQ=255)
                                                                                               aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  1:4205013/100‑1 (MQ=255)
                                                                                               aCATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGc      <  2:3453112/100‑1 (MQ=255)
                                                                                                cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  1:1994709/1‑100 (MQ=255)
                                                                                                cATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCa     >  1:459479/1‑100 (MQ=255)
                                                                                                  tGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGCTCGCCAGCGATTCGTCACCGTGCa     >  1:163027/1‑98 (MQ=255)
                                                                                                   gCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAc  >  1:3384676/1‑100 (MQ=255)
                                                                                                    |                                                                                                   
ATGATTTTATACAGATATTTTATTCTTTTACATGAAATAAACCAGCATACCATTTGAGCGATCGCCCGCTGTTCCCCCATACTGTAAAGCCGGAGACATGCTCCGGCTTTACAGTATTACTGAGAATACACCCATCAAGCATAATCGAAGGTTCAATATCGCCCCTTTTTGATCGCCAGCGATTCGTCACCGTGCAGAC  >  NC_000913/4520405‑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: