Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,497 A→C 6.8% intergenic (‑180/+35) 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,4970AC 6.8% 341.1 / 14.7 103intergenic (‑180/+35)yijO/eptCAraC family putative transcriptional activator/LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base A (44/52);  new base C (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.58e-01

Only 100 of 107 total aligned reads displayed.

CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCA  >  NC_000913/4148354‑4148644
                                                                                                                                               |                                                                                                                                                   
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)
         ttaAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCTGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATGCGACGCTTGACGCGTCTTATCCgg                                                                                                                                        <  2:389931/148‑1 (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)
             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)
                                     cTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGAt                                                                                                           <  4:764295/149‑1 (MQ=255)
                                            gAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCa                                                                                                     <  1:183155/148‑1 (MQ=255)
                                                tCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTg                                                                                                 >  2:905754/1‑148 (MQ=255)
                                                tCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTg                                                                                                 >  1:698627/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)
                                                               cGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAgt                                                                                 <  3:746900/149‑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:679367/1‑149 (MQ=255)
                                                                     ttGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCg                                                                           >  4:430553/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)
                                                                                               ccTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTtat                                                  <  1:740056/148‑1 (MQ=255)
                                                                                                cTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTtata                                                 >  1:181805/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                                >  2:390607/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                                <  2:268816/148‑1 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                                >  1:123067/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                                >  2:788562/1‑148 (MQ=255)
                                                                                                  tATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACgg                                              <  3:679367/149‑1 (MQ=255)
                                                                                                    tCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATc                                                                                              <  3:382309/99‑1 (MQ=255)
                                                                                                    tCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATc                                                                                              >  4:382309/1‑99 (MQ=255)
                                                                                                     ccAGCCTACCGCCGCGATCAGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTGTTCTTATACgg                                              <  3:541158/146‑1 (MQ=255)
                                                                                                        gCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTAc                                         >  1:486468/1‑148 (MQ=255)
                                                                                                          cTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCaa                                      <  4:138211/149‑1 (MQ=255)
                                                                                                            aCCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                                     <  2:683139/148‑1 (MQ=255)
                                                                                                            aCCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                                     <  2:932799/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                 <  2:231951/148‑1 (MQ=255)
                                                                                                                                ggATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGa                 >  1:179603/1‑148 (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:225930/148‑1 (MQ=255)
                                                                                                                                       aCGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCgg          <  1:826991/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)
                                                                                                                                              aaGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcacca  <  4:514451/149‑1 (MQ=255)
                                                                                                                                              aaGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTcacca  >  4:90087/1‑149 (MQ=255)
                                                                                                                                               aGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATc                                                                  >  2:146138/1‑84 (MQ=255)
                                                                                                                                               aGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATc                                                                  <  1:146138/84‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                   
CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGATTCACCA  >  NC_000913/4148354‑4148644

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: