Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,598,891 G→T 8.4% S241I (AGC→ATC)  hyfD → hydrogenase 4, membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,598,8910GT 8.4% 158.2 / 4.2 56S241I (AGC→ATC) hyfDhydrogenase 4, membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (28/23);  new base T (2/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.53e-01

TTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTCGCGGC  >  CP009273/2598623‑2599139
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
ttGATCAGTTATTACCAGAGCGATAATGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTATTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCNTCCTGTTTNCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCt                                                                                                                                                                                                                                                     >  1:491796/1‑274 (MQ=255)
   ttCAGTTATTACCAGAGCGTTAAAGCGAAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTCGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTAtct                                                                                                                                                                                                                                                 <  1:693605/274‑1 (MQ=255)
                  agCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCACGTGGGGGAAATCGGCCCAGCGACCGATGCAAGCGTGCCTACCGGACGCAATGGACGCGCCAACACCGATCAGCGCATATCCCCACGCCGCATCGac                                                                                                                                                                                                                                  >  1:399250/1‑274 (MQ=255)
                                     cAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGAACAGCGCCTATCTCCACGCCGACACGAAGGTGAAAATGGGGGTTTa                                                                                                                                                                                                                >  1:34028/1‑274 (MQ=255)
                                                  gTTTCTTATTACTCATTTCGGCTTGTTGGGGTTGTAGTTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGATTAGAGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATACTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCcgcg                                                                                                                                                                                                  <  1:692292/275‑1 (MQ=255)
                                                   ttACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTCCCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCgcgc                                                                                                                                                                                                 >  1:105104/1‑275 (MQ=255)
                                                           cACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAAGGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGCTGGTGAACGGGGGCGTTTACATTTTTGCCCGCGCTATTATc                                                                                                                                                                                          >  2:735371/1‑274 (MQ=255)
                                                           cACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTCTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGCCGCAATGGAAGCTCCAACACCGATCAGCGCCTAGCTCCCAGCCACAACAATGGTGAAAGAGGGCGTTTACATTTTTGACCGCGCTATTATc                                                                                                                                                                                          >  1:703913/1‑274 (MQ=255)
                                                                          gTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGATGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGATTACATTTTTGCCCTCGCTATTATCGAcctgtctcttata                                                                                                                                                                           >  1:404079/1‑262 (MQ=255)
                                                                           ttGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTCTCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGGAATGGAAGCGCACACACCGATCAg                                                                                                                                                                                                                                                          >  2:519694/1‑194 (MQ=255)
                                                                               ggTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGATCGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCAt                                                                                                                                                                     <  1:213207/275‑1 (MQ=255)
                                                                               ggTTGTATCTNGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCAGCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACACAATGGACGCGCCTACACCGATCAGCGcc                                                                                                                                                                                                                                                      <  1:519703/194‑1 (MQ=255)
                                                                                gTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAGGCGTGGCTACNGGACGCAATGGAAGCTCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCAtg                                                                                                                                                                    >  1:666152/1‑275 (MQ=255)
                                                                                           ccgccgcgACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCCTGGCTACCGGCCGCAATGGAAGCGCCCACACCGATCAGCGCCTATCCCCCCGCCGCATCCATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCCCCAGTATCGACGGCGGCACTACCCCCCAGGTGGTTggg                                                                                                                                                            >  2:739456/1‑271 (MQ=255)
                                                                                              ccgccACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGGTTGCCGCGTGGGGGGAATCGGCCCAGCTACCGCTGCAAGCGTGGCTACCGGACGGAATGGAAGCGCCAACCCCGATCAGCCCCTCTCCCCCCCCCGCATCCCTGGGGAAAGCGGGGTTTTACATTTTGGGCCGCGCTATTAGTCACGGCGGCAATTGCCCGCATGTTATTggggtggt                                                                                                                                                       >  2:252730/1‑268 (MQ=255)
                                                                                                   aggggttttttGAGGACCGGAGCGGGTGGGTTTAGCGCGATGAGGGAGTTAGAGGGCGACTCATGTGATCTGGTTTAGGGCGGGATCCGGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGAGGAAAGCGTGGATAGCGGACGCAATGGAAGCGCCAACACCGATCAGCGCATATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCAt                                                                                                                                                 <  2:666138/269‑1 (MQ=255)
                                                                                                                 gACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGGGTGGCTACCGGACGCAATGGAAGCGCCAACACCGANCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGc                                                                                                                                                                                                                     >  1:622493/1‑193 (MQ=255)
                                                                                                                   ccGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGCCGCACGTTCTCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGCAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCACCACCGCTCAGCGCCTATCTCCACGCCGCATCGACGGTGCAAGTGGGCGTTTACATTTTTGCCCGCGCTCTTATCGACGGCGGCCATCTCACGCAGGTGCTAGGCGGCGGTGGCGTGGTCATGGCACTGGTc                                                                                                                                 >  2:677155/1‑275 (MQ=255)
                                                                                                                       aaCGTTTGCGCTTAGCGCNATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGcc                                                                                                                                                                                                      <  1:573762/202‑1 (MQ=255)
                                                                                                                        acttttGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGCGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCa                                                                                                                             <  1:735385/271‑1 (MQ=255)
                                                                                                                        aCGTTTGCGCTTGGGGGGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGAGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCCCTGGTCACCa                                                                                                                             <  1:739470/274‑1 (MQ=255)
                                                                                                                              gcgcTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTACCTGGTTTATGTCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCAACGGCGGCAATACACCGCATGTGATTGGCGGGCGTGGTAGGGGGCGGGCACTGGTCTCCATTCtttt                                                                                                                      >  2:7082/1‑274 (MQ=255)
                                                                                                                                                     acacGGCGACGCACGTTAGCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCca                                                                                                <  1:360900/274‑1 (MQ=255)
                                                                                                                                                        cGGCGACGCACGTTATCTGGTTTTTGGCGGCAGTCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCAcagc                                                                                            <  1:133606/275‑1 (MQ=255)
                                                                                                                                                         ggCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCAcagc                                                                                            <  1:252733/274‑1 (MQ=255)
                                                                                                                                                           cGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGcc                                                                                                                                                                                                      >  1:303148/1‑166 (MQ=255)
                                                                                                                                                             aCGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCACCACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCCTATGATTGGCGGCGTTGGCATGGTCATGGCCCTGGTCACCATGCTTTATGGCTTTTTGATGTTTTGccctccgcatgat                                                                                      >  2:133605/1‑264 (MQ=255)
                                                                                                                                                               gCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCcgcg                                                                                                                                                                                                  <  2:303143/166‑1 (MQ=255)
                                                                                                                                                                     tATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGATACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGAGCAGCGCCTATCGCCACGCCGCATCGATGGTGAAAGTGGGAGTTTAAATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCGCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGGAGGATATGAAGCg                                                                               >  1:87333/1‑275 (MQ=255)
                                                                                                                                                                                 ggcggcATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGCGTGGAGGCGTTg                                                                                                                                                     >  2:666240/1‑193 (MQ=255)
                                                                                                                                                                                 ggcggcATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCAAAGCAGGATATGAAGCGGTTGCTGGCCTg                                                                   >  1:22925/1‑275 (MQ=255)
                                                                                                                                                                                   cgtgagcatatttGCCCCGTGGGTCAAAACGGCCCACCGACCGATGCAAGCTTGGCTTCCGGACGCAATGTCAGCGCCAACACCGCGCAGCGCCTCTCTCCACGCTGCATCGCTGGTGAAATTGGGCGTTTACATTTTTGCCCGCGCTATTATCTACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATAGTCATTGCACTTGTCCCCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGt                                                                 <  2:87332/265‑1 (MQ=255)
                                                                                                                                                                                     gcATCCTGTTTGCCGCGTGGGGGAAATCGGCCCANCTACCGATGCAGGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCAt                                                                                                                                                 <  1:666254/193‑1 (MQ=255)
                                                                                                                                                                                      cATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGcc                                                                                                                                                                                                      >  2:48391/1‑139 (MQ=255)
                                                                                                                                                                                       ttgcTGTTCGCCTCGTCGGGGAAACCGCCCCACCTACCGATGCAAGCCTGGCTACAGGACGCGATGGCGGCGCCAACAACGCTCCTGTCCTATCGCCACGCAGCATCGATGGTGAAGGTGGGCGTTTGTTTTTTTGGCCGGGCTATGATCGACGGCGGTAATATCCCGCATGTGATTGGCCGCGTTTGCATTGTCATGGCGCTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGAc                                                             <  2:22924/272‑1 (MQ=255)
                                                                                                                                                                                         ccTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAg                                                                                                                                                                                                                                                          <  1:142399/84‑1 (MQ=255)
                                                                                                                                                                                          cTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCcgcg                                                                                                                                                                                                  <  1:48392/139‑1 (MQ=255)
                                                                                                                                                                                          cTGTGTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCATCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGNGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGAc                                                             <  1:611163/272‑1 (MQ=255)
                                                                                                                                                                                             tttGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACAttttt                                                                                                                                                                                                         >  2:585392/1‑129 (MQ=255)
                                                                                                                                                                                                 ccGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGccc                                                                                                                                                                                                     <  1:585405/129‑1 (MQ=255)
                                                                                                                                                                                                  cgcgTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCATCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCAATGGTCACCATTCTTTATGGCTTTCTGATGTAGTTTGCACCGCAGGATATGAAGCGGTTGAGGGCCTGGTTGAccttgacgcta                                                   >  2:611150/1‑265 (MQ=255)
                                                                                                                                                                                                  cgcgTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGcc                                                                                                                                                                                                      >  2:519489/1‑127 (MQ=255)
                                                                                                                                                                                                      tGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGAcgg                                                                                                                                                                                     >  2:436614/1‑140 (MQ=255)
                                                                                                                                                                                                      tGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGAGCAAAGTGGGCGTTTACATTTTTGCCcgcg                                                                                                                                                                                                  <  1:519498/127‑1 (MQ=255)
                                                                                                                                                                                                          ggAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACggcggc                                                                                                                                                                                 <  1:436621/140‑1 (MQ=255)
                                                                                                                                                                                                              aTCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGACAGTGGGCGTTTACATTTTTGCCCGCGCTATTCTCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGCCTTTCTGAGGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCAGGTAGACCATCACTCAACTTGACTGGaag                                       >  1:228511/1‑272 (MQ=255)
                                                                                                                                                                                                                                gCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCg                                                                                                                                                                                         >  1:734334/1‑110 (MQ=255)
                                                                                                                                                                                                                                gCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGccc                                                                                                >  2:267194/1‑198 (MQ=255)
                                                                                                                                                                                                                                    gCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGAcgg                                                                                                                                                                                     <  2:734320/110‑1 (MQ=255)
                                                                                                                                                                                                                                    gCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCAcagc                                                                                            <  1:267198/199‑1 (MQ=255)
                                                                                                                                                                                                                                         gCTACCGGACGCAATGGAAGCGCCAACACCGATCATCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGa                                                                                                           >  2:268478/1‑179 (MQ=255)
                                                                                                                                                                                                                                             ccGGACGCAATGGAAGCGCCAACACCGATCATCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTa                                                                                                       <  1:268482/179‑1 (MQ=255)
                                                                                                                                                                                                                                                  cGCAATGGAAGCGCCAACACCGTTCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTCGCGGc  <  1:7082/275‑1 (MQ=255)
                                                                                                                                                                                                                                                               ccAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCg                                                                                                                                                                                         >  1:486605/1‑79 (MQ=255)
                                                                                                                                                                                                                                                                   cacCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGcc                                                                                                                                                                                                      >  2:634900/1‑62 (MQ=255)
                                                                                                                                                                                                                                                                   cacCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGAcgg                                                                                                                                                                                     <  2:486596/79‑1 (MQ=255)
                                                                                                                                                                                                                                                                       gATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCcgcg                                                                                                                                                                                                  <  1:634914/62‑1 (MQ=255)
                                                                                                                                                                                                                                                                        aTCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACtt                                                >  2:656730/1‑207 (MQ=255)
                                                                                                                                                                                                                                                                            gcgcCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCTCGCTACTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCt                                            <  1:656744/207‑1 (MQ=255)
                                                                                                                                                                                                                                                                            |                                                                                                                                                                                                                                                        
TTGATCAGTTATTACCAGAGCGATAAAGCGCAGCGTTCAGCACTAAAAGCGTTACTTATCACTCATATCGGCTCGTTGGGGTTGTATCTTGCCGCCGCCACGCTGTTTTTGCAGACCGGAACGTTTGCGCTTAGCGCGATGAGCGAGTTACACGGCGACGCACGTTATCTGGTTTATGGCGGCATCCTGTTTGCCGCGTGGGGGAAATCGGCCCAGCTACCGATGCAAGCGTGGCTACCGGACGCAATGGAAGCGCCAACACCGATCAGCGCCTATCTCCACGCCGCATCGATGGTGAAAGTGGGCGTTTACATTTTTGCCCGCGCTATTATCGACGGCGGCAATATCCCGCATGTGATTGGCGGCGTTGGCATGGTCATGGCACTGGTCACCATTCTTTATGGCTTTCTGATGTATTTGCCACAGCAGGATATGAAGCGGTTGCTAGCCTGGTCGACCATCACTCAACTTGGCTGGATGTTCTTCGGCTTGTCGCTCTCCATCTTCGGCTCGCGGC  >  CP009273/2598623‑2599139

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