Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,605,247 G→A 5.8% G146D (GGC→GAC)  hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,605,2470GA 5.8% 226.2 / 2.7 69G146D (GGC→GAC) hyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base A (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.61e-01

GGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCT  >  NC_000913/2605102‑2605394
                                                                                                                                                 |                                                                                                                                                   
ggTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAggcg                                                                                                                                                   <  5:46830/148‑1 (MQ=255)
     tCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGgctgc                                                                                                                                              >  1:97060/1‑148 (MQ=255)
         aCGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTcc                                                                                                                                           <  7:114496/147‑1 (MQ=255)
             aTAAAAATCGCGAACACCCGCATACCGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGc                                                                                                                                      >  1:260414/1‑148 (MQ=255)
                      gcgAACACCCGCATATCGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCGCGGGCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAg                                                                                                                             <  2:260414/148‑1 (MQ=255)
                          acacCCGCATAACGGCACGAATCGTTATTAGGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGttat                                                                                                                         >  6:251028/1‑148 (MQ=255)
                                    aaCGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGc                                                                                                                >  7:216990/1‑147 (MQ=255)
                                         cACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATa                                                                                                            >  2:17962/1‑146 (MQ=255)
                                              aTCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGTTCCTTGTACTCTCCGCCACGCTGATCTGTCAGTGGTTGAAGTTTGCATTTTCAGGCGGCTGATGCTGGGCCTTGTGCAGATATTCCAAGAGCGATACAgcgcc                                                                                                     >  6:70062/1‑147 (MQ=255)
                                                cGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGc                                                                                                   >  1:198516/1‑148 (MQ=255)
                                                  gtattaCGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGt                                                                                                 >  4:171338/2‑148 (MQ=255)
                                                     ttaCGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTc                                                                                               <  2:198516/147‑1 (MQ=255)
                                                        cGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTc                                                                                               >  5:287994/1‑144 (MQ=255)
                                                        cGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTc                                                                                               <  6:287994/144‑1 (MQ=255)
                                                        cGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCa                                                                                           >  1:19180/1‑148 (MQ=255)
                                                            tttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTaa                                                                                       <  5:251028/148‑1 (MQ=255)
                                                            tttttACTGGGGTTTATCGGCGCGATGGCGGGCCTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTaa                                                                                       <  2:97060/148‑1 (MQ=255)
                                                             ttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGt                                                                                                                                    <  3:122523/102‑1 (MQ=255)
                                                             ttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGt                                                                                                                                    >  4:122523/1‑102 (MQ=255)
                                                             ttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGtt                                                                                                                           >  2:373775/1‑111 (MQ=255)
                                                             ttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGtt                                                                                                                           <  1:373775/111‑1 (MQ=255)
                                                             ttttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTaa                                                                                       >  3:14129/1‑147 (MQ=255)
                                                              tttACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTaa                                                                                       >  7:7585/1‑146 (MQ=255)
                                                                   tgttgtTTATCGTCGCGTTTGCGTGACTGTTACTCTCCTTGACGCTGCTCGGTAAGTTGTTGTTTTTTGAAATTTCAGGCGGCTGCTACTGGGCGTTGTTCAGTTATTACCTGAGCGATAAAGCGGAGCGTTCTGCACTAAAAGCGtt                                                                                <  4:14129/145‑1 (MQ=255)
                                                                    ggTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTa                                                                               >  4:312173/1‑148 (MQ=255)
                                                                     gtgtTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGtt                                                                                <  5:108528/146‑1 (MQ=255)
                                                                     gtgtTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTAc                                                                              <  8:7585/148‑1 (MQ=255)
                                                                       gtTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGtt                                                                                >  6:251893/1‑144 (MQ=255)
                                                                       gtTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGtt                                                                                <  5:251893/144‑1 (MQ=255)
                                                                         ttATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGCTAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTAt                                                                          >  6:368829/1‑148 (MQ=255)
                                                                         ttATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTc                                                                                               >  3:284082/1‑127 (MQ=255)
                                                                         ttATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTc                                                                                               <  4:284082/127‑1 (MQ=255)
                                                                           aTCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGtt                                                                                >  2:70944/1‑140 (MQ=255)
                                                                           aTCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGtt                                                                                <  1:70944/140‑1 (MQ=255)
                                                                           aTCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATc                                                                         <  1:17962/147‑1 (MQ=255)
                                                                              ggCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTc                                                                     <  3:312173/148‑1 (MQ=255)
                                                                                  cgATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCatat                                                                 >  5:311708/1‑148 (MQ=255)
                                                                                           gACTGGTACTCTCCTCGACGCTGCTCGGTAAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGtt                                                        >  3:73389/1‑148 (MQ=255)
                                                                                               ggTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCg                                                                                  <  8:347441/118‑1 (MQ=255)
                                                                                               ggTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCg                                                                                  >  7:347441/1‑118 (MQ=255)
                                                                                               ggTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTgggg                                                    >  1:118910/1‑148 (MQ=255)
                                                                                               ggTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTgg                                                      >  2:20397/1‑146 (MQ=255)
                                                                                                gTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTggg                                                     <  1:20397/146‑1 (MQ=255)
                                                                                                     ctcctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCAc                                                                       >  5:339278/1‑123 (MQ=255)
                                                                                                     ctcctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCAc                                                                       <  6:339278/123‑1 (MQ=255)
                                                                                                     ctcctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTc                                                                     >  4:329692/1‑125 (MQ=255)
                                                                                                     ctcctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTc                                                                     <  3:329692/125‑1 (MQ=255)
                                                                                                     ctcctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTAt                                              >  1:197172/1‑148 (MQ=255)
                                                                                                        ctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTAt                                              <  2:197172/145‑1 (MQ=255)
                                                                                                        ctcGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCtt                                           <  2:118910/148‑1 (MQ=255)
                                                                                                          cGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGACGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACt                                                                      <  6:310812/119‑1 (MQ=255)
                                                                                                          cGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGACGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACt                                                                      <  6:305634/119‑1 (MQ=255)
                                                                                                          cGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGACGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACt                                                                      >  5:310812/1‑119 (MQ=255)
                                                                                                          cGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGACGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACt                                                                      >  5:305634/1‑119 (MQ=255)
                                                                                                           gACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTg                                          <  8:216990/146‑1 (MQ=255)
                                                                                                            aCGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTa                                               <  5:232896/140‑1 (MQ=255)
                                                                                                            aCGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTa                                               >  6:232896/1‑140 (MQ=255)
                                                                                                                  ctCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTgccgccgc                                   <  3:171338/146‑1 (MQ=255)
                                                                                                                     ggTCAGTTGTTGTTTTTTGAAACTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCGCGTTGGGGTTGTATCTTGCCGCCGCCACg                               <  8:177433/147‑1 (MQ=255)
                                                                                                                          gttgttgttTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTGGGGGTTGTATCTTGCCGCCGCCACGCTGttt                         <  6:311708/148‑1 (MQ=255)
                                                                                                                            tgttgttTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGttttt                       >  8:321708/1‑148 (MQ=255)
                                                                                                                              ttgttTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGc                     <  7:321708/148‑1 (MQ=255)
                                                                                                                              gtgttgtTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGAGCAGTTATTACCAGAGCGCTAAAGCGCAGCGTTCAGCACTAAAAGCGGTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGc                     <  4:73389/147‑1 (MQ=255)
                                                                                                                                 ttttttGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGa                  >  7:194905/1‑148 (MQ=255)
                                                                                                                                   ttttGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGAcc                >  3:295213/1‑148 (MQ=255)
                                                                                                                                         aaTTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCg                                                          <  7:298846/100‑1 (MQ=255)
                                                                                                                                         aaTTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCg                                                          >  8:298846/1‑100 (MQ=255)
                                                                                                                                          attagcggcggcTTCTACTGGGCGTGGTTCAGTTATTACCAGAGCGATACAGCGCAGCGTTCTGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGATGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGt         <  5:70062/142‑1 (MQ=255)
                                                                                                                                          aTTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGt         <  8:194905/148‑1 (MQ=255)
                                                                                                                                               aGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGc                     <  6:168994/131‑1 (MQ=255)
                                                                                                                                               aGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGc                     >  5:168994/1‑131 (MQ=255)
                                                                                                                                                 gcggcTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCt  <  4:295213/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GGTTATCTGACGGATAAAAATCGCGAACACCCGCATAACGGCACGAATCGTTATTACGCATTTTTACTGGTGTTTATCGGCGCGATGGCGGGACTGGTACTCTCCTCGACGCTGCTCGGTCAGTTGTTGTTTTTTGAAATTACAGGCGGCTGCTCCTGGGCGTTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCT  >  NC_000913/2605102‑2605394

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