Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 775,627 G→A 7.7% V159I (GTA→ATA)  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,6270GA 7.7% 279.1 / 7.7 91V159I (GTA→ATA) tolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base G (41/43);  new base A (2/5);  total (43/48)
Fisher's exact test for biased strand distribution p-value = 4.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.46e-01

TGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAAC  >  NC_000913/775482‑775776
                                                                                                                                                 |                                                                                                                                                     
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:1347107/150‑1 (MQ=255)
       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                                            >  2:1276454/1‑150 (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                                                                                                          <  2:1112033/150‑1 (MQ=255)
                                         aaaCGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAg                                                                                                          <  1:695082/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                                                                   >  2:1131817/1‑150 (MQ=255)
                                                                                gCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATcc                                                                   <  1:595560/150‑1 (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                          >  1:191309/1‑150 (MQ=255)
                                                                                                                         cACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTaaac                          >  2:1093663/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)
                                                                                                                                             ggCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGa      >  1:46258/1‑150 (MQ=255)
                                                                                                                                             ggCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGa      <  1:1131817/150‑1 (MQ=255)
                                                                                                                                              gCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGAc     >  2:1334507/1‑150 (MQ=255)
                                                                                                                                              gCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGAc     <  1:1082315/150‑1 (MQ=255)
                                                                                                                                               cGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACa    <  1:397611/150‑1 (MQ=255)
                                                                                                                                                 gTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAAc  >  2:1634753/1‑150 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
TGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAACTGAATTACGACAAC  >  NC_000913/775482‑775776

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: