Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 18,681 A→T 5.9% intergenic (+26/‑34) nhaA → / → nhaR sodium‑proton antiporter/transcriptional activator of nhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091318,6810AT 5.9% 309.6 / 5.9 96intergenic (+26/‑34)nhaA/nhaRsodium‑proton antiporter/transcriptional activator of nhaA
Reads supporting (aligned to +/- strand):  ref base A (50/40);  new base T (3/3);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTA  >  NC_000913/18534‑18826
                                                                                                                                                   |                                                                                                                                                 
ccTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATaaa                                                                                                                                                  >  4:280676/1‑149 (MQ=255)
ccTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATaaa                                                                                                                                                  <  7:310341/149‑1 (MQ=255)
      gTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCt                                                                                                                                             <  6:395892/148‑1 (MQ=255)
       tAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTcc                                                                                                                                           <  8:174256/149‑1 (MQ=255)
          cGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCGTTTATGGCTccccg                                                                                                                                        <  7:77459/149‑3 (MQ=255)
            tAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTAGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCtttt                                                                                                                                      >  2:358947/1‑149 (MQ=255)
              gATCCTGAACTGTTTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACGGGGGAGCCATAAACGGCTCCCTTTTc                                                                                                                                     <  8:388641/148‑1 (MQ=255)
              gATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCa                                                                                                                                    <  3:280676/149‑1 (MQ=255)
              gATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCa                                                                                                                                    <  3:323337/149‑1 (MQ=255)
               aTCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCAt                                                                                                                                   <  6:117586/149‑1 (MQ=255)
                    gAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTATGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTa                                                                                                                              >  3:323707/1‑149 (MQ=255)
                     aaCTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTAt                                                                                                                             >  1:137118/1‑149 (MQ=255)
                         gATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCtt                                                                                                                                        >  1:293174/1‑134 (MQ=255)
                         gATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCtt                                                                                                                                        <  2:293174/134‑1 (MQ=255)
                          aTTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGTATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                         >  6:437101/1‑148 (MQ=255)
                          aTTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAg                                                                                                                          <  2:271756/147‑1 (MQ=255)
                          aTTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTATTCGGCGGTAATTGGATACATCTTGTTACGCCTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCAAAAAAGGCTCCCTTAACATTGTTATCAgg                                                                                                                         >  5:24203/1‑148 (MQ=255)
                              aCTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGga                                                                                                                       >  2:215810/1‑146 (MQ=255)
                                tGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGgagaga                                                                                                                   >  8:2858/1‑148 (MQ=255)
                                tGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGgagaga                                                                                                                   >  1:112212/1‑148 (MQ=255)
                                  ggCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGaaa                                                                                                                 <  7:336333/148‑1 (MQ=255)
                                  ggCGAAACTCGGTAGCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCGCCCTTTTCATTGTTATCAGGGAGAGAAAt                                                                                                                <  6:373229/149‑1 (MQ=255)
                                      aaaCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGc                                                                                                            <  7:2858/149‑1 (MQ=255)
                                      aaaCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGc                                                                                                            <  4:12007/149‑1 (MQ=255)
                                        aCTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCAATCAGTTTGACATGACGGTTTACCGGGTAGCCATATACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCAt                                                                                                          <  1:358947/149‑1 (MQ=255)
                                         cTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATg                                                                                                         >  8:97407/1‑149 (MQ=255)
                                            ggTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGtc                                                                                                       <  1:61881/148‑1 (MQ=255)
                                                tCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGaa                                                                                                                  <  8:53432/133‑1 (MQ=255)
                                                tCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGaa                                                                                                                  >  7:53432/1‑133 (MQ=255)
                                                    ggTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAAt                                                                                                                >  7:243810/1‑131 (MQ=255)
                                                    ggTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAAt                                                                                                                <  8:243810/131‑1 (MQ=255)
                                                     gTCGGTTCTAGCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAAt                                                                                             <  2:32647/149‑1 (MQ=255)
                                                      tCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAAtt                                                                                            <  1:215810/149‑1 (MQ=255)
                                                             tATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATaaa                                                                                                                                                  >  8:404258/1‑88 (MQ=255)
                                                             tATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATaaa                                                                                                                                                  <  7:404258/88‑1 (MQ=255)
                                                               tctcTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCata                                                                                                  >  6:75734/1‑134 (MQ=255)
                                                               tctcTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTGGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCata                                                                                                  <  5:75734/134‑1 (MQ=255)
                                                                 tctTAGGAGGTAATTGGATACAGCTGGTTACGAGTTCGTTTGCTTCCATCAGTTAGACATGACGGTATACCGGGTAGACATAAAc                                                                                                                                                 >  8:260341/1‑85 (MQ=255)
                                                                   ttCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACtt                                                                                  <  5:437101/146‑1 (MQ=255)
                                                                   ttCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTa                                                                               >  7:311645/1‑149 (MQ=255)
                                                                   ttCGGCGGTAATTGGATACAGCTGGTTACCCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGCGGAGCAATAAACGGCTCCCTTTACATTGTAATCAGGGAGAGAAATGAGAATGTCTCATATCAATAAAAACCACTTGTa                                                                               >  4:79574/1‑149 (MQ=255)
                                                                     cggcggTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTAt                                                                              <  8:311645/148‑1 (MQ=255)
                                                                     cggcggTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTAt                                                                              >  2:253800/1‑148 (MQ=255)
                                                                           tAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCt                                                                       <  1:253800/149‑1 (MQ=255)
                                                                                 gATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCAt                                                                  <  2:138406/148‑1 (MQ=255)
                                                                                  aTACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGt                                                                >  5:360525/1‑149 (MQ=255)
                                                                                      aGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                         <  6:201672/88‑1 (MQ=255)
                                                                                      aGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                         >  5:201672/1‑88 (MQ=255)
                                                                                      aGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCta                                                             >  8:396347/1‑148 (MQ=255)
                                                                                      aGCTGGTTAAGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCtat                                                            <  6:24203/149‑1 (MQ=255)
                                                                                         tGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCtata                                                           >  6:75955/1‑147 (MQ=255)
                                                                                           gTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAaag                                                        >  2:425616/1‑148 (MQ=255)
                                                                                            ttACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTg                                                                                                                                 <  1:202360/74‑1 (MQ=255)
                                                                                            ttACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTg                                                                                                                                 >  2:202360/1‑74 (MQ=255)
                                                                                              aCGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAaagaag                                                     <  2:112212/148‑1 (MQ=255)
                                                                                              aCGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAgg                                                    >  5:227899/1‑149 (MQ=255)
                                                                                               cgcgTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAgg                                                    <  5:75955/148‑1 (MQ=255)
                                                                                                  gttcgttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAAc                                                                                       >  2:344441/1‑110 (MQ=255)
                                                                                                  gttcgttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTcc                                                >  7:224454/1‑149 (MQ=255)
                                                                                                  gttcgttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTcc                                                >  1:103955/1‑149 (MQ=255)
                                                                                                    tcgttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGt                                              <  2:137118/149‑1 (MQ=255)
                                                                                                      gttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                         >  4:72421/1‑72 (MQ=255)
                                                                                                      gttTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAgg                                                                                                                         <  3:72421/72‑1 (MQ=255)
                                                                                                       tttGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGt                                                                >  1:195003/1‑128 (MQ=255)
                                                                                                       tttGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGt                                                                <  2:195003/128‑1 (MQ=255)
                                                                                                             tCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTAc                                                                                          >  7:306011/1‑96 (MQ=255)
                                                                                                             tCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTAc                                                                                          <  8:306011/96‑1 (MQ=255)
                                                                                                                aTCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGc                                  >  5:137039/1‑149 (MQ=255)
                                                                                                                    gTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAg                              >  6:103434/1‑149 (MQ=255)
                                                                                                                     tttGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCg                                               <  3:70198/131‑1 (MQ=255)
                                                                                                                     tttGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAgg                             >  4:190554/1‑149 (MQ=255)
                                                                                                                     tttGACAGGACGGTTTACAGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCg                                               >  4:70198/1‑131 (MQ=255)
                                                                                                                        gACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTTTCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACGACTTGTTTTACTTGTGGCATGTCTATTAAGAAGGTTCCGTGGTTGGGTCAGCGGAGgggc                          >  6:425833/1‑149 (MQ=255)
                                                                                                                        gACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGgcgc                          >  6:167488/1‑149 (MQ=255)
                                                                                                                          cAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtt                        <  6:227899/149‑1 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGt                                                   <  4:10272/120‑1 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGt                                                   >  3:10272/1‑120 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCttta                      >  5:231039/1‑149 (MQ=255)
                                                                                                                            ggACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCttta                      >  6:296654/1‑149 (MQ=255)
                                                                                                                             gACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttat                     >  8:250828/1‑149 (MQ=255)
                                                                                                                             gACGGTTTACCGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttatt                    >  2:236353/1‑149 (MQ=255)
                                                                                                                              aCGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttatt                    <  4:323707/149‑1 (MQ=255)
                                                                                                                               cGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                   <  5:167488/149‑1 (MQ=255)
                                                                                                                               cGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                   <  5:103434/149‑1 (MQ=255)
                                                                                                                                ggTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattta                  <  7:97407/149‑1 (MQ=255)
                                                                                                                                ggTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                   >  2:374332/1‑148 (MQ=255)
                                                                                                                                  tttACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAAc                >  3:151377/1‑149 (MQ=255)
                                                                                                                                      ccGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca            >  6:176535/1‑149 (MQ=255)
                                                                                                                                       cGGGTAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCcac           <  5:425833/149‑1 (MQ=255)
                                                                                                                                       aagggAGCCGTTTATGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                   <  5:268445/139‑1 (MQ=255)
                                                                                                                                       aagggAGCCGTTTATGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCtttattt                   >  6:268445/3‑141 (MQ=255)
                                                                                                                                       aagggAGCCGTTTATGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca            >  4:314566/3‑148 (MQ=255)
                                                                                                                                       aagggAGCCGTTTATGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCca            >  3:373248/3‑148 (MQ=255)
                                                                                                                                         gggAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCcaca          >  6:199042/1‑148 (MQ=255)
                                                                                                                                              ccATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCAt    >  5:407226/1‑149 (MQ=255)
                                                                                                                                                aTAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTc              <  8:112275/137‑1 (MQ=255)
                                                                                                                                                aTAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTc              >  7:112275/1‑137 (MQ=255)
                                                                                                                                                aTAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTa  >  5:230304/1‑149 (MQ=255)
                                                                                                                                                aTAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTa  >  5:993/1‑149 (MQ=255)
                                                                                                                                                 tAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATc                                                                                                >  6:35505/1‑54 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
CCTTTGGTAGCGTAGATCCAGAACTGATTAACTGGGCGAAACTCGGTATCCTGGTCGGTTCTATCTCTTCGGCGGTAATTGGATACAGCTGGTTACGCGTTCGTTTGCGTCCATCAGTTTGACAGGACGGTTTACCGGGGAGCCATAAACGGCTCCCTTTTCATTGTTATCAGGGAGAGAAATGAGCATGTCTCATATCAATTACAACCACTTGTATTACTTCTGGCATGTCTATAAAGAAGGTTCCGTGGTTGGCGCAGCGGAGGCGCTTTATTTAACTCCACAAACCATTA  >  NC_000913/18534‑18826

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