Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,630,860 T→A 5.1% intergenic (‑258/+108) 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,8600TA 5.1% 464.4 / 8.2 154intergenic (‑258/+108)yhiI/yhiJputative membrane fusion protein (MFP) of efflux pump/DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base T (82/64);  new base A (2/6);  total (84/70)
Fisher's exact test for biased strand distribution p-value = 1.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.99e-01

Only 100 of 154 total aligned reads displayed.

AGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATA  >  NC_000913/3630713‑3631007
                                                                                                                                                   |                                                                                                                                                   
aGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTg                                                                                                                                                   <  2:925618/150‑1 (MQ=255)
 gcgcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTgg                                                                                                                                                  <  1:673345/150‑1 (MQ=255)
  cgcgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGa                                                                                                                                                 <  1:1417760/150‑1 (MQ=255)
   gcgcAGATCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGAt                                                                                                                                                <  1:1203346/150‑1 (MQ=255)
    cgcAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATa                                                                                                                                               >  1:209821/1‑150 (MQ=255)
      cAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGc                                                                                                                                             >  1:1231010/1‑150 (MQ=255)
         agCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGCTATCCAACAGGATGGATACCGtt                                                                                                                                          <  2:418853/150‑1 (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)
                     cAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTc                                                                                                                              <  1:1022689/150‑1 (MQ=255)
                      aaCGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATATCGTTTCACGCCTGTTCa                                                                                                                             >  2:1368365/1‑150 (MQ=255)
                      aaCGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCa                                                                                                                             <  2:470966/150‑1 (MQ=255)
                        cGGGGAAGGCTTAGTGAAAACGTGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGt                                                                                                                           <  2:400645/150‑1 (MQ=255)
                        cGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGt                                                                                                                           >  1:901099/1‑150 (MQ=255)
                        cGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCAAAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGt                                                                                                                           <  2:1036877/150‑1 (MQ=255)
                         ggggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGtt                                                                                                                          <  1:1058551/150‑1 (MQ=255)
                          gggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTg                                                                                                                         <  1:1059873/150‑1 (MQ=255)
                           ggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGc                                                                                                                        <  2:1231010/150‑1 (MQ=255)
                           ggAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGc                                                                                                                        <  2:514080/150‑1 (MQ=255)
                            gAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCa                                                                                                                       <  2:1622300/150‑1 (MQ=255)
                            gAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTACACGCCTGTTCAGTTGCa                                                                                                                       <  2:324714/150‑1 (MQ=255)
                             aaGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAg                                                                                                                      <  2:1057996/150‑1 (MQ=255)
                             aaGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAg                                                                                                                      <  2:1510357/150‑1 (MQ=255)
                                gCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGaa                                                                                                                   >  1:340456/1‑150 (MQ=255)
                                gCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGaa                                                                                                                   >  1:4219/1‑150 (MQ=255)
                                gCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGaa                                                                                                                   >  2:1625014/1‑150 (MQ=255)
                                 cTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAAc                                                                                                                  >  2:1559403/1‑150 (MQ=255)
                                     gTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACagag                                                                                                              <  2:345469/150‑1 (MQ=255)
                                      tGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGc                                                                                                             >  2:1404770/1‑150 (MQ=255)
                                             gggTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCtt                                                                                                      >  1:1181889/1‑150 (MQ=255)
                                                      tgttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGAATCCATACTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGttt                                                                                             >  2:810695/1‑150 (MQ=255)
                                                       gttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                            >  1:7947/1‑150 (MQ=255)
                                                        ttgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGtttt                                                                                            <  1:512955/149‑1 (MQ=255)
                                                         tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                          >  1:747808/1‑150 (MQ=255)
                                                         tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                          >  1:314403/1‑150 (MQ=255)
                                                         tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                          >  2:1453345/1‑150 (MQ=255)
                                                         tgtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTAt                                                                                          <  1:1719964/150‑1 (MQ=255)
                                                          gtCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATc                                                                                         >  2:36004/1‑150 (MQ=255)
                                                           tcatcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCa                                                                                        >  2:540100/1‑150 (MQ=255)
                                                             atcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTAGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAAc                                                                                      >  2:1636416/1‑150 (MQ=255)
                                                             atcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTAATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTAACGTTTTaaaaac                                                                                      >  1:1092559/1‑145 (MQ=255)
                                                              tcatcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTTGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACa                                                                                     <  2:1438796/150‑1 (MQ=255)
                                                               catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                    >  1:718930/1‑150 (MQ=255)
                                                               catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                    <  1:791936/150‑1 (MQ=255)
                                                               catcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAg                                                                                    <  2:1175175/150‑1 (MQ=255)
                                                                atcaGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGt                                                                                   <  2:1274036/150‑1 (MQ=255)
                                                                  caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                                 >  1:618648/1‑150 (MQ=255)
                                                                  caGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGATGCAGCAACAGAGCATTCCTTTTTACGTTTTATCAACAGtct                                                                                 >  1:228415/1‑150 (MQ=255)
                                                                    gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCa                                                                               <  1:1638667/150‑1 (MQ=255)
                                                                    gCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCa                                                                               <  1:418489/150‑1 (MQ=255)
                                                                       ttCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAAcc                                                                            <  1:154593/150‑1 (MQ=255)
                                                                        tCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCa                                                                           <  1:1402071/150‑1 (MQ=255)
                                                                           tGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGAt                                                                        >  2:260650/1‑150 (MQ=255)
                                                                            ggCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATa                                                                       >  2:436512/1‑150 (MQ=255)
                                                                               cGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTa                                                                    >  1:1708836/1‑150 (MQ=255)
                                                                                 tAAATGAAAGACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAAc                                                                  <  2:1441442/150‑1 (MQ=255)
                                                                                  aaaTGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACa                                                                 <  1:956153/150‑1 (MQ=255)
                                                                                      gAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTaa                                                             >  1:1482515/1‑150 (MQ=255)
                                                                                        aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                           <  1:1536011/150‑1 (MQ=255)
                                                                                        aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                           <  2:1393472/150‑1 (MQ=255)
                                                                                        aaCACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAata                                                           >  2:463502/1‑150 (MQ=255)
                                                                                           acGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACAGATAGTAACATTAATATTCa                                                       <  2:1409466/150‑1 (MQ=255)
                                                                                               aTACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCAtgt                                                    >  1:953741/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)
                                                                                                           gtgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTc                                        >  2:1141771/1‑150 (MQ=255)
                                                                                                            tgATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCa                                       <  2:606345/150‑1 (MQ=255)
                                                                                                             gATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCag                                      >  1:398135/1‑150 (MQ=255)
                                                                                                                cTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTTGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGc                                   <  2:73528/150‑1 (MQ=255)
                                                                                                                cTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGc                                   >  2:1395204/1‑150 (MQ=255)
                                                                                                                cTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTCTTCAGTTGCAGGAACAGAGCATACCTAATTACGTTTTATCAACAGTCTCAACAAGATAGTAACATAAATAATCATGTTGTCTATTTTTAAGAGc                                   >  1:1444936/1‑150 (MQ=255)
                                                                                                                 tataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGCTCAGAGCATTCCTTTTTACGTTTTATCATCAGGCTCAACCAGATAGTAACATTAAGATTCATGTTGTCTATGGTTCagagat                                  >  2:842625/1‑148 (MQ=255)
                                                                                                                  ataAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTg                                 <  2:209821/150‑1 (MQ=255)
                                                                                                                    aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTa                               >  2:1522960/1‑150 (MQ=255)
                                                                                                                    aaaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTa                               >  1:1651671/1‑150 (MQ=255)
                                                                                                                     aaaCTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTaa                              <  2:618648/150‑1 (MQ=255)
                                                                                                                          ggCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtg                         >  2:439713/1‑150 (MQ=255)
                                                                                                                            cTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATtgtg                       <  1:137912/150‑1 (MQ=255)
                                                                                                                             tGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgc                      >  1:1008975/1‑150 (MQ=255)
                                                                                                                             tGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgc                      <  2:747808/150‑1 (MQ=255)
                                                                                                                              gACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgcg                     <  2:1702568/150‑1 (MQ=255)
                                                                                                                              gACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTgcg                     <  1:260650/150‑1 (MQ=255)
                                                                                                                                cAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTAAGAGCTGTAACATTGTGCGcc                   >  2:1225272/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)
                                                                                                                                     cGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTa              >  2:1540789/1‑150 (MQ=255)
                                                                                                                                       gTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAAt            >  2:862509/1‑150 (MQ=255)
                                                                                                                                        tATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTATATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATc           <  1:1368365/150‑1 (MQ=255)
                                                                                                                                         atccatccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCa          <  2:428458/150‑1 (MQ=255)
                                                                                                                                          tccatccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCaa         >  1:749215/1‑150 (MQ=255)
                                                                                                                                             atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg      >  1:1612830/1‑150 (MQ=255)
                                                                                                                                             atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg      <  1:540100/150‑1 (MQ=255)
                                                                                                                                             atccTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAgg      <  2:1651671/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)
                                                                                                                                                   |                                                                                                                                                   
AGCGCGCAGAGCATTGAGCGCCAACGGGGAAGGCTTAGTGAAAACGGGTAGTAATGTTGTCATCATCAGCATTCCTGGCCGTAAATGAAACACGCATACTCATAAGTGTGATCTATAAAACTGGCTGACAAAACGGTATCCATCCTGTTGGATAGCGTTTCACGCCTGTTCAGTTGCAGGAACAGAGCATTCCTTTTTACGTTTTATCAACAGTCTCAACCAGATAGTAACATTAATATTCATGTTGTCTATGGTTCAGAGCTGTAACATTGTGCGCCTTCTAATCAAAGGCATA  >  NC_000913/3630713‑3631007

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: