Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,519,033 C→T intergenic (+231/+305) insA7 → / ← yjhU IS1 family repressor protein InsA/putative DNA‑binding transcriptional regulator YjhU

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,519,0330CT98.5% 217.8 / ‑6.5 66intergenic (+231/+305)insA7/yjhUIS1 family repressor protein InsA/putative DNA‑binding transcriptional regulator YjhU
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base T (25/40);  total (25/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTAT  >  NC_000913/4518786‑4519278
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
cAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGt                                                                                                                                                                                                                                                    >  1:450927/1‑251 (MQ=255)
         gTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCAc                                                                                                                                                                                                                                           <  1:597323/251‑1 (MQ=255)
         gTCGGAAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTAGAGATGACGGATGGCTGGCCGTTGTATGAATCAc                                                                                                                                                                                                                                           >  1:224986/1‑251 (MQ=255)
           cGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACg                                                                                                                                                                                                                                          <  1:548164/250‑1 (MQ=255)
            ggTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGcc                                                                                                                                                                                                                                        <  2:108819/251‑1 (MQ=255)
                    cgcgCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAgg                                                                                                                                                                                                                                  <  2:330315/249‑1 (MQ=255)
                              aCCGGGCAGTGATGTGATTGTCTGAGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCAGGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAAGGAAAGCTGCACg                                                                                                                                                                                                                      <  2:259940/251‑1 (MQ=255)
                                      gtgatgtgatTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCGGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGc                                                                                                                                                                                                              <  1:448066/251‑1 (MQ=255)
                                        gatgtgatTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCaa                                                                                                                                                                                                            <  1:551266/251‑1 (MQ=255)
                                             gagtGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTGTACGCGTATGACAGGATACGGCGGACGGTTGTGGCGCTCGTCTTCGGTGATCGCATTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTACACGTTATCAGCAAGCGtt                                                                                                                                                                                                       <  2:485196/248‑1 (MQ=255)
                                                 gTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCTGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGTTGATCGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGTCGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTacac                                                                                                                                                                                                   <  1:187313/251‑1 (MQ=255)
                                                 gTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGTCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTaca                                                                                                                                                                                                    <  2:450927/250‑1 (MQ=255)
                                                           aaaTGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACCCTCAGCGCAtt                                                                                                                                                                                         <  1:238609/251‑1 (MQ=255)
                                                               ggACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGc                                                                                                                                                                                     <  2:67551/251‑1 (MQ=255)
                                                                gACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGca                                                                                                                                                                                    <  2:110231/251‑2 (MQ=255)
                                                                     aCAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATa                                                                                                                                                                               <  1:40127/251‑1 (MQ=255)
                                                                           gggCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAACGCTGCACGTTATCAGCACGCGTTACACTCAGCGCATTGAGCGACATaatctg                                                                                                                                                                          >  1:317519/1‑250 (MQ=255)
                                                                                    cGGTGCTAAATCTCGTCAGCGCTGGCGGTTTTACGCGTATGACAGGATACGGAGTACGGTTGTGGCGCACGTCGTCGGTGAACGCACTCTGGCCTCACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACTGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATATTCTTAATCTgaga                                                                                                                                                                 <  2:554335/250‑1 (MQ=255)
                                                                                     ggTGCTAAATCCCGTCAGCGCTGGATGCTTTACGCGAATGCCACGATCCGGAGGACGGTTGTGGCGCACGTCTTCGGCGAACGCACTCTGGCCACACTGGCGCGCCTTCTGAGCCTGCTGACGGCCTTTGAGGTCGTGGTATAGATGACGGCTGGATCGGCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGCTATCAGCAAGCGTTACACTACGCGCACTGGGGGAGATAAACTGAATCTGGGac                                                                                                                                                                >  2:308344/1‑250 (MQ=255)
                                                                                     ggTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                                                                                                       >  2:236885/1‑227 (MQ=255)
                                                                                     ggTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGa                                                                                                                                                                                       <  1:236885/227‑1 (MQ=255)
                                                                                      gTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGaca                                                                                                                                                               >  1:275649/1‑250 (MQ=255)
                                                                                                 cGTCAGCGCTGGCTGTTTTACGCGAATGACAGGATACGGAGGACGGTTGTGGCACACGACTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCTGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCATGAAGGGAAAGTTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAACGACATAATCTGAATCTGAGACAACATCTAGCaa                                                                                                                                                    >  1:416930/1‑250 (MQ=255)
                                                                                                  gacAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGACTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTacac                                                                                                                                                                                                   <  1:158451/200‑1 (MQ=255)
                                                                                                        gctggctgTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTgggg                                                                                                                                            >  1:488257/1‑250 (MQ=255)
                                                                                                                 tgtACGCGTATGACAGGCTCCGGAAGCCGGTTGTTGCGCACGTATTCGGTGAACGAACTAGGGCGGCGCTGGGGCGTCGTATGAGCCTGCTGGCACCCTTTGACGTGGTGATATGGATGAAGGATGGCTGGCCGCTGTATGAATCACGCCGGAAGGGAAAGCTACACGTTATCAGCAAGCGTTAAACGCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGTACGGAAGtt                                                                                                                                    <  2:323766/248‑1 (MQ=255)
                                                                                                                      gcgTATGACAGGATACGGAGGACGGTTGTGGCGTACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTc                                                                                                                              <  1:188445/251‑1 (MQ=255)
                                                                                                                      gcgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTc                                                                                                                              <  1:14687/251‑1 (MQ=255)
                                                                                                                       cgTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCg                                                                                                                             >  1:85280/1‑251 (MQ=255)
                                                                                                                                ggATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCaaa                                                                                                                     >  1:269119/1‑250 (MQ=255)
                                                                                                                                ggATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCaaa                                                                                                                     >  1:489596/1‑250 (MQ=255)
                                                                                                                                    aCGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCgg                                                                                                                <  2:275649/251‑1 (MQ=255)
                                                                                                                                        aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCa                                                                                                                       >  2:58683/1‑240 (MQ=255)
                                                                                                                                        aggaCGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCa                                                                                                                       <  1:58683/240‑1 (MQ=255)
                                                                                                                                         gggcGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAg                                                                                                           <  2:568489/248‑1 (MQ=255)
                                                                                                                                                tgtgGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCCTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACAATCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATg                                                                                                    >  2:651124/1‑251 (MQ=255)
                                                                                                                                                    gcgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACa                                                                                                 >  2:367646/1‑250 (MQ=255)
                                                                                                                                                     cgcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACa                                                                                                 >  1:280181/1‑249 (MQ=255)
                                                                                                                                                     cgcACGTCGTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAg                                                                                               <  1:60712/251‑1 (MQ=255)
                                                                                                                                                      gcACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAgg                                                                                              <  1:544792/251‑1 (MQ=255)
                                                                                                                                                       cACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGt                                                                                             <  1:272977/251‑1 (MQ=255)
                                                                                                                                                              acGGTGAACGCACTCTGGCGACACTGGAGCGGCTTCTGAGCCTGCTGGCGGACTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATc                                                                                         <  2:651191/247‑1 (MQ=255)
                                                                                                                                                               cGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGc                                                                                     >  1:17297/1‑251 (MQ=255)
                                                                                                                                                                ggTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCa                                                                                    <  1:588648/251‑1 (MQ=255)
                                                                                                                                                                     aCGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATc                                                                               <  2:269119/251‑1 (MQ=255)
                                                                                                                                                                             tGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACAt                                                                        >  2:166291/1‑250 (MQ=255)
                                                                                                                                                                               gCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATaaa                                                                     <  1:348375/251‑1 (MQ=255)
                                                                                                                                                                               gCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATaa                                                                      >  2:60606/1‑250 (MQ=255)
                                                                                                                                                                                           ggggttctGGGCCTGCTGTCGGCCTTTGGGGTGGGGGTATAGGTGGGGGATGGCTGGGCGTTGGATGGATCACGCCTGAAGGGAAAGCTGGACGTTAGCAGCAAGCGTTACACTCAGGGCATTGAGCGGCATAATCTGAATCTGAGACAACATCTGGGAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACAGAAAACACTATCAGTa                                                         <  1:190445/247‑1 (MQ=255)
                                                                                                                                                                                            gTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTaa                                                        <  2:489596/251‑1 (MQ=255)
                                                                                                                                                                                              cttctGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGt                                                      >  2:442998/1‑251 (MQ=255)
                                                                                                                                                                                                      ccTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTc                                               >  2:65147/1‑250 (MQ=255)
                                                                                                                                                                                                         gctgTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCAtt                                            <  2:492480/250‑1 (MQ=255)
                                                                                                                                                                                                         gctgTCGGCCTGTGAGGTCGTGGTATAGATGACGGGTGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCTCGTGATCAGCGAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAACCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAAACGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAGCATAAAACACTATCAGTAAGGAGGAGTCAtt                                            <  2:197057/250‑1 (MQ=255)
                                                                                                                                                                                                              cGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGGAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATGGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGTGACGGAAGTGACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGGCATCGGGCATTATCTGATCATAAATCTCTATCAGTAAGTTGGAGTCATTACCGAc                                      <  2:378287/251‑1 (MQ=255)
                                                                                                                                                                                                               ggCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                                                                                                                                                                    <  1:104313/124‑1 (MQ=255)
                                                                                                                                                                                                               ggCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATaatctgaatctg                                                                                                                                                                    >  2:104313/1‑124 (MQ=255)
                                                                                                                                                                                                                        ggTCGTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTAtt                            <  1:534607/251‑1 (MQ=255)
                                                                                                                                                                                                                            gTGGTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGTAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCAt                        <  1:545416/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATAc                      <  1:36043/251‑1 (MQ=255)
                                                                                                                                                                                                                              ggTATAGATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATAc                      >  2:560361/1‑251 (MQ=255)
                                                                                                                                                                                                                                    gATGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATtgtg                >  1:137844/1‑251 (MQ=255)
                                                                                                                                                                                                                                     aTGACGGATGGCTGGCCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTgg               <  1:65147/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 tggcCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCACCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCtt    >  2:572549/1‑250 (MQ=255)
                                                                                                                                                                                                                                                  ggcCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCttat  >  1:515529/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  ggcCGTTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCtt    <  2:280181/249‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                     
CAGGCCGCAGTCGGTAACCTCGCGCATACAACCGGGCAGTGATGTGATTGTCTGCGCTGAAATGGACGAACAGTGGGGCTACGTCGGTGCTAAATCACGTCAGCGCTGGCTGTTTTACGCGTATGACAGGATACGGAGGACGGTTGTGGCGCACGTCTTCGGTGAACGCACTCTGGCCACACTGGAGCGTCTTCTGAGCCTGCTGTCGGCCTTTGAGGTCGTGGTATAGATGACGGATGGCTGGCCGCTGTATGAATCACGCCTGAAGGGAAAGCTGCACGTTATCAGCAAGCGTTACACTCAGCGCATTGAGCGACATAATCTGAATCTGAGACAACATCTGGCAAGGCTGGGACGGAAGTTACTGTCGTTCTCAAAATCGGTGGAGTTGCATGACAAGGTCATCGGGCATTATCTGAACATAAAACACTATCAGTAAGTTGGAGTCATTACCGACCATGTTTATTTCATACATTGTGGGTATTGTTCTTAT  >  NC_000913/4518786‑4519278

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: