Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,604,744 G→A 5.2% G300S (GGC→AGC)  hyfC → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,604,7440GA 5.2% 246.1 / 2.9 78G300S (GGC→AGC) hyfChydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base A (2/2);  total (38/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

CAAGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGG  >  NC_000913/2604603‑2604888
                                                                                                                                             |                                                                                                                                                
cAAGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAg                                                                                                                                           >  3:360001/1‑149 (MQ=255)
 aaGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGc                                                                                                                                          >  3:28366/1‑149 (MQ=255)
 aaGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGc                                                                                                                                          >  5:319108/1‑149 (MQ=255)
 aaGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAAAACGCTGGCACGCGGGCGTTTTTTACTCAATCACAAAGTGACCAGGCTTGGCTTCAGc                                                                                                                                          >  7:399131/1‑149 (MQ=255)
  aGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGcc                                                                                                                                         <  8:396339/149‑1 (MQ=255)
   gAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCt                                                                                                                                        <  5:61108/149‑1 (MQ=255)
    aaCTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCtt                                                                                                                                       >  5:176061/1‑149 (MQ=255)
       tttCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTtgct                                                                                                                                    >  1:106551/1‑149 (MQ=255)
         tctctcGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTtgct                                                                                                                                    <  5:202405/147‑1 (MQ=255)
           tctcGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTgcgc                                                                                                                                >  8:145755/1‑149 (MQ=255)
             tcGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGGTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTgcgc                                                                                                                                <  3:355220/147‑1 (MQ=255)
               gcctgcctGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTg                                                                                                                             >  3:352105/1‑148 (MQ=255)
               gcctgcctGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACAGGCCTCAATCGCAGAAACCACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTg                                                                                                                             >  8:197045/1‑148 (MQ=255)
                 ctgcctgcTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCAt                                                                                                                          >  5:263706/1‑149 (MQ=255)
                    cctgcTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATgg                                                                                                                        <  7:145755/148‑1 (MQ=255)
                        ctgACTTCACTTGTCGTTACGCTGCTCAAGG‑TTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCTCGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGtcttc                                                                                                                  <  5:343900/149‑1 (MQ=255)
                         tgACTTCACTTGTCGTTACGTGGCTCAATGTTTTTCTGATTTTTGTTCTTGCCTCAATCGCAGAAATCACGCTGGCACGCGGTCGTTTTTTACTCATTCTCCATGTGTCCTGGCTTGGCTTCAGCGTTGCTGCGCTTGCATGGGtcttc                                                                                                                  <  8:399131/149‑1 (MQ=255)
                             ttCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATgg                                                                                                                        >  7:80043/1‑139 (MQ=255)
                             ttCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATgg                                                                                                                        <  8:80043/139‑1 (MQ=255)
                             ttCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTAATCATTCACCATGTGACCTGGCTTGGCTTCAGCATTGCTGCGCTTGCATGGGTCTTCTGGt                                                                                                              >  2:19825/1‑149 (MQ=255)
                              tCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGt                                                                                                              >  4:6396/1‑148 (MQ=255)
                               cACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTAATCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTTCTGCGCTTGCATGGGTCTTCTGGtt                                                                                                             >  2:156181/1‑148 (MQ=255)
                               cACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATGTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGAATGGGTCTTCTGGt                                                                                                              <  5:5462/147‑1 (MQ=255)
                                aCTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTaa                                                                                                           <  5:270302/149‑1 (MQ=255)
                                aCTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTaa                                                                                                           >  1:34672/1‑149 (MQ=255)
                                 cTTGTCGTTCCGCTGCTCAAGGTTTTGCTGATGTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          <  6:73825/149‑1 (MQ=255)
                                 cTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGTCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          <  6:332240/149‑1 (MQ=255)
                                 cTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          <  3:6396/149‑1 (MQ=255)
                                 cTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          >  5:277353/1‑149 (MQ=255)
                                 cTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          <  4:360001/149‑1 (MQ=255)
                                     tCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGt                                                                                                      <  7:5862/149‑1 (MQ=255)
                                       gTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCt                                                                                                    <  3:261108/149‑1 (MQ=255)
                                       gTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCt                                                                                                    >  1:278101/1‑149 (MQ=255)
                                       gTTACGCTGCTAAAGGTTTTGATGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCt                                                                                                    >  7:2615/1‑149 (MQ=255)
                                          aCGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGt                                                                                                  <  4:352105/148‑1 (MQ=255)
                                          aCGCTCCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGt                                                                                                  <  2:382805/148‑1 (MQ=255)
                                           cGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTaa                                                                                                <  2:278101/149‑1 (MQ=255)
                                            gctgctCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTAGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAg                                                                                               >  6:213769/1‑149 (MQ=255)
                                            gctgctCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTAGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAg                                                                                               >  6:208532/1‑149 (MQ=255)
                                                   aaGGTTTTGATGATTATTGTACTGGCCTCAATCGCAGAAAAAACGCTGGCACGCGGGCGATTTTTACTCAATCACAATGTCAAATGGCTTGGCTTAATCCTTGCTGAGCTTGCATTGGTATTCTGGTTAAAAGGTATTTAAGGAGCACt                                                                                        >  6:179205/1‑149 (MQ=255)
                                                    aGGTTTTGCGGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACt                                                                                        >  5:110976/1‑148 (MQ=255)
                                                               aTTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAatat                                                                              <  3:347302/147‑1 (MQ=255)
                                                                      tACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCt                                                                     >  1:218602/1‑149 (MQ=255)
                                                                      tACTGGCCTAAATCGCAGAAAACACGCTGGCAAGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCt                                                                     <  4:207644/149‑1 (MQ=255)
                                                                         tGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTACCCGGTCTGTAAGGAGCACTGACGGATGATGGAAAATCTTGc                                                                  >  6:286586/1‑149 (MQ=255)
                                                                         tGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGc                                                                  >  6:281330/1‑149 (MQ=255)
                                                                         tGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGc                                                                  >  8:162034/1‑149 (MQ=255)
                                                                          ggCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          <  2:46929/108‑1 (MQ=255)
                                                                          ggCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAc                                                                                                          >  1:46929/1‑108 (MQ=255)
                                                                          ggCCTCAATCGAAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTTGCTTCCGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAAACGGACTGTAAGGAGCACTGACTGAATATGGAAAATCTTGc                                                                  >  8:393733/1‑148 (MQ=255)
                                                                              tCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTga                                                             >  4:68646/1‑149 (MQ=255)
                                                                              tCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTga                                                             >  7:327897/1‑149 (MQ=255)
                                                                                aaTCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTgac                                                            <  2:34672/148‑1 (MQ=255)
                                                                                 aTCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTgacga                                                          >  1:63433/1‑149 (MQ=255)
                                                                                  tCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTgacgac                                                         >  5:341561/1‑149 (MQ=255)
                                                                                    gCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGt                                                       <  3:125143/149‑1 (MQ=255)
                                                                                        aaaaCACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTAtt                                                   >  7:38978/1‑149 (MQ=255)
                                                                                             acGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATtgctgc                                              >  4:337900/1‑149 (MQ=255)
                                                                                                    cACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTAGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTAt                                       <  5:208532/149‑1 (MQ=255)
                                                                                                    cACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTAGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGGCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTAt                                       <  5:213769/149‑1 (MQ=255)
                                                                                                     aCGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATc                                      <  4:28366/149‑1 (MQ=255)
                                                                                                     aCGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATc                                      <  6:187677/149‑1 (MQ=255)
                                                                                                        cgGGCGTTTTTTACTCATTCACCATGTTACCTGGCTTTGATTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGgc                                   <  7:393733/149‑1 (MQ=255)
                                                                                                            cGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTGCAGCCTTGCTGGGCTTGCATGGGTCTTCTGGTTAACGGGTCTTTAAGGATCACTGACGGAATAATGAAAATCTTTCTCTGACGACGTTATTGCTGGCTTTTATCGGCGCAc                               <  5:179205/149‑1 (MQ=255)
                                                                                                             gTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACt                              >  7:342135/1‑149 (MQ=255)
                                                                                                                ttttACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTgg                            <  6:277353/148‑1 (MQ=255)
                                                                                                                  ttACTCATTCACAATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGtc                          <  6:110976/148‑1 (MQ=255)
                                                                                                                    aCTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTaa                                                                                                >  6:200213/1‑76 (MQ=255)
                                                                                                                    aCTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTaa                                                                                                <  5:200213/76‑1 (MQ=255)
                                                                                                                        aTTCACCATGTGCCCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGtcgtttc                     <  5:286586/147‑1 (MQ=255)
                                                                                                                        aTTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATtgct                                                >  4:7131/1‑120 (MQ=255)
                                                                                                                        aTTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATtgct                                                <  3:7131/120‑1 (MQ=255)
                                                                                                                        aTTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGtcgtttc                     <  2:106551/147‑1 (MQ=255)
                                                                                                                        aTTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGtcgtttc                     <  5:281330/147‑1 (MQ=255)
                                                                                                                            aCCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTc               <  7:162034/149‑1 (MQ=255)
                                                                                                                             ccATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCg              >  1:213733/1‑149 (MQ=255)
                                                                                                                              cATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGc             <  6:176061/149‑1 (MQ=255)
                                                                                                                                    aCCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCTCTGACGGAATATGGAAAATCTTGCGCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAAc       <  6:263706/149‑1 (MQ=255)
                                                                                                                                       tggcttggctTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTc    >  2:306293/1‑149 (MQ=255)
                                                                                                                                         gcttggcttCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTc    <  6:341561/147‑1 (MQ=255)
                                                                                                                                         gcttggcttCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCgg  <  2:218602/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                
CAAGAACTTTCTCTCGCCTGCCTGCTGACTTCACTTGTCGTTACGCTGCTCAAGGTTTTGCTGATTTTTGTACTGGCCTCAATCGCAGAAAACACGCTGGCACGCGGGCGTTTTTTACTCATTCACCATGTGACCTGGCTTGGCTTCAGCCTTGCTGCGCTTGCATGGGTCTTCTGGTTAACCGGTCTGTAAGGAGCACTGACGGAATATGGAAAATCTTGCTCTGACGACGTTATTGCTGCCTTTTATCGGCGCACTGGTCGTTTCGTTTTCGCCACAACGTCGG  >  NC_000913/2604603‑2604888

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