Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,857 C→A 5.2% intergenic (‑255/+111) 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,8570CA 5.2% 455.6 / 8.7 146intergenic (‑255/+111)yhiI/yhiJputative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base C (79/59);  new base A (2/6);  total (81/65)
Fisher's exact test for biased strand distribution p-value = 1.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 158 total aligned reads displayed.

CGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAT  >  NC_000913/3630708‑3631006
                                                                                                                                                     |                                                                                                                                                     
cGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTatccatcc                                                                                                                                                       <  2:1233438/150‑1 (MQ=255)
cGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTatccatcc                                                                                                                                                       <  1:1275404/150‑1 (MQ=255)
  aCCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTg                                                                                                                                                     <  1:544129/150‑1 (MQ=255)
  aCCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTg                                                                                                                                                     >  1:1438796/1‑150 (MQ=255)
     aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAATCACGCATACTCATAAGTGGGATCGATAAAACTGGCTGACTAAACGGTATCCATCCTGTTg                                                                                                                                                  >  2:302415/1‑150 (MQ=255)
     aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTg                                                                                                                                                  <  2:925618/150‑1 (MQ=255)
      gcgcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTgg                                                                                                                                                 <  1:673345/150‑1 (MQ=255)
             gagCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGt                                                                                                                                          >  1:345469/1‑150 (MQ=255)
                cATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGtttt                                                                                                                                       <  2:656136/150‑2 (MQ=255)
                    gAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGc                                                                                                                                   >  1:936669/1‑150 (MQ=255)
                    gAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttttgtca                                                                                                                                   <  1:1563950/150‑6 (MQ=255)
                     aGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGcc                                                                                                                                  <  2:928208/150‑1 (MQ=255)
                         ccAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGtt                                                                                                                              >  2:1175223/1‑150 (MQ=255)
                         ccAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGttttgtcag                                                                                                                                  <  1:784791/146‑7 (MQ=38)
                          cAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTc                                                                                                                             <  1:1022689/150‑1 (MQ=255)
                          cAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTc                                                                                                                             <  2:1607726/150‑1 (MQ=255)
                           aaCGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCa                                                                                                                            <  2:470966/150‑1 (MQ=255)
                             cGGGGAAGGCTTAGTGAAAACGTGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGt                                                                                                                          <  2:400645/150‑1 (MQ=255)
                              ggggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGtt                                                                                                                         <  1:1058551/150‑1 (MQ=255)
                               gggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTg                                                                                                                        <  1:1059873/150‑1 (MQ=255)
                               gggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGATGGAAAGCGTTTCACGCCTGTTCAGTTg                                                                                                                        <  1:227996/150‑1 (MQ=255)
                                ggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGc                                                                                                                       <  2:514080/150‑1 (MQ=255)
                                ggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGc                                                                                                                       <  2:1231010/150‑1 (MQ=255)
                                  aaGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAg                                                                                                                     <  2:1057996/150‑1 (MQ=255)
                                  aaGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAg                                                                                                                     <  2:1510357/150‑1 (MQ=255)
                                    ggCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGa                                                                                                                   >  2:1447271/1‑150 (MQ=255)
                                     gCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGaa                                                                                                                  >  2:1625014/1‑150 (MQ=255)
                                     gCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGaa                                                                                                                  >  1:4219/1‑150 (MQ=255)
                                      cTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAAc                                                                                                                 >  1:1702568/1‑150 (MQ=255)
                                         aGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACaga                                                                                                              >  1:1393472/1‑150 (MQ=255)
                                         aGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACaga                                                                                                              >  1:606345/1‑150 (MQ=255)
                                          gTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACagag                                                                                                             <  2:345469/150‑1 (MQ=255)
                                           tGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGc                                                                                                            >  2:1404770/1‑150 (MQ=255)
                                           tGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGc                                                                                                            >  1:1409466/1‑150 (MQ=255)
                                                aCGGGTAGTAATGTTGTCATCATCAGCAATCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTcc                                                                                                       >  2:1430549/1‑150 (MQ=255)
                                                  gggTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtt                                                                                                     >  1:1181889/1‑150 (MQ=255)
                                                   ggTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCttt                                                                                                    >  1:428458/1‑150 (MQ=255)
                                                      agtaATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTa                                                                                                 >  2:588955/1‑150 (MQ=255)
                                                        taATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACg                                                                                               <  1:1625014/150‑1 (MQ=255)
                                                          aTGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtt                                                                                             <  2:349772/150‑1 (MQ=255)
                                                           tgttgtCATCCTCATCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGttt                                                                                            <  1:302415/150‑1 (MQ=255)
                                                            gttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTTACTAAACGGTATCCATCCTGTTGTATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                           <  1:701380/150‑1 (MQ=255)
                                                              tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                         >  1:747808/1‑150 (MQ=255)
                                                              tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                         >  1:314403/1‑150 (MQ=255)
                                                                tcatcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCa                                                                                       >  2:1402071/1‑150 (MQ=255)
                                                                tcatcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCa                                                                                       <  1:94623/150‑1 (MQ=255)
                                                                  atcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTAATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTAACGTTTTaaaaac                                                                                     >  1:1092559/1‑145 (MQ=255)
                                                                   tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTTGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACa                                                                                    <  2:1438796/150‑1 (MQ=255)
                                                                    catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAtgcacgc                                                                                   >  2:1627407/1‑144 (MQ=255)
                                                                    catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                   <  2:1175175/150‑1 (MQ=255)
                                                                    catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                   <  1:791936/150‑1 (MQ=255)
                                                                     atcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGt                                                                                  <  2:1274036/150‑1 (MQ=255)
                                                                       caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGATGCAGCAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                                >  1:228415/1‑150 (MQ=255)
                                                                         gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCa                                                                              <  1:1353140/150‑1 (MQ=255)
                                                                         gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCa                                                                              <  1:1638667/150‑1 (MQ=255)
                                                                                tGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGAt                                                                       >  2:260650/1‑150 (MQ=255)
                                                                                 ggCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATa                                                                      >  2:436512/1‑150 (MQ=255)
                                                                                  gCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAg                                                                     >  2:1536011/1‑150 (MQ=255)
                                                                                    cGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTa                                                                   >  1:1708836/1‑150 (MQ=255)
                                                                                    cGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCATTTTTACGTTTTATCAACAGTCTCAACCAGATAGTa                                                                   >  1:827777/1‑150 (MQ=255)
                                                                                       aaaTGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACa                                                                <  1:956153/150‑1 (MQ=255)
                                                                                           gAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTaa                                                            >  1:1482515/1‑150 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGTATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          <  2:1468786/150‑1 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          >  1:727934/1‑150 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          <  2:1393472/150‑1 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          >  2:463502/1‑150 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                          <  1:1536011/150‑1 (MQ=255)
                                                                                             aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATAAata                                                          <  1:1559403/150‑1 (MQ=255)
                                                                                                     tACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCAtgtt                                                  >  1:193151/1‑150 (MQ=255)
                                                                                                        tCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTc                                               <  1:36004/150‑1 (MQ=255)
                                                                                                           tAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTAt                                            >  1:1090736/1‑150 (MQ=255)
                                                                                                              gtggtaTCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGt                                         <  2:1592126/145‑1 (MQ=255)
                                                                                                                 tgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCa                                      <  2:606345/150‑1 (MQ=255)
                                                                                                                 tgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCa                                      >  1:1221988/1‑150 (MQ=255)
                                                                                                                     cTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTTGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGc                                  <  2:73528/150‑1 (MQ=255)
                                                                                                                     cTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTCTTCAGTTGCAGGAACAGAGCATACCTAATTACGTTTTATCAACAGTCTCAACAAGATAGTAACATAAATAATCATGTTGTCTATTTTTAAGAGc                                  >  1:1444936/1‑150 (MQ=255)
                                                                                                                      tataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGCTCAGAGCATTCCTTTTTACGTTTTATCATCAGGCTCAACCAGATAGTAACATTAAGATTCATGTTGTCTATGGTTCagagat                                 >  2:842625/1‑148 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                <  2:209821/150‑1 (MQ=255)
                                                                                                                       ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                >  2:407311/1‑150 (MQ=255)
                                                                                                                         aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTa                              <  1:410371/150‑1 (MQ=255)
                                                                                                                          aaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTaa                             <  2:618648/150‑1 (MQ=255)
                                                                                                                               ggCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtg                        >  2:439713/1‑150 (MQ=255)
                                                                                                                               ggCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtg                        >  1:960172/1‑150 (MQ=255)
                                                                                                                                 cTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtgtg                      <  1:137912/150‑1 (MQ=255)
                                                                                                                                  tGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgc                     >  2:1684024/1‑150 (MQ=255)
                                                                                                                                  tGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgc                     >  1:1008975/1‑150 (MQ=255)
                                                                                                                                  tGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgc                     <  2:747808/150‑1 (MQ=255)
                                                                                                                                     cAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATCAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGcc                  >  2:1384741/1‑150 (MQ=255)
                                                                                                                                       aaaCGCTATCCAACAGGATGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCtt                >  1:1225566/1‑150 (MQ=255)
                                                                                                                                        aaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCACCCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTc               >  2:80935/1‑150 (MQ=255)
                                                                                                                                        aaCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTc               >  2:1538021/1‑150 (MQ=255)
                                                                                                                                         aCGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCt              >  2:1124340/1‑150 (MQ=255)
                                                                                                                                            gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAAt           >  2:1542110/1‑150 (MQ=255)
                                                                                                                                            gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAAt           >  2:862509/1‑150 (MQ=255)
                                                                                                                                              atccatccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCa         <  2:428458/150‑1 (MQ=255)
                                                                                                                                                  atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg     <  2:1651671/150‑1 (MQ=255)
                                                                                                                                                  atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg     <  1:540100/150‑1 (MQ=255)
                                                                                                                                                   tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCTGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGc    >  2:1113144/1‑150 (MQ=255)
                                                                                                                                                   tccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGc    <  1:145901/150‑1 (MQ=255)
                                                                                                                                                     cTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAt  <  2:953741/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGACCAGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCAT  >  NC_000913/3630708‑3631006

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