Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,291,339 T→G I288M (ATT→ATG nrfD → formate‑dependent nitrite reductase, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,291,3390TG98.3% 201.7 / NA 60I288M (ATT→ATGnrfDformate‑dependent nitrite reductase, membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base G (32/27);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 4.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACACCGTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCAT  >  NC_000913/4291125‑4291564
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 
tACACCGGTTGGAAATCCCTGGGGTATGGGTTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               <  2:688562/235‑1 (MQ=255)
      gTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGc                                                                                                                                                                                                                              <  1:670105/214‑1 (MQ=255)
      gTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGc                                                                                                                                                                                                                              >  2:670105/1‑214 (MQ=255)
      gTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               >  1:568691/1‑229 (MQ=255)
      gTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               <  2:568691/229‑1 (MQ=255)
      gTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGtgtg                                                                                                                                                                                                             >  1:268639/1‑231 (MQ=255)
        aTGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGTTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCTTGCCTGCCGTGCTGGCGTCCTGTGGGGc                                                                                                                                                                                                         >  2:313717/1‑233 (MQ=255)
         tGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTgggg                                                                                                                                                                                                          <  2:589803/231‑1 (MQ=255)
               tCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               <  1:233774/220‑1 (MQ=255)
               tCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               >  2:233774/1‑220 (MQ=255)
                 ccTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTgg                                                                                                                                                                                                 >  1:436489/1‑232 (MQ=255)
                 ccTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCG‑‑GGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATGGCAGTTCCGGCATGCCGTCCGTGCTGGCGTGGTGTGGGGCCAGTCTGATTGGCg                                                                                                                                                                                           >  2:635082/1‑236 (MQ=255)
                   tgtgGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTACTGGCGGCGTGTGGGGCCAGTCTGGTCGGCg                                                                                                                                                                                           >  1:537883/1‑236 (MQ=255)
                                       tCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGGGTTGATGCTGCGCctttt                                                                                                                                                                        >  1:126546/1‑235 (MQ=255)
                                           gctgGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGGGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGGGCctt                                                                                                                                                                          >  2:381901/1‑229 (MQ=255)
                                            ctgGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGACTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               <  1:696845/191‑1 (MQ=255)
                                            ctgGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCCGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGACTGCCGTGCTGGCGGCgtg                                                                                                                                                                                                               >  2:696845/1‑191 (MQ=255)
                                              ggtggCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCGGCCGTGCTGGCGGCGGGTGGGGCCAGGCTGGTCGGCGGGTTGTGGCTGCGGCTTTTCAtgact                                                                                                                                                                 >  2:92672/1‑231 (MQ=255)
                                                tggCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGGGTTGATGCTGCGCTTTTTCATTCTCTAc                                                                                                                                                              >  1:303783/1‑236 (MQ=255)
                                                                tCTGGCGGTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACgg                                                                                                                                              <  1:92672/236‑1 (MQ=255)
                                                                        tGGGCGATGACGGTAAAGTTCGTGGGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTTGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTa                                                                                                                                       <  1:326497/235‑1 (MQ=255)
                                                                          ggCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCg                                                                                                                                         <  2:104978/231‑1 (MQ=255)
                                                                          ggCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCg                                                                                                                                         >  1:104978/1‑231 (MQ=255)
                                                                          ggCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGGTGATGCTTCGCTTTTTCATTCTCTACGCCGGCCAGGTGACGGGGGCGTATGc                                                                                                                                    >  2:21484/1‑236 (MQ=255)
                                                                          ggCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGCGTGGGGCCAGTCGGGTCGGCGTGTTGATGCTGCGCTTTTTCATGCTCTACCCCCGACAGTTGCCgtgtgtgtaagc                                                                                                                                    >  2:377677/1‑225 (MQ=255)
                                                                           gCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGc                                                                                                                                    <  1:431877/235‑1 (MQ=255)
                                                                                    gTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTggcgg                                                                                                                                                                                                                   <  2:406542/147‑1 (MQ=255)
                                                                                    gTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTggcgg                                                                                                                                                                                                                   >  1:406542/1‑147 (MQ=255)
                                                                                      aaaGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAgg                                                                                                                         <  2:436489/235‑1 (MQ=255)
                                                                                            cGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGCAAAGAGGTGGTTTc                                                                                                                  >  1:4437/1‑236 (MQ=255)
                                                                                            cGTGCGCGGGTGGCGGCATTAGGGGGCGGTTTCTGGCCGTGGTGTTTGTGGCTTGGTTTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTACTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTc                                                                                                                  <  2:537883/236‑1 (MQ=255)
                                                                                              tGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCt                                                                                                                 <  1:518257/235‑1 (MQ=255)
                                                                                                                       ggTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGtt                                                                                       >  2:12579/1‑236 (MQ=255)
                                                                                                                                 cGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCgtt                                                                             <  1:293899/236‑1 (MQ=255)
                                                                                                                                   tggtggTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTa                                                                                                          >  1:102061/1‑205 (MQ=255)
                                                                                                                                   tggtggTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTa                                                                                                          <  2:102061/205‑1 (MQ=255)
                                                                                                                                   tggtggTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCgttg                                                                            <  1:381901/235‑1 (MQ=255)
                                                                                                                                    ggtggtTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGGGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCgttgtt                                                                          >  1:8828/1‑236 (MQ=255)
                                                                                                                                     gtggtTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGGTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTAGTCCTTACTGAAGTCGGTGTTCTTGGgttgttg                                                                         >  2:262504/1‑236 (MQ=255)
                                                                                                                                        gtTCTGGCTTGGGGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCgttgttgtt                                                                       <  1:313717/235‑1 (MQ=255)
                                                                                                                                              gCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGtctc                                                                >  1:111903/1‑236 (MQ=255)
                                                                                                                                              gCTTGGTGTCGCAGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGtct                                                                 <  2:211997/235‑1 (MQ=255)
                                                                                                                                               cTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCg                                                               <  1:618878/236‑1 (MQ=255)
                                                                                                                                                   gtttCGCCGGGCTGGGCCTTATTGTGCCAATGCTGCTCAAGCCGGGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCgggg                                                            <  1:421436/232‑1 (MQ=255)
                                                                                                                                                    tgtCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTc                                                          <  1:49054/236‑1 (MQ=255)
                                                                                                                                                     gtCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCa                                                         <  1:17125/236‑1 (MQ=255)
                                                                                                                                                        gCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGCAGTCGGTTTTCTGGCGTTGTTGGTAAGTCTCGTGGTCAACg                                                      >  2:259574/1‑236 (MQ=255)
                                                                                                                                                                         tgtgCCAATGCTGCGCAATGCGTGGGTCAATCGCAGTGCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACg                                     <  1:721833/236‑1 (MQ=255)
                                                                                                                                                                            gCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTTACCCCGTTGACGGcc                                  >  2:685131/1‑236 (MQ=255)
                                                                                                                                                                            gCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGc                                   <  1:259574/235‑1 (MQ=255)
                                                                                                                                                                                     gctCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGTCGTAAGCCAGAAAAGGGGTGGTTTCTGGACGGATTCCTTCCTGAAGTCGGTTTTCTGTAGTTGTTGCTAAGGCCAGGGGTCAACGTGTTGAACCCGTTGTCGGCCTTcgatgc                          >  2:392365/1‑231 (MQ=255)
                                                                                                                                                                                      ctCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTAtt                                                                                                        <  2:333347/156‑1 (MQ=255)
                                                                                                                                                                                      ctCAAACCGTGGGTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTAtt                                                                                                        >  1:333347/1‑156 (MQ=255)
                                                                                                                                                                                                gggTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCTGGAGTGCGGt                  >  1:72999/1‑232 (MQ=255)
                                                                                                                                                                                                gggTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTGGTTGTTAAGTCTCGGGGTAAACGCGTTGAACCCTTTTAAGGACTTCGCGTGTGTGCggg                  >  2:518650/1‑231 (MQ=255)
                                                                                                                                                                                                gggTCAATCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGCCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGt                  >  1:417358/1‑232 (MQ=255)
                                                                                                                                                                                                      gtCGCAGTCCCGTCATGCCTGCCGTGCCGGCGGCGTGTGGGCCCAGTATGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACATTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGc         <  2:8828/234‑1 (MQ=255)
                                                                                                                                                                                                      aTCGCAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTgcgc             >  1:531786/1‑231 (MQ=255)
                                                                                                                                                                                                          cAGTTCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGcc    <  2:268639/236‑1 (MQ=255)
                                                                                                                                                                                                             ttCCGGCATGCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCat  >  1:135808/1‑235 (MQ=255)
                                                                                                                                                                                                                      |                                                                                                                                                                                                                                 
TACACCGTATGGAAATCCCTGTGGTATGGGGTGAAATCTTCCTGCTGGTGGCGTTTTTTGTCGGTCTGGCGCTGGGCGATGACGGTAAAGTTCGTGCGCTGGTGGCGGCATTAGGTGGCGGTTTCTGGACGTGGTGGTTCTGGCTTGGTGTCGCCGGGCTGGGGCTGATTGTGCCAATGCTGCTCAAACCGTGGGTCAATCGCAGTTCCGGCATTCCTGCCGTGCTGGCGGCGTGTGGGGCCAGTCTGGTCGGCGTGTTGATGCTGCGCTTTTTCATTCTCTACGCCGGACAGTTGACGGTGGCGTAAGCCAGAAAAGAGGTGGTTTCTGGACGTATTCCTTCCTGAAGTCGGTTTTCTGGCGTTGTTGTTAAGTCTCGGGGTCAACGTGTTGACCCCGTTGACGGCCTTCGCGGGAGTGCGGTTGCGCTGGCCTGCCAT  >  NC_000913/4291125‑4291564

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