Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_007779 2,606,887 A→G 24.7% 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,8870AG24.7% 63.4 / 7.8 47M439V (ATG→GTG) hyfFhydrogenase 4 subunit F
Reads supporting (aligned to +/- strand):  ref base A (17/18);  new base G (10/2);  total (27/20)
Fisher's exact test for biased strand distribution p-value = 4.67e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAT  >  NC_007779/2606653‑2607124
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
gtccTGCTCAGTTACGGCCGGGGCGATCTCATGGTGGTGGGTGGGCTTGTCAAAATCAGCCCATTTACCGCCGTGCGGTTTGGCGGTGGTGCGCTGGCGCTGGCAGGGAGGCCGCCGTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTACTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGAt                                                                                                                                                                                                                                                <  2:589318/233‑1 (MQ=255)
gTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGCATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGACTGGTACGGATggg                                                                                                                                                                                                                                             >  1:1517292/1‑238 (MQ=255)
    tGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTGCCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                         <  1:325776/239‑1 (MQ=255)
    ggataaaGTCCGGCACGCGGGATCTCAACGTCGTCTGTGGGATGCCAAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                         <  2:406146/234‑1 (MQ=255)
     gtccaatgccGGCACGCGCGATCTCACCGTCGTCTGTGGGAGGCTCAAAATCGTGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTgcg                                                                                                                                                                                                                                         <  2:159321/229‑1 (MQ=255)
     gCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTACCGGTGGCTgggc                                                                                                                                                                                                                                        >  1:2266437/1‑239 (MQ=255)
      cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGAGCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCAGGGACCGGATGGCTccgc                                                                                                                                                                                                                                        >  1:9166/1‑239 (MQ=255)
      cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTACGGGTGGCTGcgcg                                                                                                                                                                                                                                       >  1:2103411/1‑239 (MQ=255)
      cTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGcgcg                                                                                                                                                                                                                                       >  1:823515/1‑239 (MQ=255)
           gTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGAAGGATGGCTGCGCgggggg                                                                                                                                                                                                                                  >  1:794172/1‑236 (MQ=255)
              cGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCCCGTAATCACCTGCTTATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTTCGGATGGCTGCGCGGGGGTTat                                                                                                                                                                                                                               >  1:1035337/1‑238 (MQ=255)
                gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGGTGCGCggggggttatg                                                                                                                                                                                                                             >  1:360154/1‑231 (MQ=255)
                gCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGAAGGGTGGCTGGGGGGGGGGTaa                                                                                                                                                                                                                               >  1:159321/1‑237 (MQ=255)
                     cgcgATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCTGTACGGGTGGCTTCGCGGGGGGTTATggggga                                                                                                                                                                                                                        >  1:1573646/1‑235 (MQ=255)
                                    gtcTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCCTCAACATTTTTCTTAGCGAATTTATTACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGTACCGGGGGGCGGGCGGGGGGTTAAGGGGGAAACCCCCCCAGGccc                                                                                                                                                                                                          >  2:1641909/1‑238 (MQ=255)
                                                          tGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGAAGGATGGGTTCGCGGGGGTTTATGGCGAAACCGCCGCCGGCCGGTAACCGGGGTGGTTGCGGCTgg                                                                                                                                                                                   >  1:460270/1‑239 (MQ=255)
                                                               tttACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCCCCGGCCGTTCACCGGGGGGGTCTTGGCCggggga                                                                                                                                                                               >  1:1688654/1‑234 (MQ=255)
                                                                 tACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGGTGGCTGCGGGGGGGTTAATGGCGAAACCGCCGCCGGCCGGTTACCGGGGGGAGCTCGggtgggtggac                                                                                                                                                                             >  1:545388/1‑235 (MQ=255)
                                                                      ccgTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCCGGTACGGATGGCTGGGCGGGTGGTTATGGCGAAACCGCCGCACGGCGTTAACCCGGGGGAGCTCGGCTGGTTGaccccccc                                                                                                                                                                        >  1:531999/1‑236 (MQ=255)
                                                                              tttGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGGCTGGGACGGATGGCTGGGCGGGTGGTAATGGCGGAACCGCCGCCGGGCCGTAAACGGGGGGGTGTCGGGTGGGTGGACACCCCGCCCATgg                                                                                                                                                              >  1:208407/1‑238 (MQ=255)
                                                                                                                  cgagtgaccgtttttCTTGGCTATTTTATGCCCTTGCCCGCGGGCTGGGCAGGTCTTCCCCTGTGGTTTATCGTCCTGCTGTTTTTGCTGTTAAGGTGGTTGCTGGGGGGCCGGTTGCGGAGGGCGGCGCGGTTGTGAATGGCAAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                         <  2:1963302/229‑1 (MQ=255)
                                                                                                                  cccTTCATATTTTTTCTTGGCGTTTTTATGACCTTTCCCGCGGGGTGGGCCCGTATTCCCCTGCTGATTTTCGTCCTGTTGTTGTTGTTGTTAACGGGGGTGCTGGCGGCCCTGGTGGGGAGGGCTGCGGGGGTGTTATGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCa                                                                                                                         <  2:431732/239‑1 (MQ=255)
                                                                                                                    cTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGATGGCTGCGCGGGTGTTAATGGGGCAACCGCCGCCGGGCGTTAACCGGGGGGGTCTCGGCTGGGTTGCCACCCCCCCCCTGGGGAGTCTGTCGGGCATGGTGCTGGGGAGGGGGACGCaca                                                                                                                       >  2:833955/1‑237 (MQ=255)
                                                                                                                                 cTTAGCGAATTTATGACCTTTACCGCGGGCCGGGCGCGTATTCACCTGTTGATTATCGTCCGGTTGTTATTGCTGTTAACTCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGt                                                                                                          <  1:486081/239‑1 (MQ=255)
                                                                                                                                                   aTTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGGTAATGGCGAAAACGCCGCCGGCCGGTAAACGGGGGGATCTCGGGGGGGTGACCCCCCCCCCCAACGGGGTTTTGGCGGGCATGGTTCTTGCGATGGGGACCGATATTACCCAAACCGTCATCAAGACACTgggggg                                                                                        >  2:1471804/1‑239 (MQ=255)
                                                                                                                                                     tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCCGGCCGTTAACCGGGGGGGTCGCGGGTGGGTGGGCACCACGCCCACGGGGATTCTGGCGGGCATGATGGCGGGGGGGGGGACCCACATATCCCAAACTGTCCTCAGGATCCCggcgggcg                                                                                      >  2:54818/1‑239 (MQ=255)
                                                                                                                                                     tACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTCGTTCGGATGGCTTCGCGGGTGTTAATGGCGAAAACGCCGCAGGCCGGTAACCGGGGGGATCTCGGCTGGGTGACCACCCACCCAATGGGTGATCTTGCGGGCATGATGCTGGCGGTGGGGACGCAAATTTCAAAACCTTTCCCCAGGAGCACggcggggg                                                                                      >  2:2044672/1‑237 (MQ=255)
                                                                                                                                                         gccgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGGTAACCGGGGGGGTCTCGGCTGGGTGACCACCCCCCCCATGGTGGATCTGCTGGGCATGATGCTGGGGGGGGGGAAGCCTATTCCCCAACCTGTCCTCCGGGTCCTgggccggtctct                                                                                  >  2:547853/1‑229 (MQ=255)
                                                                                                                                                           cgGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGGGTTAATGGCGGAAACGCCGCAGGGCGCTAACCGGGGGGGGCTCGGCTGGGTGGCCACCCCGCCCAAGGGGGTTCTGCTGGGCATGATGCTGGGGGTGGGGACACATATTCCACAACCCGGCATCAGGATCCTCGCGGgggccttc                                                                                 >  2:757974/1‑234 (MQ=255)
                                                                                                                                                             gACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGGTGGCTGCGCGGGGGTTTATGGCGAAACCGCCCCCGGCCGTTGACCGGGGGGGGCCCGGCTGGGTGACCCCCCCCCCCACGGGGGATCTGCTGGGGCTGATTCTGGCGGTGGGGAACCCTCTTCCACAACCCTTCCTCAAGGTCCTcgcggggcggtccc                                                                                >  2:1608897/1‑229 (MQ=255)
                                                                                                                                                             gACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGAAGGAAGGCTGCGCGGGGGGTAATGGCGCAACCGCCGCAGGCCGTTAAACGGGGGGATCTCGGCTGGGTTCCCAACCCCCCACATGGGTGTCTGCTGGGCATGATGCTGGTGTAGGGGAAGCATAATTCACAACCTGTACTCAGGAGTCTGGGGGGGGCCTCCACt                                                                              >  2:1608907/1‑239 (MQ=255)
                                                                                                                                                                  gCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGGACGGGTGGCTGCGCGGGGGGTAATGGCGGAACCGCCGCCGGGCGCTAACCCGGGGGATCTCGGCTGGGTGGCCCCCCCCCCCAAGGTGGATCTTCCGGGTCTGATGGTGGGGGGGGGGACGCCTATTTCACAACACCTCCTCACGAACCTGGCGGGGGCGTCCCCcataat                                                                         >  2:1369248/1‑237 (MQ=255)
                                                                                                                                                                   cgcttatcccccTGCTGATTTTCGTCTTGTGTTTTTGGGGTTTACGGGGGGTGTGGGCGGGCCTGGTCGGGAGGGCGGCGCGGGTGTAAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                         <  2:1573646/228‑1 (MQ=255)
                                                                                                                                                                   cACGTAATCCCGTGCTGTTTATCGTCTTGCTGTTATTGCTGTTAGCGCGGGTGCTGGCGGGCGTGGTGCGGATGGGTGGGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGt                                                                         <  1:948575/238‑1 (MQ=255)
                                                                                                                                                                   cACGTAACCGCCTGCTGATTATCGGCCTGCTTTTTTTGCTGTTACCGCGGGTGGGGGGGGGCCTGGGAGGGAGGGCTGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGaccac                                                                                                                                                                           <  1:1750191/142‑1 (MQ=255)
                                                                                                                                                                                 ctgTTTATCGTCCTGTTTTTTTGGCTGTTAACGGTGGTGGTGGCGGGCCTGGTACGGATGGCGGCGCGGGTTTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATggg                                                                                                                               <  1:431491/170‑1 (MQ=255)
                                                                                                                                                                                    atatgaaacgggcttTTATTGTTGTTAACGCTGGTGCTGGCGGCCGGGGTGCGGATGCCTGCGGGGTTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGa                                                       <  1:1550082/228‑1 (MQ=255)
                                                                                                                                                                                         cgcccTGCTGTTTTTGCTGTTAACGCTGGTGCGGGCGGGCCTGGTCGGGAGGGCGGGGCGGGTGTTAAGGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                                                  <  1:461512/236‑1 (MQ=255)
                                                                                                                                                                                          gTCCTGGTGTTTTGGCGGTTAGCGGGGGGGCTGGCGGGCCGGGGAGGGGGGGTTGGGCGGGGGTTAATGGGGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATctgc                                                  <  1:1221911/238‑1 (MQ=255)
                                                                                                                                                                                              tgctgTTCTTGCTGTTACCGCTGGTGCTGGCGGCCCGGGTCGGGATGGCGGCGCGGTTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  1:2106930/239‑1 (MQ=255)
                                                                                                                                                                                               gctgTTGTTGTGGTTACCGTGGGTGGGGGCGGGCTTGGTGGGGATGGCGGCGCGGGTGTGAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGca                                             <  2:753619/238‑1 (MQ=255)
                                                                                                                                                                                                 tgTTTTTGCTGTTAACGCTGGTGCGGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACta                                                                             <  1:1333999/204‑1 (MQ=255)
                                                                                                                                                                                                  tttATTGCGGTTAACGCTGGGGCGGGCGGGCCGGGTCCGGATGGCTGCGTGGGTGTAAATGGCAAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAAc                                         <  1:1265765/238‑1 (MQ=255)
                                                                                                                                                                                                             tAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCCCCGGCCGGTAACCCGGGGGATCTCGGCTGGGTGGCCACCCCCCCCAAGGGGATTCTGCTTGTCATGATGCTTGCGATGGGAACGCATATTCCCCAACCTGTCCTCAGGATCCTGGCGGGGGCTTCCCCCACAGTCCCCCCCGGGACGCCCGACCAGCCTGCACaaactataccctc                               >  2:2145910/1‑227 (MQ=255)
                                                                                                                                                                                                                  gcgggTGCTGGCGGGCCGGGTACGGGTGGCTGCGGGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTc                                                                                  <  1:749919/180‑1 (MQ=255)
                                                                                                                                                                                                                      gTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGGGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCCATgggggtt                                                                                                                                                         >  1:1285232/1‑103 (MQ=255)
                                                                                                                                                                                                                                   ctgttccGGAGGGCGGCGCGGGTGTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAgg          <  1:97915/233‑1 (MQ=255)
                                                                                                                                                                                                                                            aTGGCTGGGGGGTTTTTAATGGCGAACCCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAACACCTGGCATGATTTTTTGCCTTCAGGCACCGCAt  <  2:616944/238‑1 (MQ=255)
                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                               
GTACTGCTCAAGTACGGCACGCGCGATCTCAACGTCGTCTGTGGGATGCTCAAAATCATGCCATTTACCGCCGTGCTGTTTGGCGGCGGTGCGCTGGCGCTGGCAGGGATGCCGCCCTTCAACATTTTTCTTAGCGAATTTATGACCATTACCGCCGGACTGGCACGTAATCACCTGCTGATTATCGTCCTGCTGTTATTGCTGTTAACGCTGGTGCTGGCGGGCCTGGTACGGATGGCTGCGCGGGTGTTAATGGCGAAACCGCCGCAGGCCGTTAACCGGGGTGATCTCGGCTGGTTGACCACCTCGCCAATGGTGATTCTGCTGGTCATGATGCTGGCGATGGGAACGCATATTCCACAACCTGTCATCAGGATCCTGGCGGGCGCTTCCACTATAGTCCTCTCAGGGACGCACGATCTGCCTGCACAACGTAGCACCTGGCATGATTTTTTGCCTTCAGGCACCGCAT  >  NC_007779/2606653‑2607124

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