Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,075,234 C→A 100% D296Y (GAC→TAC)  fdoH ← formate dehydrogenase‑O, Fe‑S subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,075,2340CA96.1% 844.4 / 21.6 285D296Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (2/9);  new base A (110/164);  total (112/173)
Fisher's exact test for biased strand distribution p-value = 2.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.36e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 285 total aligned reads displayed.

ATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075090‑4075378
                                                                                                                                                |                                                                                                                                                
aTCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtc                                                                                                                                              <  1:2240482/149‑1 (MQ=255)
       aGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGt                                                                                                                                       <  2:424642/149‑1 (MQ=255)
        gTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                      <  2:3370605/149‑1 (MQ=255)
        gTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                      >  2:1488602/1‑149 (MQ=255)
            aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  <  1:361376/149‑1 (MQ=255)
            aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  <  2:1117392/149‑1 (MQ=255)
            aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  >  1:1813610/1‑149 (MQ=255)
             aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                 >  2:1922416/1‑149 (MQ=255)
             aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                 >  2:3566261/1‑149 (MQ=255)
                aCGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattat                                                                                                                              <  2:4059264/149‑1 (MQ=255)
                  gggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtc                                                                                                                            <  2:1287813/149‑1 (MQ=255)
                   ggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           >  2:1599696/1‑149 (MQ=255)
                     aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTCCTCGTATTTCTCTTCGTGCAGATTATTctct                                                                                                                         >  1:1113550/1‑149 (MQ=255)
                     aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTctct                                                                                                                         <  1:2645126/149‑1 (MQ=255)
                      aCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGAATATTCTCtt                                                                                                                        >  2:1699437/1‑149 (MQ=255)
                       cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTc                                                                                                                       >  2:2765373/1‑149 (MQ=255)
                        aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTcc                                                                                                                      <  1:3194870/149‑1 (MQ=255)
                         aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     <  1:3270725/149‑1 (MQ=255)
                         aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     <  2:2663877/149‑1 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                   >  2:1212991/1‑149 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                   >  2:1163195/1‑149 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                   <  1:2181573/149‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 >  2:2087261/1‑149 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCGCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 <  2:3009732/149‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTAATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 >  1:3842146/1‑149 (MQ=255)
                              cccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATcc                                                                                                                <  1:3482631/149‑1 (MQ=255)
                              cccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCGCGTATTTCTCTTCGTGCAGTGTATTCTCTTCCTCATcc                                                                                                                >  2:1354720/1‑149 (MQ=255)
                                cagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGc                                                                                                              >  1:2241244/1‑149 (MQ=255)
                                 agcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCa                                                                                                             <  1:617976/149‑1 (MQ=255)
                                  gcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCAc                                                                                                            <  2:2144234/149‑1 (MQ=255)
                                      gCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                        <  1:886999/149‑1 (MQ=255)
                                      gCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                        >  2:1773325/1‑149 (MQ=255)
                                       cTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGAGCTCCTAGTATTTATAATAGTGGAGATTATTCTCATACTCATCCGCAcagatc                                                                                                       >  1:1635926/1‑144 (MQ=255)
                                           ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGat                                                                                                   <  1:2560132/149‑2 (MQ=255)
                                           ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAc                                                                                                   <  1:1293078/149‑1 (MQ=255)
                                           ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAc                                                                                                   <  2:3976653/149‑1 (MQ=255)
                                            cgccgcCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAcc                                                                                                  <  1:673720/149‑1 (MQ=255)
                                            cgccgcCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAcc                                                                                                  <  1:1014555/149‑1 (MQ=255)
                                                gccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGaca                                                                                              <  2:2777011/149‑1 (MQ=255)
                                                 ccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGacac                                                                                             >  1:3266374/1‑149 (MQ=255)
                                                   aGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACAccc                                                                                           <  2:2869630/149‑1 (MQ=255)
                                                    ggcggcagcagcAGGCGGTGATCCAGGGGTTGATACGTGCCGGCGCGGTGTAGCGCACGATGGGGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          <  1:2152389/143‑1 (MQ=255)
                                                    ggATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          <  1:1237802/149‑1 (MQ=255)
                                                       tGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGt                                                                                       <  2:866507/149‑1 (MQ=255)
                                                        gAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTa                                                                                      <  2:1663724/149‑1 (MQ=255)
                                                         aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCGCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAg                                                                                     >  1:2993218/1‑149 (MQ=255)
                                                           gCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAGCCGCACGGTTCGGACCGACACCCACGTAGTg                                                                                   <  2:2197716/149‑1 (MQ=255)
                                                            cAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGTCACCCACGTAGTgg                                                                                  <  1:3347068/149‑1 (MQ=255)
                                                             aGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGa                                                                                 <  2:1264132/149‑1 (MQ=255)
                                                                   cGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATAc                                                                           <  1:1733637/149‑1 (MQ=255)
                                                                    ggTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACt                                                                          >  1:3927684/1‑149 (MQ=255)
                                                                      tGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTgg                                                                        >  2:3773248/1‑149 (MQ=255)
                                                                              tGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                                <  1:816510/149‑1 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGt                                                              <  1:289351/149‑1 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGt                                                              <  2:3684603/149‑1 (MQ=255)
                                                                                 ttGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTa                                                             <  1:2745191/149‑1 (MQ=255)
                                                                                  tGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAg                                                            <  1:711517/149‑1 (MQ=255)
                                                                                    aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGcc                                                          <  1:2922065/149‑1 (MQ=255)
                                                                                       cGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCGCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCa                                                       <  2:2304122/149‑1 (MQ=255)
                                                                                           ccGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGc                                                   <  1:2377437/149‑1 (MQ=255)
                                                                                            cGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGcc                                                  <  1:187577/149‑1 (MQ=255)
                                                                                            cGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGcc                                                  <  2:3184034/149‑1 (MQ=255)
                                                                                                 cgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCCAccg                                             >  2:3860918/1‑147 (MQ=255)
                                                                                                  ggTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGCAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGc                                            >  2:921364/1‑149 (MQ=255)
                                                                                                  ggTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGc                                            <  2:3634643/149‑1 (MQ=255)
                                                                                                     gtAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTCCTCGTATTTCTCTTCGTGCCGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcgc                                         >  1:2483098/1‑149 (MQ=255)
                                                                                                     gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcgc                                         <  1:1011087/149‑1 (MQ=255)
                                                                                                            aCGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGt                                  <  1:3406570/149‑1 (MQ=255)
                                                                                                              gATGGTGTCACGTCGTTTCATTTGCGCGCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                <  2:22384/149‑1 (MQ=255)
                                                                                                                tGGTGTCTCGTCGTTTCATTTGCGCTCCTCTTATTTCTCTTGGTGCAGTTTATTCTCTTCCTCATCCGCTCGTTTCGGACCGACACGAAGGTAGTGTAAGATAGTGGCCGCTAAGGTAGCCTCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  1:1175093/149‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  2:1418749/149‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  2:536827/149‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  2:3927680/149‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              >  2:3365800/1‑149 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCAACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  1:1531691/149‑1 (MQ=255)
                                                                                                                         gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCtt                     >  2:2983482/1‑149 (MQ=255)
                                                                                                                          tcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    >  2:2018985/1‑149 (MQ=255)
                                                                                                                          tcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    >  2:210624/1‑149 (MQ=255)
                                                                                                                          tcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    <  1:735052/149‑1 (MQ=255)
                                                                                                                          tcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTATCTTACTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    >  1:3649403/1‑149 (MQ=255)
                                                                                                                            gtTTCATTTGCGCTCCTAGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  >  1:2781557/1‑149 (MQ=255)
                                                                                                                             tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  1:492398/1‑149 (MQ=255)
                                                                                                                             tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  1:609965/1‑149 (MQ=255)
                                                                                                                             tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  1:2007871/1‑149 (MQ=255)
                                                                                                                              ttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAg                >  2:746460/1‑149 (MQ=255)
                                                                                                                                cATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGaa              <  2:2823329/149‑1 (MQ=255)
                                                                                                                                cATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGaa              <  1:2631074/149‑1 (MQ=255)
                                                                                                                                 aTTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAt             >  2:3194687/1‑149 (MQ=255)
                                                                                                                                  tttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            <  2:1813610/149‑1 (MQ=255)
                                                                                                                                  tttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            <  1:453772/149‑1 (MQ=255)
                                                                                                                                    tGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGNGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTc          <  1:2087261/149‑1 (MQ=255)
                                                                                                                                     gcgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa         <  2:3842146/149‑1 (MQ=255)
                                                                                                                                     gcgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa         >  2:2999177/1‑149 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACTGTTCGGACCGACTCCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATTCCTTTCCAGAATTTCACg       <  2:461256/149‑1 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       <  2:2529321/149‑1 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       <  1:1694164/149‑1 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       <  1:139260/149‑1 (MQ=255)
                                                                                                                                           ctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt   <  1:947643/149‑1 (MQ=255)
                                                                                                                                           ctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt   <  1:3389394/149‑1 (MQ=255)
                                                                                                                                            tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTc  <  1:737309/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075090‑4075378

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