Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 775,613 T→C 6.9% I154T (ATC→ACC)  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,6130TC 6.9% 245.9 / 7.5 88I154T (ATC→ACC) tolQmembrane spanning protein in TolA‑TolQ‑TolR complex
Reads supporting (aligned to +/- strand):  ref base T (39/42);  new base C (2/4);  total (41/47)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

TGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAA  >  NC_000913/775466‑775761
                                                                                                                                                   |                                                                                                                                                    
tGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCg                                                                                                                                                    <  1:601583/150‑1 (MQ=255)
 ggAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGc                                                                                                                                                   >  1:410061/1‑150 (MQ=255)
   aaGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGccc                                                                                                                                                 >  1:733547/1‑150 (MQ=255)
     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                                                                                                                             <  2:1347107/150‑1 (MQ=255)
                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                             >  2:1276454/1‑150 (MQ=255)
                       cGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCa                                                                                                                             >  1:1419606/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                                                                                           <  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                                                    >  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:1023789/1‑150 (MQ=255)
                                                                                                                gTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTa                                    >  2:1158850/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              >  1:575931/1‑150 (MQ=255)
                                                                                                                                      aTGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTa              >  2:889997/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)
                                                                                                                                                   |                                                                                                                                                    
TGGAAGGGGCGTCGCGTGCTATGCGTATCTCCATGAACCGTGAACTTGAAAATCTGGAAACGCACATTCCGTTCCTCGGTACGGTTGGCTCCATCAGCCCGTATATTGGTCTGTTTGGTACGGTCTGGGGGATCATGCACGCCTTTATCGCCCTCGGGGCGGTAAAACAAGCAACACTGCAAATGGTTGCGCCCGGTATCGCAGAAGCGTTGATTGCGACTGCAATTGGTCTGTTTGCCGCTATCCCGGCAGTTATGGCCTACAACCGCCTCAACCAGCGCGTAAACAAACTGGAA  >  NC_000913/775466‑775761

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: