Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,848 G→C 7.4% intergenic (‑246/+120) 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,8480GC 7.4% 319.0 / 14.2 95intergenic (‑246/+120)yhiI/yhiJputative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base G (41/47);  new base C (3/4);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.86e-01

CAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAA  >  NC_000913/3630702‑3630996
                                                                                                                                                  |                                                                                                                                                    
cAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTa                                                                                                                                                    <  3:849173/149‑1 (MQ=255)
  aaTACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTat                                                                                                                                                   >  1:480361/1‑148 (MQ=255)
    tACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTatcc                                                                                                                                                 <  1:308797/148‑1 (MQ=255)
          cAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTg                                                                                                                                           >  1:910239/1‑148 (MQ=255)
          cAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGt                                                                                                                                          <  4:340669/149‑1 (MQ=255)
           aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACACATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                         <  4:815277/149‑1 (MQ=255)
              gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                         >  3:798615/1‑146 (MQ=255)
              gcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGtt                                                                                                                                         <  4:798615/146‑1 (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)
                                                   aaaaCGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCa                                                                                                 >  4:95008/1‑149 (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)
                                                                 tgttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGt                                                                                                                                                     >  2:308618/1‑83 (MQ=255)
                                                                 tgttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGt                                                                                                                                                     <  1:308618/83‑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                                                                                <  2:910239/148‑1 (MQ=255)
                                                                     gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTa                                                                                <  1:361156/148‑1 (MQ=255)
                                                                         tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAAc                                                                           >  4:302384/1‑149 (MQ=255)
                                                                         tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAAc                                                                           <  3:703662/149‑1 (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                                                                                 >  3:446793/1‑131 (MQ=255)
                                                                                     cTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                 <  4:446793/131‑1 (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                                                                            <  3:806382/95‑1 (MQ=255)
                                                                                                                              taAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCaa                                                                            >  4:806382/1‑95 (MQ=255)
                                                                                                                              taAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                      >  3:821628/1‑149 (MQ=255)
                                                                                                                               aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                      <  2:480361/148‑1 (MQ=255)
                                                                                                                               aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                      >  1:422061/1‑148 (MQ=255)
                                                                                                                                  aCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTCACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                            <  1:675759/123‑1 (MQ=255)
                                                                                                                                  aCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTCACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCa                                            >  2:675759/1‑123 (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:95008/149‑1 (MQ=255)
                                                                                                                                                  gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTaa  <  3:127680/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CAAATACGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAA  >  NC_000913/3630702‑3630996

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: