Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,855 T→A 7.3% intergenic (‑253/+113) yhiI ← / ← yhiJ putative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,630,8550TA 7.3% 320.8 / 10.8 94intergenic (‑253/+113)yhiI/yhiJputative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base T (39/48);  new base A (3/4);  total (42/52)
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

Only 100 of 101 total aligned reads displayed.

CAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGG  >  NC_000913/3630712‑3631003
                                                                                                                                               |                                                                                                                                                    
cAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTg                                                                                                                                                  >  1:910239/1‑148 (MQ=255)
cAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGt                                                                                                                                                 <  4:340669/149‑1 (MQ=255)
 aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACACATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                                <  4:815277/149‑1 (MQ=255)
    gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                                <  4:798615/146‑1 (MQ=255)
    gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                                >  3:798615/1‑146 (MQ=255)
      gcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATa                                                                                                                                           <  4:838642/149‑1 (MQ=255)
      gcAGAGCAGTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGAGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAATAACGCATACTCATAAGTGTGATCTATAAAACGGGCTGACAAACCGGTATCCATCCTGTTGGATa                                                                                                                                           <  4:618234/149‑1 (MQ=255)
         gagCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCg                                                                                                                                        >  3:586257/1‑149 (MQ=255)
          agCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGt                                                                                                                                       <  4:468013/149‑1 (MQ=255)
             aTTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttt                                                                                                                                     <  2:66259/148‑1 (MQ=255)
               tGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATac                                                                                                                                          >  2:105927/1‑140 (MQ=255)
                gAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAATCTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttttgtc                                                                                                                                 <  3:528528/149‑5 (MQ=255)
                 aGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACg                                                                                                                                 >  2:460721/1‑148 (MQ=255)
                 aGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACg                                                                                                                                 <  2:722039/148‑1 (MQ=255)
                  gcgcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGc                                                                                                                                >  1:410907/1‑148 (MQ=255)
                   cgcCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCt                                                                                                                              <  4:730197/149‑1 (MQ=255)
                      cAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGtt                                                                                                                           <  3:623751/149‑1 (MQ=255)
                            ggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATAATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTg                                                                                                                     >  4:855638/1‑149 (MQ=255)
                              aaGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCa                                                                                                                   <  4:888851/149‑1 (MQ=255)
                               aGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAg                                                                                                                  <  3:855638/149‑1 (MQ=255)
                                   ttAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAAc                                                                                                              >  4:337144/1‑149 (MQ=255)
                                     aGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACag                                                                                                            <  4:811415/149‑1 (MQ=255)
                                      gTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACaga                                                                                                           <  4:676948/149‑1 (MQ=255)
                                      gTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACag                                                                                                            <  1:460721/148‑1 (MQ=255)
                                        gAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACagag                                                                                                          <  1:547773/148‑1 (MQ=255)
                                        gAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACagag                                                                                                          >  1:37357/1‑148 (MQ=255)
                                           aaCGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCAt                                                                                                       >  1:547462/1‑148 (MQ=255)
                                            aCGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTc                                                                                                     <  4:728686/149‑1 (MQ=255)
                                                 tagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtttt                                                                                                >  3:22176/1‑149 (MQ=255)
                                                 tagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCttt                                                                                                 <  2:135270/148‑1 (MQ=255)
                                                 tagtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGGTTCACGCCTGTTCAGTTGCAGGAACAGAGCTTTCCtttt                                                                                                <  4:226602/149‑1 (MQ=255)
                                                 tagtaATGTTGTCATCATCAGCATACCTGGCCGTAAATGAAACACGCATACTCATCAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtttt                                                                                                <  4:58976/149‑1 (MQ=255)
                                                  agtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtttt                                                                                                <  2:525671/148‑1 (MQ=255)
                                                        gttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATa                                                                                                                                           <  3:670284/99‑1 (MQ=255)
                                                        gttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATa                                                                                                                                           >  4:670284/1‑99 (MQ=255)
                                                           gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACaga                                                                                                           >  4:161335/1‑128 (MQ=255)
                                                           gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACaga                                                                                                           <  3:161335/128‑1 (MQ=255)
                                                           gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTa                                                                                       <  1:361156/148‑1 (MQ=255)
                                                           gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTa                                                                                       <  2:910239/148‑1 (MQ=255)
                                                               tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAAc                                                                                  <  3:703662/149‑1 (MQ=255)
                                                               tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAAc                                                                                  >  4:302384/1‑149 (MQ=255)
                                                                catcaGCATTCCTGGCCGTAAATGAAACACGCATACTAATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACa                                                                                 <  3:586969/149‑1 (MQ=255)
                                                                  tcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                >  2:868352/1‑148 (MQ=255)
                                                                   caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGTATTCCTTTCTACGTTTTATCAACAGtc                                                                              >  3:488784/1‑149 (MQ=255)
                                                                   caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttttgtc                                                                                                                                 >  1:632606/1‑94 (MQ=25)
                                                                   caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttttgtc                                                                                                                                 <  2:632606/98‑5 (MQ=25)
                                                                     gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtctc                                                                            >  3:287798/1‑149 (MQ=255)
                                                                     gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtctc                                                                            >  4:361973/1‑149 (MQ=255)
                                                                      cATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGATCATTCCTTTTTACGTTTTATCAACAGtctc                                                                            >  1:592115/1‑148 (MQ=255)
                                                                        ttCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCaa                                                                          <  2:735221/148‑1 (MQ=255)
                                                                        ttCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCaa                                                                          <  1:811989/148‑1 (MQ=255)
                                                                           cTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                        <  4:446793/131‑1 (MQ=255)
                                                                           cTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                        >  3:446793/1‑131 (MQ=255)
                                                                           cTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAg                                                                      >  4:724596/1‑149 (MQ=255)
                                                                              gCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATa                                                                   <  3:296290/149‑1 (MQ=255)
                                                                                 gTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTa                                                                >  3:257222/1‑149 (MQ=255)
                                                                                  tAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTaa                                                               <  4:22176/149‑1 (MQ=255)
                                                                                     aTGAAACACACATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                             <  1:490250/132‑1 (MQ=255)
                                                                                     aTGAAACACACATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                             >  2:490250/1‑132 (MQ=255)
                                                                                       gAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACAtt                                                           >  1:520937/1‑148 (MQ=255)
                                                                                        aaaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTa                                                          <  1:394583/148‑1 (MQ=255)
                                                                                          acacGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAat                                                        <  2:37357/148‑1 (MQ=255)
                                                                                          acacGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAat                                                        <  1:237194/148‑1 (MQ=255)
                                                                                          acacGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAat                                                        >  1:108093/1‑148 (MQ=255)
                                                                                            acGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAatat                                                      <  2:592115/148‑1 (MQ=255)
                                                                                            acGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATAtt                                                     <  4:287798/149‑1 (MQ=255)
                                                                                             cGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATAtt                                                     >  2:894563/1‑148 (MQ=255)
                                                                                              gCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                   <  4:257222/149‑1 (MQ=255)
                                                                                               cATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                   <  1:97446/148‑1 (MQ=255)
                                                                                                      aTAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCt                                           >  4:127680/1‑149 (MQ=255)
                                                                                                          gtgtGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATgg                                       >  3:268563/1‑149 (MQ=255)
                                                                                                             tgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGtt                                     >  1:314601/1‑148 (MQ=255)
                                                                                                                  tataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGc                               >  3:653453/1‑149 (MQ=255)
                                                                                                                   ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTa                                                                                       >  1:368052/1‑92 (MQ=255)
                                                                                                                   ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTa                                                                                       <  2:368052/92‑1 (MQ=255)
                                                                                                                    taAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCaa                                                                                   >  4:806382/1‑95 (MQ=255)
                                                                                                                    taAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCaa                                                                                   <  3:806382/95‑1 (MQ=255)
                                                                                                                    taAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                             >  3:821628/1‑149 (MQ=255)
                                                                                                                     aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                             >  1:422061/1‑148 (MQ=255)
                                                                                                                     aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                             <  2:480361/148‑1 (MQ=255)
                                                                                                                        aCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTCACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                   >  2:675759/1‑123 (MQ=255)
                                                                                                                        aCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTCACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                                   <  1:675759/123‑1 (MQ=255)
                                                                                                                           ggCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACAt                       >  1:421827/1‑148 (MQ=255)
                                                                                                                             cTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtgt                    <  4:239248/149‑1 (MQ=255)
                                                                                                                             cTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtg                     >  1:641277/1‑148 (MQ=255)
                                                                                                                                  aaaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGcc               >  4:656989/1‑149 (MQ=255)
                                                                                                                                   aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGcc               >  1:507687/1‑148 (MQ=255)
                                                                                                                                   aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGcc               <  2:520937/148‑1 (MQ=255)
                                                                                                                                   aaaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCt              >  4:748907/1‑149 (MQ=255)
                                                                                                                                       ggTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCt           <  2:314601/148‑1 (MQ=255)
                                                                                                                                        gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTaa         <  3:127680/149‑1 (MQ=255)
                                                                                                                                        gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTaa         <  3:95008/149‑1 (MQ=255)
                                                                                                                                         tATCCATCCTGTTGGGTAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTaa         <  2:832303/148‑1 (MQ=255)
                                                                                                                                             catccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCaaa    >  3:568165/1‑149 (MQ=255)
                                                                                                                                             catccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCaa     <  2:410907/148‑1 (MQ=255)
                                                                                                                                              atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTAt                                         >  2:391013/1‑111 (MQ=255)
                                                                                                                                              atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTAt                                         <  1:391013/111‑1 (MQ=255)
                                                                                                                                               tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg  <  3:748907/149‑1 (MQ=255)
                                                                                                                                               tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAg   <  2:641277/148‑1 (MQ=255)
                                                                                                                                               tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAg   >  2:18155/1‑148 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
CAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGG  >  NC_000913/3630712‑3631003

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: