Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,495 C→G 6.8% intergenic (‑178/+37) yijO ← / ← eptC AraC family putative transcriptional activator/LPS heptose I phosphoethanolamine transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,148,4950CG 6.8% 357.3 / 14.2 103intergenic (‑178/+37)yijO/eptCAraC family putative transcriptional activator/LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base C (44/52);  new base G (2/5);  total (46/57)
Fisher's exact test for biased strand distribution p-value = 4.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-01

Only 100 of 103 total aligned reads displayed.

CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCAC  >  NC_000913/4148354‑4148642
                                                                                                                                             |                                                                                                                                                   
cTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGt                                                                                                                                               <  2:306065/148‑1 (MQ=255)
cTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGt                                                                                                                                               <  1:386236/148‑1 (MQ=255)
cTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTc                                                                                                                                              <  4:306424/149‑1 (MQ=255)
     attattAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTcttatc                                                                                                                                         <  3:748540/149‑6 (MQ=255)
         ttaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGc                                                                                                                                     >  4:906189/1‑149 (MQ=255)
          taAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATCCGGc                                                                                                                                     <  2:645810/148‑1 (MQ=255)
           aaaaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATcc                                                                                                                                        >  3:436784/1‑144 (MQ=39)
           aaaaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATcc                                                                                                                                        <  4:436784/144‑1 (MQ=39)
             aaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTCTCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCTTGACGCGTct                                                                                                                                             <  2:828852/137‑2 (MQ=255)
             aaaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCTTGACGCGTct                                                                                                                                             >  1:828852/1‑136 (MQ=255)
              aaaGGCCACCTGAATAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGGAGGCCGGATAATACGCGTCAAGCGTCGCATCCGGCAAAt                                                                                                                                 >  2:631236/1‑148 (MQ=255)
              aaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGGGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAAt                                                                                                                                 >  1:285777/1‑148 (MQ=255)
              aaaGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATa                                                                                                                                <  3:906189/149‑1 (MQ=255)
                aGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAg                                                                                                                               >  2:322851/1‑148 (MQ=255)
                    cACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCt                                                                                                                          <  3:137609/149‑1 (MQ=255)
                        tGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCt                                                                                                                      >  3:485947/1‑149 (MQ=255)
                         gAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTg                                                                                                                     >  4:362507/1‑149 (MQ=255)
                          aaaaGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGa                                                                                                                    <  4:570944/149‑1 (MQ=255)
                             aGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATaa                                                                                                                 <  4:51293/149‑1 (MQ=255)
                             aGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATaa                                                                                                                 <  3:294737/149‑1 (MQ=255)
                             aGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATa                                                                                                                  >  1:231951/1‑148 (MQ=255)
                                gtgGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAAtt                                                                                                               <  1:300182/148‑1 (MQ=255)
                                 tgGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTa                                                                                                              >  2:224811/1‑148 (MQ=255)
                                   gCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTg                                                                                                           >  4:156381/1‑149 (MQ=255)
                                     cTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAt                                                                                                         <  4:764295/149‑1 (MQ=255)
                                       gAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAtt                                                                                                        >  2:225930/1‑148 (MQ=255)
                                            gAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCa                                                                                                   <  1:183155/148‑1 (MQ=255)
                                                tCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTg                                                                                               >  1:698627/1‑148 (MQ=255)
                                                tCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTg                                                                                               >  2:905754/1‑148 (MQ=255)
                                                  gATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATc                                                                                            <  3:136459/149‑1 (MQ=255)
                                                    tGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGc                                                                                          >  4:784806/1‑149 (MQ=255)
                                                        aaCATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCAt                                                                                       <  2:477237/148‑1 (MQ=255)
                                                         aCATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAg                                                                                     <  4:929156/149‑1 (MQ=255)
                                                           tttGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCTTCCCATCCTGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAgga                                                                                   >  4:206000/2‑148 (MQ=255)
                                                             tGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCa                                                                                 <  3:795413/149‑1 (MQ=255)
                                                             tGCGGCAGATGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGc                                                                                  <  2:834146/148‑1 (MQ=255)
                                                                ggCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtg                                                                              >  4:679374/1‑149 (MQ=255)
                                                                 gCAGTTGTATTGCCGGATGTTGCGTAAACGCCTTATCCAGCCTACCGCTGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtgt                                                                             <  3:784806/149‑1 (MQ=255)
                                                                 gCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtgt                                                                             <  4:494930/149‑1 (MQ=255)
                                                                    gTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTc                                                                          <  3:694525/149‑1 (MQ=255)
                                                                    gTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTc                                                                          <  4:798115/149‑1 (MQ=255)
                                                                     ttGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCg                                                                         >  4:430553/1‑149 (MQ=255)
                                                                     ttGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCg                                                                         >  4:679367/1‑149 (MQ=255)
                                                                       gTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTa                                                                       <  3:362507/149‑1 (MQ=255)
                                                                          ttGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATc                                                                    >  3:567999/1‑149 (MQ=255)
                                                                          ttGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATc                                                                    >  3:138211/1‑149 (MQ=255)
                                                                            gCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCg                                                                   <  1:224811/148‑1 (MQ=255)
                                                                            gCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGa                                                                  <  3:701854/149‑1 (MQ=255)
                                                                              cGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGAGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGAt                                                                 <  1:355324/148‑1 (MQ=255)
                                                                                gATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCa                                                               >  1:268816/1‑148 (MQ=255)
                                                                                 aTGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGt                                                             <  4:466676/149‑1 (MQ=255)
                                                                                  tgtgGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTg                                                            <  3:696077/149‑1 (MQ=255)
                                                                                      gCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGtt                                                        >  3:139654/1‑149 (MQ=255)
                                                                                      gCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGAAAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGtt                                                        >  3:545144/1‑149 (MQ=255)
                                                                                       cGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGtt                                                        >  2:932439/1‑148 (MQ=255)
                                                                                          aaaCGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGa                                                                  >  2:23381/1‑135 (MQ=255)
                                                                                          aaaCGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGa                                                                  <  1:23381/135‑1 (MQ=255)
                                                                                             cGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCtt                                                  <  1:3799/148‑1 (MQ=255)
                                                                                               ccTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTtat                                                <  1:740056/148‑1 (MQ=255)
                                                                                                cTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTtata                                               >  1:181805/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                              >  1:123067/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                              >  2:390607/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                              <  2:268816/148‑1 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                              >  2:788562/1‑148 (MQ=255)
                                                                                                  tATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACgg                                            <  3:679367/149‑1 (MQ=255)
                                                                                                    tCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATc                                                                                            >  4:382309/1‑99 (MQ=255)
                                                                                                    tCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATc                                                                                            <  3:382309/99‑1 (MQ=255)
                                                                                                     ccAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGCTTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACgg                                            >  4:541158/1‑146 (MQ=255)
                                                                                                     ccAGCCTACCGCCGCGATCAGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTGTTCTTATACgg                                            <  3:541158/146‑1 (MQ=255)
                                                                                                      cAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTa                                        >  3:514451/1‑149 (MQ=255)
                                                                                                        gCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTAc                                       >  1:486468/1‑148 (MQ=255)
                                                                                                          cTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCaa                                    <  4:138211/149‑1 (MQ=255)
                                                                                                            aCCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                                   <  2:932799/148‑1 (MQ=255)
                                                                                                            aCCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                                   <  2:683139/148‑1 (MQ=255)
                                                                                                                   cgATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCgg                            <  1:905754/148‑1 (MQ=255)
                                                                                                                     aTCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTa                         <  4:602626/149‑1 (MQ=255)
                                                                                                                     aTCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTa                         <  4:567999/149‑1 (MQ=255)
                                                                                                                      tCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCTTCCGGCAAATAGCGCCCAGCTGATAATTACTAATTACCCACCTGATCGCCATACGGCAGTGTGTCGTAATCGATCAGCGCGTTTTTCTTATACGGGTTACAAATCCAGCGGGTAg                        >  4:524770/1‑149 (MQ=255)
                                                                                                                         gTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtgt                                                                             >  1:610807/1‑93 (MQ=255)
                                                                                                                         gTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgtgt                                                                             <  2:610807/93‑1 (MQ=255)
                                                                                                                           aGGGCGGATAAGACGCGTCAATCGTCGCATACGGCAAATAGCGCCTGGCTGATAATTACTTATTACACATCTGATCGCCATAGGACAGTGTTTCGTAATCGTTTAGTGTGTTTTTCTTATACGGGTTACGAATCCATTGGGTAGTTTCt                   <  4:662808/149‑1 (MQ=255)
                                                                                                                           aGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTc                    <  2:181805/148‑1 (MQ=255)
                                                                                                                            ggCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCtt                  <  4:808811/149‑1 (MQ=255)
                                                                                                                            ggCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCt                   <  2:909246/148‑1 (MQ=255)
                                                                                                                            ggCCAGATAAGACGCGTCAAGCGTCGCATCAGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCt                   <  1:631236/148‑1 (MQ=255)
                                                                                                                              ccGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCttt                 >  1:291917/1‑148 (MQ=255)
                                                                                                                                ggATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGaa              <  4:545144/149‑1 (MQ=255)
                                                                                                                                ggATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGa               >  1:179603/1‑148 (MQ=255)
                                                                                                                                ggATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGa               <  2:231951/148‑1 (MQ=255)
                                                                                                                                  aTAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACt            >  3:570928/1‑149 (MQ=255)
                                                                                                                                  aTAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACt            >  3:373251/1‑149 (MQ=255)
                                                                                                                                   tAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACt            >  1:775688/1‑148 (MQ=255)
                                                                                                                                    aaGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTg           >  1:194745/1‑148 (MQ=255)
                                                                                                                                    aaGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGc          <  4:310537/149‑1 (MQ=255)
                                                                                                                                       aCGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCgg        <  1:826991/148‑1 (MQ=255)
                                                                                                                                       aCGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCgg        <  1:225930/148‑1 (MQ=255)
                                                                                                                                       aCGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGa       >  4:16163/1‑149 (MQ=255)
                                                                                                                                        cgcgTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGa       >  1:540089/1‑148 (MQ=255)
                                                                                                                                        cgcgTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGa       >  2:693511/1‑148 (MQ=255)
                                                                                                                                            tCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcac  >  4:453468/1‑149 (MQ=255)
                                                                                                                                             |                                                                                                                                                   
CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCAC  >  NC_000913/4148354‑4148642

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: