Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,356,737 T→C 5.4% intergenic (+26/‑167) glpC → / → yfaD anaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster/transposase_31 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,356,7370TC 5.4% 308.9 / 5.7 93intergenic (+26/‑167)glpC/yfaDanaerobic sn‑glycerol‑3‑phosphate dehydrogenase, C subunit, 4Fe‑4S iron‑sulfur cluster/transposase_31 family protein
Reads supporting (aligned to +/- strand):  ref base T (49/39);  new base C (2/3);  total (51/42)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.31e-01

ATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATCATAAA  >  NC_000913/2356589‑2356884
                                                                                                                                                    |                                                                                                                                                   
aTAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAt                                                                                                                                                     <  2:693298/149‑1 (MQ=255)
     agaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                                                                  >  2:49066/1‑147 (MQ=255)
     agaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgtg                                                                                                                                                  <  1:49066/147‑1 (MQ=255)
     agaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGt                                                                                                                                                >  4:684177/1‑149 (MQ=255)
      gaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                                                                   >  4:342340/1‑145 (MQ=255)
      gaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                                                                   <  3:342340/145‑1 (MQ=255)
      gaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGtt                                                                                                                                               >  1:509971/1‑149 (MQ=255)
           cggcgCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                                                                          <  1:45904/149‑1 (MQ=255)
               gcAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAt                                                                                                                                      >  2:450262/1‑149 (MQ=255)
                 aGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCGGGCAtca                                                                                                                                    <  2:181116/149‑3 (MQ=255)
                   aTCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                                  <  2:509971/149‑1 (MQ=255)
                   aTCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                                  <  4:638503/149‑1 (MQ=255)
                     cTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGc                                                                                                                                <  4:393058/149‑1 (MQ=255)
                     cTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGc                                                                                                                                <  3:387588/149‑1 (MQ=255)
                          ggtCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTaca                                                                                                                           >  4:704689/1‑149 (MQ=255)
                           gtCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTacac                                                                                                                          >  2:13944/1‑149 (MQ=255)
                             cACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGTATTTTTGCGTACACGa                                                                                                                        <  4:85018/149‑1 (MQ=255)
                              aCCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGAt                                                                                                                       >  2:777550/1‑149 (MQ=255)
                              aCCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGAt                                                                                                                       >  1:527892/1‑149 (MQ=255)
                              aCCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGAt                                                                                                                       >  3:877844/1‑149 (MQ=255)
                                   cTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGt                                                                                                                  <  3:97527/149‑1 (MQ=255)
                                   cTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGt                                                                                                                  <  2:527892/149‑1 (MQ=255)
                                             tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTa                                                                                                        >  2:35211/1‑149 (MQ=255)
                                             tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTa                                                                                                        >  4:770664/1‑149 (MQ=255)
                                             tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTa                                                                                                        >  2:43227/1‑149 (MQ=255)
                                                aaaTGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                                                                     <  3:353759/149‑1 (MQ=255)
                                                  aTGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGt                                                                                                   <  1:528713/149‑1 (MQ=255)
                                                     gCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGc                                                                                                >  4:497394/1‑149 (MQ=255)
                                                         aTTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACg                                                                                            >  3:339016/1‑149 (MQ=255)
                                                           tGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCa                                                                                          <  4:293593/149‑1 (MQ=255)
                                                             agaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGt                                                                                        <  1:13944/149‑1 (MQ=255)
                                                                 gTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCatacatt                                                                                    <  4:307432/149‑7 (MQ=255)
                                                                  tCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaaaa                                                                                   <  1:454018/149‑1 (MQ=255)
                                                                  tCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaaaa                                                                                   >  3:704762/1‑149 (MQ=255)
                                                                    cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAtt                                                                                 <  3:770664/149‑1 (MQ=255)
                                                                     acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTg                                                                                                 >  4:772956/1‑132 (MQ=255)
                                                                     acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTg                                                                                                 <  3:772956/132‑1 (MQ=255)
                                                                     acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                                >  2:658669/1‑149 (MQ=255)
                                                                     acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                                >  3:24443/1‑149 (MQ=255)
                                                                      caAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTcc                                                                               <  2:456735/149‑1 (MQ=255)
                                                                         gTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCctct                                                                            <  1:571779/149‑1 (MQ=255)
                                                                          tCTTCGCTGCGAACATCCGATTACGCTACTGGCCCATGAGCTGGCTTAAACTCCTTTCTGATGTCCGGTAAGCATGTGGTTGCTGGGCATTTGTTCGTACACGATCCCGTGCCCAATGTATGCGTTGCAACGCAGTGCAAATTCCTCTg                                                                           >  2:137657/1‑149 (MQ=255)
                                                                          tCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTg                                                                           >  3:431399/1‑149 (MQ=255)
                                                                               gCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAAc                                                                      <  3:684177/149‑1 (MQ=255)
                                                                                 tGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTCGCAACGCAGTGAAAATTCCTCTGAAAACGt                                                                    >  1:493273/1‑149 (MQ=255)
                                                                                        aTCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaa                                                             >  2:676836/1‑149 (MQ=255)
                                                                                            gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGc                                                         >  4:862613/1‑149 (MQ=255)
                                                                                            gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGc                                                         >  4:249391/1‑149 (MQ=255)
                                                                                            gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGc                                                         >  4:441336/1‑149 (MQ=255)
                                                                                                aCGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaa                                                     <  4:859329/149‑1 (MQ=255)
                                                                                                 cGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaaa                                                    >  2:264358/1‑149 (MQ=255)
                                                                                                 cGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGaaa                                                    <  2:222027/149‑1 (MQ=255)
                                                                                                    tACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAt                                                                                                                                                     <  1:98879/49‑1 (MQ=255)
                                                                                                    tACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAt                                                                                                                                                     >  2:98879/1‑49 (MQ=255)
                                                                                                     aCTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTgg                                                <  3:887983/149‑1 (MQ=255)
                                                                                                      cTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGc                                               <  4:704762/149‑1 (MQ=255)
                                                                                                       tGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCa                                              <  4:877844/149‑1 (MQ=255)
                                                                                                       tGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCa                                              >  3:243513/1‑149 (MQ=255)
                                                                                                       tGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCa                                              >  2:310034/1‑149 (MQ=255)
                                                                                                         gCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCa                                              <  1:414715/147‑1 (MQ=255)
                                                                                                         gCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCa                                              >  2:414715/1‑147 (MQ=255)
                                                                                                         gCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGa                                            <  3:497394/149‑1 (MQ=255)
                                                                                                          cccAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGAt                                           <  1:43227/149‑1 (MQ=255)
                                                                                                          cccAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGAt                                           <  1:264358/149‑1 (MQ=255)
                                                                                                            cAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGt                                         <  4:431399/149‑1 (MQ=255)
                                                                                                             aGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTc                                        >  4:81594/1‑149 (MQ=255)
                                                                                                             aGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTc                                        >  2:402079/1‑149 (MQ=255)
                                                                                                                     gCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCt                                >  3:777392/1‑149 (MQ=255)
                                                                                                                     gCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCt                                >  1:227501/1‑149 (MQ=255)
                                                                                                                       ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                              <  1:402079/149‑1 (MQ=255)
                                                                                                                       ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                              >  4:413940/1‑149 (MQ=255)
                                                                                                                       ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                              >  4:134480/1‑149 (MQ=255)
                                                                                                                       ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                              >  1:421641/1‑149 (MQ=255)
                                                                                                                           aCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATa                          >  1:935882/1‑149 (MQ=255)
                                                                                                                            cTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              >  1:768711/1‑96 (MQ=255)
                                                                                                                            cTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              <  2:768711/96‑1 (MQ=255)
                                                                                                                            cTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATaa                         >  1:504295/1‑149 (MQ=255)
                                                                                                                             tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                        >  1:556535/1‑149 (MQ=255)
                                                                                                                             tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                        >  2:232495/1‑149 (MQ=255)
                                                                                                                             tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                        <  1:232495/149‑1 (MQ=255)
                                                                                                                             tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                        >  3:344225/1‑149 (MQ=255)
                                                                                                                              ccTTTCTGATGCCCGGTAAGCATGTGGTTACCGGTCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGc                                  <  2:452068/138‑1 (MQ=255)
                                                                                                                              ccTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGc                                  >  1:452068/1‑138 (MQ=255)
                                                                                                                              ccTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAcc                       >  3:128387/1‑149 (MQ=255)
                                                                                                                                tttCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              >  1:881416/1‑92 (MQ=255)
                                                                                                                                tttCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              <  2:881416/92‑1 (MQ=255)
                                                                                                                                  tCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTaa                   <  4:158689/149‑1 (MQ=255)
                                                                                                                                       tGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGc                                               <  4:506466/116‑1 (MQ=39)
                                                                                                                                       tGCCCGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGc                                               >  3:506466/1‑116 (MQ=39)
                                                                                                                                        gCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCATCGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTcgc             >  4:774732/1‑149 (MQ=255)
                                                                                                                                        gCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTc               >  1:893322/1‑147 (MQ=255)
                                                                                                                                        gCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTc               <  2:893322/147‑1 (MQ=255)
                                                                                                                                         cccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              <  2:102085/83‑1 (MQ=255)
                                                                                                                                         cccGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCct                                                                              >  1:102085/1‑83 (MQ=255)
                                                                                                                                          ccGGTAACCACATGCTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGt           >  3:382158/1‑149 (MQ=255)
                                                                                                                                           cGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGtc          <  1:658669/149‑1 (MQ=255)
                                                                                                                                              tAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGtcatc       <  1:35211/149‑1 (MQ=255)
                                                                                                                                                 gCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATCATa    <  4:783772/149‑1 (MQ=255)
                                                                                                                                                   aTGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATCATaaa  >  1:879701/1‑149 (MQ=255)
                                                                                                                                                   aTGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATCATaaa  >  2:414920/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
ATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATCATAAA  >  NC_000913/2356589‑2356884

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