Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 775,620 C→G 7.7% L156L (CTC→CTG 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,6200CG 7.7% 276.9 / 10.9 89L156L (CTC→CTGtolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base C (40/42);  new base G (2/5);  total (42/47)
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.53e-01

GGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATT  >  NC_000913/775471‑775768
                                                                                                                                                     |                                                                                                                                                    
ggggCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTc                                                                                                                                                      >  1:1019328/1‑150 (MQ=255)
  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)
                                                                                                                                                     |                                                                                                                                                    
GGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATT  >  NC_000913/775471‑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: