Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,148,493 G→T 6.8% intergenic (‑176/+39) 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,4930GT 6.8% 342.5 / 8.6 100intergenic (‑176/+39)yijO/eptCAraC family putative transcriptional activator/LPS heptose I phosphoethanolamine transferase
Reads supporting (aligned to +/- strand):  ref base G (41/52);  new base T (2/5);  total (43/57)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 104 total aligned reads displayed.

CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGA  >  NC_000913/4148354‑4148637
                                                                                                                                           |                                                                                                                                                
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)
             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)
                                                               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)
                                                                    gTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGAcgcg                                                                                                                                                  >  3:167632/1‑72 (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)
                                                                            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                                                             <  1:23381/135‑1 (MQ=255)
                                                                                          aaaCGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGa                                                             >  2:23381/1‑135 (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:268816/148‑1 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                         >  2:788562/1‑148 (MQ=255)
                                                                                                 ttATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATAc                                         >  2:390607/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)
                                                                                                     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:683139/148‑1 (MQ=255)
                                                                                                            aCCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAAt                              <  2:932799/148‑1 (MQ=255)
                                                                                                                     aTCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTa                    <  4:567999/149‑1 (MQ=255)
                                                                                                                     aTCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTa                    <  4:602626/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:373251/1‑149 (MQ=255)
                                                                                                                                  aTAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACt       >  3:570928/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)
                                                                                                                                           |                                                                                                                                                
CTGGCATTATTAAAAAAGGCCACCTGAAAAGTGTGGCCTGAGGGGAGTTCGATGTAAACATTGCGGCAGTTGTATTGCCGGATGTGGCGTAAACGCCTTATCCAGCCTACCGCCGCGATCTGTAGGCCGGATAAGACGCGTCAAGCGTCGCATCCGGCAAATAGCGCCTGGCTGATAATTACTGATTACCCACCTGATCGCCATAGGGCAGTGTGTCGTAATCGATCAGTGCGTTTTTCTTATACGGGTTACCAATCCAGCGGGTAGTTTCTTTGAACTGCGGA  >  NC_000913/4148354‑4148637

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: