Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,849,674 C→G 6.0% C59S (TGC→TCC)  hycB ← hydrogenase 3, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,849,6740CG 6.0% 219.5 / 4.3 67C59S (TGC→TCC) hycBhydrogenase 3, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (2/2);  total (32/35)
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.34e-01

CGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATA  >  NC_000913/2849526‑2849814
                                                                                                                                                    |                                                                                                                                            
cGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCgcgc                                                                                                                                              <  6:231824/149‑1 (MQ=255)
   aaTATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAgg                                                                                                                                           <  6:82467/149‑1 (MQ=255)
       tCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCCCGGCCCCATCGACGCGGGTGATGGCGTTACCCGGGCAGACCACCGCGCAGGGTg                                                                                                                                        >  1:237976/1‑148 (MQ=255)
         cAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCAt                                                                                                                                     >  1:199669/1‑149 (MQ=255)
         cAGCGGACGGCGGCCGGAAAAGTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCAt                                                                                                                                     <  6:128267/149‑1 (MQ=255)
               aCGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTca                                                                                                                                >  4:154077/1‑148 (MQ=255)
                cGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGGGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTcaca                                                                                                                              <  2:167578/149‑1 (MQ=255)
                   cTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAgtg                                                                                                                           <  1:112392/149‑1 (MQ=255)
                      ccGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAgtgg                                                                                                                          <  5:115693/147‑1 (MQ=255)
                       cGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGa                                                                                                                       >  3:166903/1‑149 (MQ=255)
                          aaaaTTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGa                                                                                                                       >  6:172121/1‑146 (MQ=255)
                               tCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTg                                                                                                               >  8:136666/1‑149 (MQ=255)
                               tCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGAGCATCTTCACAGTGGTGACAGAGCTg                                                                                                               >  1:150669/1‑149 (MQ=255)
                               tCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTACACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTg                                                                                                               >  7:43212/1‑149 (MQ=255)
                                 aaTTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTgct                                                                                                             >  6:13784/1‑148 (MQ=255)
                                  aTTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTgcgg                                                                                                            >  3:167090/1‑149 (MQ=255)
                                     gcgcCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTgcggcg                                                                                                          >  8:3031/1‑148 (MQ=255)
                                         cAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGAt                                                                                                     >  7:71038/1‑149 (MQ=255)
                                         cAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGAt                                                                                                     >  3:132523/1‑149 (MQ=255)
                                               ggCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCtttt                                                                                               <  1:2763/149‑1 (MQ=255)
                                                  aGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTttca                                                                                            >  2:73436/1‑149 (MQ=255)
                                                  aGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTttca                                                                                            <  3:233046/149‑1 (MQ=255)
                                                    gCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACTGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGATCTGCGGCGCAGATTCTTTttcatt                                                                                          <  4:166903/149‑1 (MQ=255)
                                                    gCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTttcatt                                                                                          <  4:167090/149‑1 (MQ=255)
                                                    gCGATGCCGCACAGCTAGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATAGACGCGGGTGATGGCGTTAACCGGGAAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCTGATACTTTtacatt                                                                                          >  4:124738/1‑149 (MQ=255)
                                                       aTGCCGCACAGCTTGTAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTcaca                                                                                                                              <  2:118233/110‑1 (MQ=255)
                                                       aTGCCGCACAGCTTGTAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTcaca                                                                                                                              >  1:118233/1‑110 (MQ=255)
                                                       aTGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAg                                                                                       >  5:241069/1‑149 (MQ=255)
                                                           cGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCAt                                                                                    >  2:99553/1‑148 (MQ=255)
                                                                  cTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCActct                                                                              >  1:132547/1‑147 (MQ=255)
                                                                    tGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCCCAGTGGTGACAGAGCTGCGGCGCAGCTTCGTTTTCATTCAGCATCACTCTCAgg                                                                          >  3:185652/1‑149 (MQ=255)
                                                                          ttACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCgg                                                                     >  1:73994/1‑148 (MQ=255)
                                                                            aCGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCa                                                                   <  1:73436/148‑1 (MQ=255)
                                                                                 caGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCattgatt                                                             <  1:99553/149‑1 (MQ=255)
                                                                                  aGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCattgattg                                                            >  4:214206/1‑149 (MQ=255)
                                                                                  aGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCattgattg                                                            <  3:158269/149‑1 (MQ=255)
                                                                                   ggatttCATTCAACTGCACGGCCCCATCGACGCGGGTTATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTATTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGc                                                           <  5:167333/146‑1 (MQ=255)
                                                                                   ggCTTTCATTCAACTGCACGGCCCGATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGc                                                           <  5:170745/149‑1 (MQ=255)
                                                                                     cTTTCATGCAACTGCCCGGCCCGATCGCCGCGGGTGATGGCGTTAACGGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGTTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTGTCAGGCGCGGCATTGATTGCAg                                                         <  2:172591/149‑1 (MQ=255)
                                                                                          attcaACTGCACTGCCCCATCGTCGCTGGTGATGGCGTTAACCGGGCAGACGACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCg                                                     <  6:241069/148‑1 (MQ=255)
                                                                                          attcaACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGCTTCCAGGCCg                                                     <  6:179773/148‑1 (MQ=255)
                                                                                          attcaACTGCACGGCCCCATCGACGCGGGGGATGGCGGTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGt                                                    <  1:108350/149‑1 (MQ=255)
                                                                                                   cACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGAt                                           <  5:13784/149‑1 (MQ=255)
                                                                                                   cACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGAt                                           <  6:73633/149‑1 (MQ=255)
                                                                                                       gCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATGCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGt                                       <  7:3031/149‑1 (MQ=255)
                                                                                                         cccATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGtct                                     <  3:214206/149‑1 (MQ=255)
                                                                                                          ccATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGtctc                                    >  1:143607/1‑149 (MQ=255)
                                                                                                           cATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTctct                                   <  2:199669/149‑1 (MQ=255)
                                                                                                               gACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAAc                               >  5:112611/1‑149 (MQ=255)
                                                                                                                 cgcgGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAg                             <  8:33428/149‑1 (MQ=255)
                                                                                                                  gcgGGTGCGGGCGTTAACCGGGCAGCCCACCGCGGAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAgg                            <  6:10117/149‑1 (MQ=255)
                                                                                                                  gcgGGTGATGGCGTTAACCGGGCAGACCACCGCGGAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCCTCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAgg                            >  5:10117/1‑149 (MQ=255)
                                                                                                                  gcgGGTGATGGCGTTAACCGGGCAGACCACCGCGGAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGAGTCTCTGAACAgg                            >  5:7030/1‑149 (MQ=255)
                                                                                                                  gaggttGATGGCGTTAACCGGGCAGACCACCGCGGATGGTGCATCTTCACAGTGCTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCACCATCACTCTCAGGCGCGTCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAgg                            <  6:7030/144‑1 (MQ=255)
                                                                                                                      gggATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAggcggc                        <  7:136666/147‑1 (MQ=255)
                                                                                                                       tGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGcc                       >  1:73403/1‑149 (MQ=255)
                                                                                                                            gCGTTAACCGGGCAGACCACCGCGCAGGGGGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTcac                   <  3:124738/148‑1 (MQ=255)
                                                                                                                                  aCCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtgt             <  4:132523/148‑1 (MQ=255)
                                                                                                                                  aCCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAgtgt             <  5:172121/148‑1 (MQ=255)
                                                                                                                                  aCCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGtgtg            <  8:71038/149‑1 (MQ=255)
                                                                                                                                   ccGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAgg                                                        <  6:175907/104‑1 (MQ=255)
                                                                                                                                   ccGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAgg                                                        >  5:175907/1‑104 (MQ=255)
                                                                                                                                   ccGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTgg           >  2:247286/1‑149 (MQ=255)
                                                                                                                                     gggCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTgg           >  2:123120/1‑147 (MQ=255)
                                                                                                                                       gCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGc       <  6:112611/149‑1 (MQ=255)
                                                                                                                                             caccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCt                                                 <  7:218164/101‑1 (MQ=255)
                                                                                                                                             caccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCt                                                 >  8:218164/1‑101 (MQ=255)
                                                                                                                                             caccGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATa  >  5:193598/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                            
CGGAATATCCAGCGGACGGCTGCCGGAAAATTCAATTGCGCCAAACGGGCAGGCGATGCCGCACAGCTTGCAGCTTACGCACAGGCTTTCATTCAACTGCACGGCCCCATCGACGCGGGTGATGGCGTTAACCGGGCAGACCACCGCGCAGGGTGCATCTTCACAGTGGTGACAGAGCTGCGGCGCAGATTCTTTTTCATTCAGCATCACTCTCAGGCGCGGCATTGATTGCAGGCCGTGCTGGCGATGCGTCTCTGAACAGGCGGCCTCACAAGTGTGGCAGCCGATA  >  NC_000913/2849526‑2849814

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