Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,068,924 G→A 8.6% intergenic (‑116/+23) argO ← / ← mscS L‑arginine exporter/small conductance mechanosensitive channel MscS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,068,9240GA 8.6% 308.5 / 11.6 103intergenic (‑116/+23)argO/mscSL‑arginine exporter/small conductance mechanosensitive channel MscS
Reads supporting (aligned to +/- strand):  ref base G (51/43);  new base A (5/4);  total (56/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 107 total aligned reads displayed.

CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068775‑3069071
                                                                                                                                                     |                                                                                                                                                   
cAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagagag                                                                                                                                                     <  2:799112/150‑1 (MQ=255)
cAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCagagag                                                                                                                                                     <  1:181841/150‑1 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                <  2:1180501/150‑1 (MQ=255)
     gCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtatta                                                                                                                                                >  1:1167452/1‑150 (MQ=255)
      cAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGtattat                                                                                                                                               >  1:1320345/1‑150 (MQ=255)
       aaGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATc                                                                                                                                              >  2:467836/1‑150 (MQ=255)
        aGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCt                                                                                                                                             >  2:298334/1‑150 (MQ=255)
          aCCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGATAGAGGCCAGAGAGTATTATCTgg                                                                                                                                           <  1:374450/150‑1 (MQ=255)
           ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGc                                                                                                                                          <  2:441613/150‑1 (MQ=255)
           ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGc                                                                                                                                          <  1:354515/150‑1 (MQ=255)
           ccTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGc                                                                                                                                          >  1:326455/1‑150 (MQ=255)
            cTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGcc                                                                                                                                         <  2:665642/150‑1 (MQ=255)
             ttGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGATAATACTCTCTGGCCt                                                                                                                                        <  1:254056/150‑1 (MQ=255)
             ttGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCt                                                                                                                                        >  1:907936/1‑150 (MQ=255)
                     aataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTg                                                                                                                                <  1:1583923/150‑1 (MQ=255)
                     aataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCCGGATTCTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCATCAGCGTTg                                                                                                                                >  1:1311409/1‑150 (MQ=255)
                      ataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGa                                                                                                                               <  2:1561913/150‑1 (MQ=255)
                      ataaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGa                                                                                                                               >  1:1550744/1‑150 (MQ=255)
                       taaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAt                                                                                                                              <  1:1114884/150‑1 (MQ=255)
                        aaGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGAtt                                                                                                                             >  1:1395359/1‑150 (MQ=255)
                             aaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                                                                        <  2:373803/150‑1 (MQ=255)
                             aaCACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCa                                                                                                                        <  2:326455/150‑1 (MQ=255)
                                 cGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCtt                                                                                                                    <  1:1544951/150‑1 (MQ=255)
                                 cGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCtt                                                                                                                    <  2:1304664/150‑1 (MQ=255)
                                 cGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCtt                                                                                                                    <  2:1316279/150‑1 (MQ=255)
                                         atTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGAATGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGtcttct                                                                                                            >  1:1194442/1‑150 (MQ=255)
                                          ttCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTcttctt                                                                                                           <  1:765863/150‑1 (MQ=255)
                                           tCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCttt                                                                                                          <  1:469689/150‑1 (MQ=255)
                                            ccGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTc                                                                                                         >  2:1522340/1‑150 (MQ=255)
                                             cGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCa                                                                                                        <  2:901777/150‑1 (MQ=255)
                                                   aaGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGGATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCt                                                                                                  <  2:1548186/150‑1 (MQ=255)
                                                          tGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGtt                                                                                           >  2:523594/1‑150 (MQ=255)
                                                              gCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTcaca                                                                                       >  2:1474579/1‑150 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGATGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                      <  2:709253/150‑1 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                      >  1:1233465/1‑150 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                      >  1:94973/1‑150 (MQ=255)
                                                               cAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACAt                                                                                      >  1:878065/1‑150 (MQ=255)
                                                                 gagaCTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATcc                                                                                    <  2:387438/150‑1 (MQ=255)
                                                                      tataTCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTg                                                                               <  1:1382182/150‑1 (MQ=255)
                                                                       atatCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGc                                                                              >  2:1001537/1‑150 (MQ=255)
                                                                        tatCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCg                                                                             >  1:1124328/1‑150 (MQ=255)
                                                                         atCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCgg                                                                            <  1:953128/150‑1 (MQ=255)
                                                                            gcgGGATTGTTATTAGCGGAAATTGATAATATTAAGTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTa                                                                         >  2:811405/1‑150 (MQ=255)
                                                                              gggATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACg                                                                       >  2:179564/1‑150 (MQ=255)
                                                                                gATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACggg                                                                     <  2:907936/150‑1 (MQ=255)
                                                                                 aTTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGa                                                                    <  1:467836/150‑1 (MQ=255)
                                                                                  ttGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCg                                                                             >  1:696482/1‑140 (MQ=255)
                                                                                  ttGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCg                                                                             <  2:696482/140‑1 (MQ=255)
                                                                                   tGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAg                                                                  >  2:253787/1‑150 (MQ=255)
                                                                                      tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                                                               >  2:590411/1‑150 (MQ=255)
                                                                                      tattaGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTg                                                               >  2:782688/1‑150 (MQ=255)
                                                                                         taGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATa                                                            <  1:1416104/150‑1 (MQ=255)
                                                                                          aGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCTTTTGCGGGTACGGGAAGCTGATAc                                                           >  1:1011636/1‑150 (MQ=255)
                                                                                             ggAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcgg                                                        >  2:119391/1‑150 (MQ=255)
                                                                                              gAAATTGATAATATTAATTGGCATTCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGGCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggc                                                       >  2:869756/1‑150 (MQ=255)
                                                                                              gAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggc                                                       <  2:545939/150‑1 (MQ=255)
                                                                                              gAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggc                                                       >  1:1240526/1‑150 (MQ=255)
                                                                                               aaaTTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACcggcg                                                      <  1:590411/150‑1 (MQ=255)
                                                                                                  ttGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCa                                                   <  2:1124328/150‑1 (MQ=255)
                                                                                                   tGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCATCTTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAt                                                  <  1:559388/150‑1 (MQ=255)
                                                                                                   tGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAt                                                  >  2:1548526/1‑150 (MQ=255)
                                                                                                   tGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTCTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCAt                                                  <  2:1395359/150‑1 (MQ=255)
                                                                                                     ataataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCa                                                <  2:1447411/150‑1 (MQ=255)
                                                                                                      taataTTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGCGAGTATTATCTGGCCGCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCggcggcgaaag                                               <  2:686514/150‑6 (MQ=255)
                                                                                                        atatTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAAt                                             <  2:1503477/150‑1 (MQ=255)
                                                                                                             aaTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACg                                        >  1:479442/1‑150 (MQ=255)
                                                                                                             aaTTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACg                                        <  1:782688/150‑1 (MQ=255)
                                                                                                                 ggCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCTTTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTa                                    <  2:1011636/150‑1 (MQ=255)
                                                                                                                 ggCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTa                                    >  1:467899/1‑150 (MQ=255)
                                                                                                                  gCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTaa                                   <  1:87964/150‑1 (MQ=255)
                                                                                                                   cAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAt                                  >  1:824467/1‑150 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGc                              >  2:1054370/1‑150 (MQ=255)
                                                                                                                       cAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGc                              >  2:2170/1‑150 (MQ=255)
                                                                                                                        aGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCt                             >  1:377555/1‑150 (MQ=255)
                                                                                                                         gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTc                            >  2:388840/1‑150 (MQ=255)
                                                                                                                         gAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGACTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTc                            >  1:829953/1‑150 (MQ=255)
                                                                                                                            aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAg                         <  1:173045/150‑1 (MQ=255)
                                                                                                                            aaaCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAg                         <  1:179564/150‑1 (MQ=255)
                                                                                                                              aCTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGca                       <  1:298334/150‑1 (MQ=255)
                                                                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGcac                      >  1:230022/1‑150 (MQ=255)
                                                                                                                               cTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGcac                      <  2:1167452/150‑1 (MQ=255)
                                                                                                                                 aataaGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACAt                    >  2:1066398/1‑150 (MQ=255)
                                                                                                                                   caacgctgagGCCCGATAATACTCTCTGGCCTCAGCGTTGAGTCCGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCAGTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAAGACGCTCCAGCACATcc                  >  2:796993/8‑150 (MQ=255)
                                                                                                                                      gagagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg               >  1:623977/1‑150 (MQ=255)
                                                                                                                                      cgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg               >  1:1208376/5‑150 (MQ=255)
                                                                                                                                      cgctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAg               >  1:162830/5‑150 (MQ=255)
                                                                                                                                       gctgagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGt              >  1:1279881/4‑150 (MQ=255)
                                                                                                                                       agagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGt              <  2:364028/150‑1 (MQ=255)
                                                                                                                                        gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa             >  2:340200/1‑150 (MQ=255)
                                                                                                                                        gagagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTa             >  1:177267/1‑150 (MQ=255)
                                                                                                                                          gagGCCAGATAATACTCTCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTaca           >  2:161801/1‑150 (MQ=255)
                                                                                                                                          gagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                     >  2:1388592/1‑124 (MQ=255)
                                                                                                                                          gagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGttt                                     <  1:1388592/124‑1 (MQ=255)
                                                                                                                                          gagGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTaca           <  2:467899/150‑1 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg         >  1:199468/1‑150 (MQ=255)
                                                                                                                                            ggCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACg         >  2:347581/1‑150 (MQ=255)
                                                                                                                                              ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGtt       >  2:1117816/1‑150 (MQ=255)
                                                                                                                                              ccAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGtt       >  2:1624372/1‑150 (MQ=255)
                                                                                                                                                   gagTATTATCTGGCCTGAGGGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACAGGTTTTGCa  <  2:1239917/150‑1 (MQ=255)
                                                                                                                                                   gagTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCa  <  1:523594/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CAAGTGCAAGACCTTGAAAGTAATAAGAAAACACGTTGGATATTCCGAATTAAGTATCTGATGCAGAGACTATATCGCGGGATTGTTATTAGCGGAAATTGATAATATTAATTGGCAATCAGAAAAACTAATAAGAGAGAGGCCAGAGAGTATTATCTGGCCTCAGCGTTGATTACGCAGCTTTGTCTTCTTTCACCCGCTTAAAGTTCACATCCATTTGCGGGTACGGGAAGCTGATACCGGCGGCATCAAATTCACGTTTAATACGCTCCAGCACATCCCAGTACACGTTTTGCA  >  NC_000913/3068775‑3069071

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