Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,356,734 G→C 7.0% intergenic (+23/‑170) 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,7340GC 7.0% 343.4 / 12.9 100intergenic (+23/‑170)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 G (47/46);  new base C (2/5);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AGATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATC  >  NC_000913/2356587‑2356879
                                                                                                                                                   |                                                                                                                                                 
aGATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAGATGGGAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTACGc                                                                                                                                                  <  2:113225/149‑1 (MQ=255)
aGATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGc                                                                                                                                                  >  4:198878/1‑149 (MQ=255)
    agaagaAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCAtgt                                                                                                                                              >  4:1059135/1‑149 (MQ=255)
      aagaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTgg                                                                                                                                            <  1:373643/149‑1 (MQ=255)
       agaaAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGt                                                                                                                                           >  3:248065/1‑149 (MQ=255)
            gcggcgCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCg                                                                                                                                      <  3:364613/149‑1 (MQ=255)
            gcggcgCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCg                                                                                                                                      <  2:653054/149‑1 (MQ=255)
             cggcgCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                                                                     <  1:783839/149‑1 (MQ=255)
             cggcgCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGCCCCAGGCGCTGGCTTAATCTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                                                                     <  4:616429/149‑1 (MQ=255)
                 gcAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCg                                                                                                                                      >  3:117956/1‑144 (MQ=255)
                 gcAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCg                                                                                                                                      <  4:117956/144‑1 (MQ=255)
                  cAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                                                                     <  1:905219/144‑1 (MQ=255)
                  cAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAACCACATGCTTACCgg                                                                                                                                     >  2:905219/1‑144 (MQ=255)
                   aGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttt                                                                                                                               <  2:86338/149‑1 (MQ=255)
                     aTCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                             >  2:649521/1‑149 (MQ=255)
                     aTCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCAttttt                                                                                                                             >  4:782424/1‑149 (MQ=255)
                               cACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGa                                                                                                                   >  2:257015/1‑149 (MQ=255)
                                  cGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTc                                                                                                                >  1:448986/1‑149 (MQ=255)
                                   gACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTcc                                                                                                               <  4:667841/149‑1 (MQ=255)
                                     cTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGt                                                                                                             >  2:119229/1‑149 (MQ=255)
                                      tGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTg                                                                                                            >  4:417554/1‑149 (MQ=255)
                                         gAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGccc                                                                                                         >  3:895300/1‑149 (MQ=255)
                                         gAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGccc                                                                                                         >  3:496399/1‑149 (MQ=255)
                                            aCCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTTATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAAt                                                                                                      <  2:22230/149‑1 (MQ=255)
                                            aCCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAAt                                                                                                      <  2:909201/149‑1 (MQ=255)
                                            aCCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAAt                                                                                                      <  2:692278/149‑1 (MQ=255)
                                            aCCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAAt                                                                                                      <  4:770661/149‑1 (MQ=255)
                                              cTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGt                                                                                                    <  3:782424/149‑1 (MQ=255)
                                               tGTAATTGGCAGATTGAGATGTCCACAAGTCTTAGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCTCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTacac                                                                                                                     >  1:315152/1‑131 (MQ=255)
                                               tGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTacac                                                                                                                     <  2:315152/131‑1 (MQ=255)
                                                gTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTAt                                                                                                  <  4:894942/149‑1 (MQ=255)
                                                  aaaTGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                                                                <  3:586806/149‑1 (MQ=255)
                                                  aaaTGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGc                                                                                                <  3:398929/149‑1 (MQ=255)
                                                     tGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGtt                                                                                             <  1:119229/149‑1 (MQ=255)
                                                      ggCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTg                                                                                            <  3:943448/149‑1 (MQ=255)
                                                       gctgATTGAGATGTCCACAAGTCTGCGCAGTGAACATCCTATTACGCTACTGGCCCACGCGCTGTCTTGCACAGCTTATTGACTCCCGCTACGCGTGCGGTTACCGGGCGTTATTGCGTACACGTTTCCGGGCCGCATGTATGTGTTGc                                                                                           >  4:122297/4‑149 (MQ=255)
                                                       gCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGc                                                                                           <  1:620560/149‑1 (MQ=255)
                                                       gCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGc                                                                                           <  4:686241/149‑1 (MQ=255)
                                                           aTTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACg                                                                                       >  2:190055/1‑149 (MQ=255)
                                                              gagaTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAg                                                                                    >  3:1003443/1‑149 (MQ=255)
                                                                 aTGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGa                                                                                 >  1:1091162/1‑149 (MQ=255)
                                                                      cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaaa                                                                               <  2:147279/146‑1 (MQ=255)
                                                                      cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGaaa                                                                               >  1:147279/1‑146 (MQ=255)
                                                                      cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAtt                                                                            >  2:680299/1‑149 (MQ=255)
                                                                      cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAtt                                                                            <  1:834613/149‑1 (MQ=255)
                                                                      cacaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAAtt                                                                            >  4:937323/1‑149 (MQ=255)
                                                                       acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                           <  3:809314/149‑1 (MQ=255)
                                                                       acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                           <  1:613699/149‑1 (MQ=255)
                                                                       acaAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTc                                                                           >  3:337782/1‑149 (MQ=255)
                                                                            tCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTg                                                                      >  2:731549/1‑149 (MQ=255)
                                                                                     cGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtct                                                             <  3:568113/149‑1 (MQ=255)
                                                                                      gAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtctc                                                            >  1:794526/1‑149 (MQ=255)
                                                                                      gAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGtctc                                                            <  1:317104/149‑1 (MQ=255)
                                                                                          aTCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaa                                                        >  3:244375/1‑149 (MQ=255)
                                                                                           tCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa                                                       <  2:1104921/149‑1 (MQ=255)
                                                                                           tCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa                                                       >  4:704910/1‑149 (MQ=255)
                                                                                           tCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTACGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCaaa                                                       <  2:908192/149‑1 (MQ=255)
                                                                                             cGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAgg                                                     >  3:195454/1‑149 (MQ=255)
                                                                                              gATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGc                                                    >  2:701518/1‑149 (MQ=255)
                                                                                               aTTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCt                                                   <  4:496399/149‑1 (MQ=255)
                                                                                                ttACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTg                                                  <  3:895598/149‑1 (MQ=255)
                                                                                                    gCTACTGGCCCAGGCGCTGGCTTATACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCGTCTGGAAACGTCTCGCAAGGGCTGAAAc                                              >  3:912252/1‑149 (MQ=255)
                                                                                                    gCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAAc                                              >  4:236853/1‑149 (MQ=255)
                                                                                                       aCTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTgg                                           <  4:98155/149‑1 (MQ=255)
                                                                                                             ccAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATg                                     >  4:996198/1‑149 (MQ=255)
                                                                                                              cAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGt                                    >  1:957327/1‑149 (MQ=255)
                                                                                                               aGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTc                                   >  2:907884/1‑149 (MQ=255)
                                                                                                                  cgcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCaaa                                <  3:198878/149‑1 (MQ=255)
                                                                                                                  cgcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCaaa                                <  2:18158/149‑1 (MQ=255)
                                                                                                                  cgcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCaaa                                <  2:1059663/149‑1 (MQ=255)
                                                                                                                   gcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAg                               >  1:606712/1‑149 (MQ=255)
                                                                                                                   gcTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAg                               >  4:225808/1‑149 (MQ=255)
                                                                                                                     tGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGc                             <  1:257015/149‑1 (MQ=255)
                                                                                                                     tGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGc                             >  2:381865/1‑149 (MQ=255)
                                                                                                                        cTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTg                          <  1:649521/149‑1 (MQ=255)
                                                                                                                        cTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTg                          <  3:417554/149‑1 (MQ=255)
                                                                                                                         ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                         >  4:466466/1‑149 (MQ=255)
                                                                                                                         ttAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTgg                         <  4:248065/149‑1 (MQ=255)
                                                                                                                            aaCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGAt                      >  2:636093/1‑149 (MQ=255)
                                                                                                                             aCTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAATCTGGCAGATGTCAAAGGCCTGGGATa                     >  3:785069/1‑149 (MQ=255)
                                                                                                                              cTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATaa                    >  1:975766/1‑149 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGa                                       >  1:1088389/1‑129 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGa                                       <  2:1088389/129‑1 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAg                               <  1:983661/137‑1 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAg                               >  2:983661/1‑137 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTg                          >  2:770931/1‑142 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTg                          <  1:770931/142‑1 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                   <  3:790393/149‑1 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                   >  4:869594/1‑149 (MQ=255)
                                                                                                                               tCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAc                   >  3:958415/1‑149 (MQ=255)
                                                                                                                                 cTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCg                 >  1:639440/1‑149 (MQ=255)
                                                                                                                                  tttCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAAGCGt                >  2:862804/1‑149 (MQ=255)
                                                                                                                                  tttCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGt                >  1:491760/1‑149 (MQ=255)
                                                                                                                                        aTGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTc          <  3:705251/149‑1 (MQ=255)
                                                                                                                                        aTGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATCCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAGTTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTc          <  1:429733/149‑1 (MQ=255)
                                                                                                                                          gCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTcgc        <  3:1078415/149‑1 (MQ=255)
                                                                                                                                               gTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGtcat   <  3:508897/149‑1 (MQ=255)
                                                                                                                                                tAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCt                           >  2:340559/1‑124 (MQ=255)
                                                                                                                                                tAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCt                           <  1:340559/124‑1 (MQ=255)
                                                                                                                                                tAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGtcatc  <  4:195454/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
AGATAGAAGAAAGCGGCGCAGATCTGGTGGTCACCGACTGCGAAACCTGTAAATGGCAGATTGAGATGTCCACAAGTCTTCGCTGCGAACATCCGATTACGCTACTGGCCCAGGCGCTGGCTTAAACTCCTTTCTGATGCCCGGTAAGCATGTGGTTACCGGGCATTTTTGCGTACACGATTCCGTGCCCAATGTATGCGTTGCAACGCAGTGAAAATTCCTCTGAAAACGTCTCGCAAAGGCTGAAACTGGCAGATGTCAAAGGCCTGGGATAACCGTAATGTCGCGTCATC  >  NC_000913/2356587‑2356879

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