Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,606,887 A→G 21.0% M439V (ATG→GTG)  hyfF → hydrogenase 4 subunit F

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0077792,606,8870AG21.0% 129.9 / 7.9 79M439V (ATG→GTG) hyfFhydrogenase 4 subunit F
Reads supporting (aligned to +/- strand):  ref base A (23/36);  new base G (15/2);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 6.15e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCA  >  NC_007779/2606657‑2607123
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                        <  1:193696/239‑1 (MQ=255)
tGCGCAAGTACGGCACGCGGGATCTCAACGTCGTCGGTGGGATGCTCAAAATCATGCCATTTCCCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                        <  1:923774/239‑1 (MQ=255)
  cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGGCGGGTGGCTGcgcg                                                                                                                                                                                                                                      >  1:89156/1‑239 (MQ=255)
       gTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCCGGGACTGGGGGGTGCGCgggggg                                                                                                                                                                                                                                 >  1:2067887/1‑236 (MQ=255)
          cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTGTTAGCGAATTTATGACCATTACCGCCGGACTGGAACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCGGGTGCTGGCGGGTTAGGGACGGCTtgctgccgggggggtta                                                                                                                                                                                                                               >  1:1286741/1‑228 (MQ=255)
          cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTaa                                                                                                                                                                                                                              >  1:1670027/1‑239 (MQ=255)
          cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGGTaa                                                                                                                                                                                                                              >  1:936655/1‑239 (MQ=255)
            gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCCGGGAAGGGTGGCTGCGCGGGTGTTAATg                                                                                                                                                                                                                            >  1:535345/1‑239 (MQ=255)
            gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGCGCTGGCGGGGCCGGGTCCGATGGGCGCGCCGgtgtgt                                                                                                                                                                                                                               >  2:186656/1‑234 (MQ=255)
                 cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGGAAGGGTTGGTTCGCGGGGGGTTATGGTGaa                                                                                                                                                                                                                       >  1:1067878/1‑239 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTAAGGGTGGCCGGGCGGGGGTTAATGGCGGAAACGGCCCAGGGCgg                                                                                                                                                                                                       >  1:884602/1‑238 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGGTGGCTGCGCGGGGGGTAATGGCGAAACCGCCGCCGGGCGt                                                                                                                                                                                                       >  1:327838/1‑239 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGTTAAAGGCGCAACCGCCgccggccgg                                                                                                                                                                                                       >  1:816277/1‑238 (MQ=255)
                                gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTACGGGTGGGTGCGCGGGGGGTAATGGCGGAACCGCCGCCGGCCGt                                                                                                                                                                                                       >  1:843706/1‑239 (MQ=255)
                                 tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATgg                                                                                                                                                                                                                                              <  1:1658054/201‑1 (MQ=255)
                                 tcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCTGTACGGGTGGCTGGGCGGGGGTTAATTGCGAAACCCCCCCCGGCCggt                                                                                                                                                                                                      >  1:520756/1‑237 (MQ=255)
                                      gggATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTTATGGGGGAACCGCCGCAGGCCGGTAAACg                                                                                                                                                                                                 >  1:919190/1‑239 (MQ=255)
                                          tGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGGTAATGGCGGAACCGCCGCAGGCCGGTAAACggggg                                                                                                                                                                                             >  1:498181/1‑238 (MQ=255)
                                           gCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTAAGGGTGGGTGCGCGGGGGGTAATGGCGAAACCGCCGCCGGCCCGTTACCCggggg                                                                                                                                                                                             >  1:1174737/1‑238 (MQ=255)
                                            cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGGGTTAATTGGGGAAACGCCGCAGGCCGGTAAACGGGGGGa                                                                                                                                                                                           >  1:609466/1‑239 (MQ=255)
                                            cTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTACGGGTGGGTGCGCGGGGGTTAATGGGGAAAACGCCGCCGGCCGGTTACCggggggg                                                                                                                                                                                           >  1:1820528/1‑236 (MQ=255)
                                                   tCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAt                                                                                                                                                                                                                                                <  1:1819495/181‑1 (MQ=255)
                                                      tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGGCGTTAACCGGGGGGGTCTCGGCTgg                                                                                                                                                                                >  1:1996242/1‑239 (MQ=255)
                                                           tttACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTAAGGGTGGGCTCGCGGGGGGTTAAGGCGGAACCCCCCCAGGGCCTTTACCCGGGGGGGCTCCGCTGGGTgg                                                                                                                                                                            >  1:1774715/1‑237 (MQ=255)
                                                                  ccgTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTACGGGTTGCTGCGCGGGGGTTAATGGGGGAACCCCCGCAAGCCGTTAACCCGGGGGGGCTCGGCTGGGTTaccacccc                                                                                                                                                                     >  1:1340566/1‑236 (MQ=255)
                                                                     tGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAAGGCTGCGCGGGTGGTAATGGGGGAAACCCCGCCGGGCCGTAACCGGGGGGGTCTCGGCGGGGTGGACCACCTcccc                                                                                                                                                                  >  1:434906/1‑239 (MQ=255)
                                                                       ttGTTGGGGGGGGGTGGGCTGGCGCTGGCGGGGACCCCGCCCTTCAACTTTTTTTTGGGCTAATTTATGACCATTCCCGCCGGACTGGCACGTATTCCCCTGCTGATTAGCGTCCTGCTGTTTTTGCTGTTAACGGTGGGGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAAt                                                                                                                                                               <  2:1996242/238‑1 (MQ=255)
                                                                            tGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCCCCCTTCAACATTTTTCTTAGCGATTTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGa                                                                                                                                                          <  1:234146/239‑1 (MQ=255)
                                                                                               ttggcGGGGATGCCGCCCTTAACCATTTTTCTTAGCGATTTTATGCCCTTTCCCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGac                                                                                                                                                                           <  1:1472706/202‑1 (MQ=255)
                                                                                                  gcAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAt                                                                                                                                                                                                                                                <  1:42317/134‑1 (MQ=255)
                                                                                                            cgccCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                        >  1:268315/1‑131 (MQ=255)
                                                                                                            cgccCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                        <  2:268315/131‑1 (MQ=255)
                                                                                                              cccTTCAACATTTTTTTTAGCAATTTTATGCCCTTTCCCGCCGGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                        <  1:996459/239‑1 (MQ=255)
                                                                                                                   aaataTTTTTCTGAGCGAATTTAGGACTGTCCCCGCCGGACGGGCACGTAATCCCCTGCTGTTTGTCGTCCGGCTGTTTTTGCTGTAACCGCTGGTGCGGGCGGGCGTGGTACGGATGGCTGGGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTc                                                                                                                   <  2:936655/235‑1 (MQ=255)
                                                                                                                         ttttCTTAGCGATTTTATGACCATTCCCGCCGGCCTGGCACGTAATCCCCTGCTGATTATCGTCCTGTTGTTATTGCTGTTAACGCTGGTGCGGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAAc                                                                                                             <  1:1922538/239‑1 (MQ=255)
                                                                                                                         ttttCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGGTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCCCGCCAATGGGGATTCTGCTGGTCATGATGCTGGCGAGGGGGACGCCTATTCCCCaa                                                                                                              >  1:1407625/1‑238 (MQ=255)
                                                                                                                          tttCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGTCATGGTGCTGGCGATGGGGACGCATATTCCACAAc                                                                                                             >  1:105792/1‑238 (MQ=255)
                                                                                                                                aGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCCGtca                                                                                                       >  1:223587/1‑238 (MQ=255)
                                                                                                                                aGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAAACGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCTCGCCCATGGTGAtt                                                                                                                                                        >  1:1532281/1‑189 (MQ=255)
                                                                                                                                            ccctttACCGCCGGACGGGCCCGTAACCCCCTGCTGATTTTCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTCCGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTgg                                                                                           <  1:708915/234‑1 (MQ=255)
                                                                                                                                               aTTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGGTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTggcggg                                                                                       >  1:433073/1‑239 (MQ=255)
                                                                                                                                                 tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGGTGACCACCTCGCCCATGGTGATTCTGGTGGTCATGATGCTGGCGGTGGGAACGCATATTCCAAAAACTGTCATCAGGATCCTggcgggcg                                                                                     >  1:1304054/1‑239 (MQ=255)
                                                                                                                                                  cccgccgGCCTGGCCCGTAATCCCCTGCTGTTTATCGTCCTGCTGTTATGGTTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                                                    <  1:648139/238‑1 (MQ=255)
                                                                                                                                                  ccccccgGACTGGCGCGTAATCACCTGCTGATTATGGTCCTGCTGTTATTGCTTTTAACGCTGGTGCTGGCGGGCCTGGTGCGGATGGCTGCGCGGGTGTAAAGGGGAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGgcgc                                                                                    <  1:808651/235‑1 (MQ=255)
                                                                                                                                                    cgccgGACTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTACCGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                        <  1:1243113/201‑1 (MQ=255)
                                                                                                                                                    cgccgGACGGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTACCGCGGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCtt                                                                                  <  1:276278/239‑1 (MQ=255)
                                                                                                                                                       cgGCCTGGCACGTAATCCCCTGCTGATTATCGTCTTGCTTTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCa                                                                               <  1:1904590/239‑1 (MQ=255)
                                                                                                                                                       cgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGACTGGTACGGGTGGCTGCGCGGGTGTTAATGGCGAAAccg                                                                                                                                                                                                                   >  1:2028037/1‑109 (MQ=255)
                                                                                                                                                          ccTGGCACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCGGTTCCGGAGGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTgg                                                                                                                                                                                <  1:272808/138‑1 (MQ=255)
                                                                                                                                                          aCTGGCACGTAACCCCCTGCGGATTATTGTCCTGCTGTTATTGTTGTTAACGCTGGGCTGGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACt                                                                             <  1:1032935/238‑1 (MQ=255)
                                                                                                                                                               cACGTAATCCCCTGCTGTTTATCGTCCTGTTGTTATTGCGGTTAACGCGGGTGCGGGCGGGCCTGGTACGGAGGGCTGCGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                        <  1:719238/238‑1 (MQ=255)
                                                                                                                                                               cACGTAATCCCCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                        <  1:186656/238‑1 (MQ=255)
                                                                                                                                                                  ttAATCACCGGCTGATTATCGTCCTGCTGTTAGTGTTTTTAGCGCTGGTGCTGGCGGGCCTGGTACGGAGGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTc                                                                       <  1:868088/235‑1 (MQ=255)
                                                                                                                                                                  gTACTCCCCTGCTGATTATCGTCCGGCTGTTTTTGCTGTTAACGCTGGTGCGGGCGGGCCGGGTAGGGATGGCTGCGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                 <  1:1107348/226‑1 (MQ=255)
                                                                                                                                                                           tgctgATTATCGTCCTGCTGTTATTGCTGTTACCGCTGGTGCGGGCGGGCCTGGTACGGATGGCTGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACg                                                           <  1:157901/239‑1 (MQ=255)
                                                                                                                                                                            ggtgATTATCGCCCTGCTGTTATTGCTGTTAACGCGGGTGCTGGGGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACg                                                           <  1:1600378/236‑1 (MQ=255)
                                                                                                                                                                            gctgATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                        >  1:1495550/1‑177 (MQ=255)
                                                                                                                                                                             ctgTTTATCGTCCTGTTGTTATTGCTGTTACCGCTGGTGCTGGCGGGCCTGGTACGGAGGGCTGCGCGGTTGTTAAGGTCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGa                                                                                                                                                                            <  1:1154253/124‑1 (MQ=255)
                                                                                                                                                                             ctgATTATCGTCCTGCTTTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTctg                                                                                                                                                     <  1:994907/147‑1 (MQ=255)
                                                                                                                                                                              tgATTACCTTCCGGTTGTTATTGCTGTTAAGGCTGGAGATGGGGGGCGTGGTACGGGTGGGTGCGTGGGTGTTAATGGGGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCa                                                         <  1:697885/238‑1 (MQ=255)
                                                                                                                                                                                aTTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCTAATGGGGAtc                                                                                                                                                        >  1:1794868/1‑140 (MQ=255)
                                                                                                                                                                                      gTCCTTCTTTTTTTGCTGTACAGGCTGGTGGTGGGGCGCGGGGTACGGGTGGCTGCGGGGGTGTTAAGGGCAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                                                 <  2:1407625/238‑1 (MQ=255)
                                                                                                                                                                                       tCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGGGATTCTGCTGGGCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCCTCCACTATAGTCCTCTCAggg                                                              >  1:1245925/1‑224 (MQ=255)
                                                                                                                                                                                          tgctgTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTCCGGATGGCTGCGGGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                            <  1:2109449/239‑1 (MQ=255)
                                                                                                                                                                                             ttttATGGCTGTTAGCGCTGGTGCGGGGGGGCCTGGTACGGATGGCGGCGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgcctg                                              <  1:737926/232‑1 (MQ=255)
                                                                                                                                                                                                tATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                        <  1:210278/157‑1 (MQ=255)
                                                                                                                                                                                                tATTGCTGTTAACGCTGGTGCGGGGGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACg                                       <  1:792341/238‑1 (MQ=255)
                                                                                                                                                                                                     cTGTTAGCGCGGGTGCGGGCGGGCTGGGTAGGGAGGGCTGCGCGGGTGTAAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCa                                  <  1:1877243/238‑1 (MQ=255)
                                                                                                                                                                                                      ttttAACGCTGGTGCGGGCGGGCCTGGTACGGAGGGCTGCGCGGGTTTTAATGGCAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATAGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCa                                                                               <  1:238593/190‑1 (MQ=255)
                                                                                                                                                                                                      tGTTAACGCTGGGGCCCGGCGGGCTTGGTCGGCTGGCGGCGGCTGGTTGTAATGCGAAAACCGCGGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGCCCACCTCGCCAATGGTGATTCTGCTGGTCAGGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGAGACGCACGATCTGCCTGCACAACGTAGCa                                  <  2:834642/238‑1 (MQ=255)
                                                                                                                                                                                                         tAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATATGCCTGCACAACGTAGCAcccg                              >  1:1292671/1‑236 (MQ=255)
                                                                                                                                                                                                                 gTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGGGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCCCTATAGTCCTCTCAGGGACGCACGATctgcctg                                              >  1:416248/1‑214 (MQ=255)
                                                                                                                                                                                                                                   tACGGATGGGTGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                            <  1:1335915/200‑1 (MQ=255)
                                                                                                                                                                                                                                   tACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCAcc     >  1:318574/1‑239 (MQ=255)
                                                                                                                                                                                                                                    ccGGTTGGCTGGGCGGGTTTTAAGGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCAcc     <  1:308627/237‑1 (MQ=255)
                                                                                                                                                                                                                                       gATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCat                                                                                                                                             >  1:1133654/1‑99 (MQ=255)
                                                                                                                                                                                                                                       gATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                 >  2:1076172/1‑159 (MQ=255)
                                                                                                                                                                                                                                       gATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                 <  1:1076172/159‑1 (MQ=255)
                                                                                                                                                                                                                                       gATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAg                                    >  1:30277/1‑204 (MQ=255)
                                                                                                                                                                                                                                        aTGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTTCACCATAGGCCTCCCAGGGACGCACCACCTGCCTGCACAACGGAGCACCTGTCATGCTTTTTTTCCTTCACGCCCCGCa  >  2:810469/1‑237 (MQ=255)
                                                                                                                                                                                                                                        aTGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACt                                                                             >  1:853483/1‑162 (MQ=255)
                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                               
TGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCA  >  NC_007779/2606657‑2607123

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