Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 775,621 G→A 8.3% G157R (GGG→AGG)  tolQ → membrane spanning protein in TolA‑TolQ‑TolR complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913775,6210GA 8.3% 251.8 / 8.3 86G157R (GGG→AGG) tolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base A (2/5);  total (41/45)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.65e-01

GGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATT  >  NC_000913/775473‑775768
                                                                                                                                                    |                                                                                                                                                   
ggcggcGCGAGCGATGCGTATCGCCATGACCCGTGACCTAGAAAAGCTGGACACGCACATGCCGGTCCTCGGGACGGTTGGCTCCCTCAGCCCGTATATTGGTCTGTGTGGTACGGGCTGGGGGATCATGCACGCCTTTATCGCCCTCgg                                                                                                                                                    <  1:325634/150‑1 (MQ=255)
 gCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCggg                                                                                                                                                   >  2:1541890/1‑150 (MQ=255)
 gCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCCGGGGGATCATGCACGCCTTTATCGCCCTCggg                                                                                                                                                   <  2:1495839/150‑1 (MQ=255)
         tGCTATGCTTATCTCCATGAACCGTGAACTTGAAAATCTGCATACGCACATGCCGTTCCGCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTaaa                                                                                                                                           <  1:1141877/150‑1 (MQ=255)
           cTATGCGTATCTCCATTAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTTGTCTGTTTGGTACGGTGTGGGGGATCATGCACGCCTTTATCGCCCCGAGGGCGATaaagg                                                                                                                                         <  2:849339/150‑3 (MQ=255)
           cTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTTGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaagg                                                                                                                                         <  2:551308/150‑3 (MQ=255)
           cTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaagg                                                                                                                                         <  1:722800/150‑3 (MQ=255)
            tATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACa                                                                                                                                        >  1:1112033/1‑150 (MQ=255)
             aTGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACaa                                                                                                                                       >  1:758139/1‑150 (MQ=255)
                cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                                    >  1:1419606/1‑150 (MQ=255)
                cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                                    >  2:1276454/1‑150 (MQ=255)
                cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                                    <  2:1347107/150‑1 (MQ=255)
                cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaagg                                                                                                                                         >  2:617591/1‑143 (MQ=255)
                cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTACCGCCCCGAGGGCGATaaagg                                                                                                                                         <  1:617591/145‑3 (MQ=255)
                  tATCTCCATGACCCGTGAACTTGAAAATCTGGAAACGCTCATTCCGTGCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACCGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAac                                                                                                                                  <  1:86964/150‑1 (MQ=255)
                   aTCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAaca                                                                                                                                 >  2:13651/1‑150 (MQ=255)
                      tcCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTg                                                                                                                              >  2:1238928/1‑150 (MQ=255)
                       ccATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGc                                                                                                                             <  2:1419606/150‑1 (MQ=255)
                         aTGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCaa                                                                                                                           <  2:1271191/150‑1 (MQ=255)
                          tGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCaaa                                                                                                                          >  1:1615876/1‑150 (MQ=255)
                            aaCCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATg                                                                                                                        >  2:1276213/1‑150 (MQ=255)
                              ccGTGAACTTGAAATTCTGCAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCGCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGTGCGGTAAAACAAGCAACACTGCAAATCGt                                                                                                                      <  2:536849/150‑1 (MQ=255)
                                gTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTg                                                                                                                    <  1:811572/150‑1 (MQ=255)
                                 tGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTgc                                                                                                                   <  2:410061/150‑1 (MQ=255)
                                  taaCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACTCTGCAAATGGTTgcg                                                                                                                  <  1:814566/149‑1 (MQ=255)
                                    aCTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGcc                                                                                                                >  1:1266501/1‑150 (MQ=255)
                                       tGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCgg                                                                                                             >  1:542429/1‑150 (MQ=255)
                                           aaTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATc                                                                                                         >  1:497333/1‑150 (MQ=255)
                                            aTCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCg                                                                                                        <  1:1432799/150‑1 (MQ=255)
                                              cTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCa                                                                                                      <  1:1443128/150‑1 (MQ=255)
                                                  attcGCTCATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCGGTATATTGGTCTGTTTGGTTCTGTCTGGGGGATCAGGCACGCCTTTATCGCCCTCGGGGCGGTATAAGAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAg                                                                                                  <  1:701260/147‑1 (MQ=255)
                                                  aaaCGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAg                                                                                                  <  1:695082/150‑1 (MQ=255)
                                                  aaaCGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAg                                                                                                  <  2:1112033/150‑1 (MQ=255)
                                                        acaTTCCGTTCCTCGGTACGTGTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGa                                                                                            <  1:942642/150‑1 (MQ=255)
                                                         caTTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGAt                                                                                           >  1:1644151/1‑150 (MQ=255)
                                                           ttCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTg                                                                                         <  1:1707890/150‑1 (MQ=255)
                                                            tCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGGACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGc                                                                                        <  2:1136825/150‑1 (MQ=255)
                                                             ccGTTCCTCGGTACGGTTTGCTCCATCAGCCCGTATATTGGTCTGTTTGGTTCGGTCTGGGGGATCATGCACGCCTTTATCGCACTCGGGGCGGTAAAACAAGCAACACTGCAAAAGGTTGCGACCGGTATCCCAGAAGCGATGATTGcc                                                                                       >  1:471714/1‑149 (MQ=255)
                                                             ccGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTTGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCAGGTATCGCAGAAGCGTTGATTGCg                                                                                       <  1:1016363/150‑1 (MQ=255)
                                                                   cTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTATATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCa                                                                                 >  2:397611/1‑150 (MQ=255)
                                                                      ggTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAAtt                                                                              <  2:497333/150‑1 (MQ=255)
                                                                          cGGTTGTCTCCATCAGCCCGTATAGTGGTCTGTGTGGTACGGGCTGGGGGAGCATGCACGCCTTAAACGCCCACGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCTCAGAAGCGTTGATTGCGACTGCAATTGGTc                                                                          <  2:51660/150‑1 (MQ=255)
                                                                          cGGTTGGCTCCATCAGCCAGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACCAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTc                                                                          >  2:215927/1‑150 (MQ=255)
                                                                            gTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTg                                                                        <  1:1653689/150‑1 (MQ=255)
                                                                              tGGCTCCATCAGCCCGGATATTGGTCTGTTTGGTACGGTCAGGGGGATCCTGCACGCCTTTATCGACCTCCGGGCGGTAAAACAAGCAACCCTGCAGATGGTGTCGCCCGGTATCGAAGAAGCGTTGTGTGCGAATGCATTTGCTCGGtt                                                                      >  1:800534/1‑150 (MQ=255)
                                                                               ggCGCCATCAGCCCGTATAGTGGTCTGTGTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGttt                                                                     <  2:327422/150‑1 (MQ=255)
                                                                                gCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTg                                                                    <  1:1276454/150‑1 (MQ=255)
                                                                                 ctacatcatCCCGTATATTGGGCTGTTTCGTACGGTGGGGGGGATAATGCACTCCTTTTTCGCCCTCGGGGCGGTAAAACAACCAACACTGCAAAGGGTTGCGCCAGGTATAGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGc                                                                   <  2:800534/147‑1 (MQ=255)
                                                                                 cTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAACCGTTGATTGCGACTGCAATTGGTCTGTTTGc                                                                   >  2:1507909/1‑150 (MQ=255)
                                                                                  tCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGcc                                                                  <  1:1276213/150‑1 (MQ=255)
                                                                                   ccATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCg                                                                 <  1:1541890/150‑1 (MQ=255)
                                                                                    cATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGc                                                                >  1:670803/1‑150 (MQ=255)
                                                                                         gCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATcc                                                           <  1:595560/150‑1 (MQ=255)
                                                                                         gCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATcc                                                           >  2:1131817/1‑150 (MQ=255)
                                                                                            cGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCgg                                                        >  2:914047/1‑150 (MQ=255)
                                                                                             gTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGc                                                       <  2:758139/150‑1 (MQ=255)
                                                                                               atatTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAg                                                     <  2:1419700/150‑1 (MQ=255)
                                                                                                    ggTCTGTTTGGTACGGTCTGGGTGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTCTCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATg                                                <  1:888870/150‑1 (MQ=255)
                                                                                                    ggTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGGGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATg                                                <  1:810853/150‑1 (MQ=255)
                                                                                                    ggTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATg                                                >  2:319750/1‑150 (MQ=255)
                                                                                                     gTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATgg                                               >  2:1712039/1‑150 (MQ=255)
                                                                                                       cTGTTTGGTACGGTCTGGGGTATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGcc                                             >  2:196969/1‑150 (MQ=255)
                                                                                                         gTTTGGTTCGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                           <  1:1013516/150‑1 (MQ=255)
                                                                                                         gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                           >  2:1158850/1‑150 (MQ=255)
                                                                                                         gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                           >  2:1023789/1‑150 (MQ=255)
                                                                                                            tGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGACGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCAACac                                        >  1:78901/1‑149 (MQ=255)
                                                                                                             ggTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAAc                                       >  2:1082315/1‑150 (MQ=255)
                                                                                                               tACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCACATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCg                                     <  2:1019328/150‑1 (MQ=255)
                                                                                                                 cGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGcc                                   <  1:1238928/150‑1 (MQ=255)
                                                                                                                  ggTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCt                                  >  2:528370/1‑150 (MQ=255)
                                                                                                                    gcTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCa                                <  1:1658553/149‑1 (MQ=255)
                                                                                                                       gggggATCATGCACGCCTGTATCGCCCTCGGGGCGGTAAAAAAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGTGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAAcc                             <  2:1278623/150‑1 (MQ=255)
                                                                                                                          ggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGCCTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAgc                          >  1:859362/1‑150 (MQ=255)
                                                                                                                          ggATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAgc                          >  1:1448122/1‑150 (MQ=255)
                                                                                                                           gATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACATGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAgcg                         <  2:1463429/150‑1 (MQ=255)
                                                                                                                               aTGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTa                     >  2:889997/1‑150 (MQ=255)
                                                                                                                               aTGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTa                     >  1:575931/1‑150 (MQ=255)
                                                                                                                                  cACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaac                  >  2:1093663/1‑150 (MQ=255)
                                                                                                                                  cACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaac                  >  1:191309/1‑150 (MQ=255)
                                                                                                                                   cttgttttACCGCCCCGAGGGCGATAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCACTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaaca                 >  1:89704/6‑150 (MQ=255)
                                                                                                                                   aCGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaaca                 <  1:1712039/150‑1 (MQ=255)
                                                                                                                                        tttACCGCCCCGAGGGCGATAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCACTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTg            <  2:89704/150‑1 (MQ=255)
                                                                                                                                         ttATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTgg           <  2:1266501/150‑1 (MQ=255)
                                                                                                                                           aTCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACACCCGCCTCAACCAGCGCGTAAACAAACTGGaa         >  1:902123/1‑150 (MQ=255)
                                                                                                                                             cGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACt       >  1:1395501/1‑150 (MQ=255)
                                                                                                                                              gCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACACCCGCCTCACCCAGCTCGTAAACAAACTGGAACtt      >  2:1083523/1‑149 (MQ=255)
                                                                                                                                                ccTCGGGGCGGTAAAACAAGCAACCCTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGaa    <  1:73844/150‑1 (MQ=255)
                                                                                                                                                ccTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGaa    >  1:1556398/1‑150 (MQ=255)
                                                                                                                                                  tCGGGGCGGTAATACAAGCAACACTGCAAATGGTTGCGCCAGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAGCCGCCTCAACCAGCGCGTAAACAAACTGGACCTGAAtt  >  2:120247/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATT  >  NC_000913/775473‑775768

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