Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,075,234 C→A 91.6% 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,2340CA91.6% 486.2 / 33.6 190D296Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (5/12);  new base A (68/105);  total (73/117)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

TGCATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTT  >  CP009273/4075087‑4075377
                                                                                                                                                   |                                                                                                                                               
     cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                            <  1:2148191/148‑1 (MQ=255)
     cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                            <  1:2600916/148‑1 (MQ=255)
      aaCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtctc                                                                                                                                           >  2:892087/1‑148 (MQ=255)
       aCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCGCCGCGTATTTatct                                                                                                                                          >  1:291712/1‑148 (MQ=255)
        ccAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCtt                                                                                                                                         <  1:2303765/148‑1 (MQ=255)
         cAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTc                                                                                                                                        >  1:1493618/1‑148 (MQ=255)
            ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                     >  2:1416435/1‑148 (MQ=255)
            ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                     <  1:1581423/148‑1 (MQ=255)
            ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                     >  1:1942443/1‑148 (MQ=255)
            ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                     >  2:1416900/1‑148 (MQ=255)
               aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                  <  1:2656706/148‑1 (MQ=255)
               aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                  <  1:1003597/148‑1 (MQ=255)
               aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                  >  1:1652246/1‑148 (MQ=255)
               aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                  >  1:1812163/1‑148 (MQ=255)
                aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                 <  1:816440/148‑1 (MQ=255)
                aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                 <  1:1744772/148‑1 (MQ=255)
                 ggACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                <  2:163254/148‑1 (MQ=255)
                 ggACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                <  2:325897/148‑1 (MQ=255)
                    cGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattat                                                                                                                             <  1:1470756/148‑1 (MQ=255)
                       gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                          >  2:2295804/1‑148 (MQ=255)
                          cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCtt                                                                                                                       >  1:1810158/1‑148 (MQ=255)
                          cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCtt                                                                                                                       <  2:3038792/148‑1 (MQ=255)
                           aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTc                                                                                                                      <  1:1364826/148‑1 (MQ=255)
                            aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGTAGGCGTTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGTTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTcc                                                                                                                     <  2:1846833/148‑1 (MQ=255)
                            aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTcc                                                                                                                     <  1:2566125/148‑1 (MQ=255)
                             aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                    <  1:980607/148‑1 (MQ=255)
                             aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                    >  2:2357433/1‑148 (MQ=255)
                             aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                    <  2:315971/148‑1 (MQ=255)
                             aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                    >  2:34812/1‑148 (MQ=255)
                              aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                   >  2:2399961/1‑148 (MQ=255)
                              aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                   >  2:335494/1‑148 (MQ=255)
                              aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                   >  2:738947/1‑148 (MQ=255)
                               agcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                  <  2:2266561/148‑1 (MQ=255)
                               agcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                  >  1:3193350/1‑148 (MQ=255)
                                gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCGTTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                 <  2:76882/148‑1 (MQ=255)
                                gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                 <  1:1010051/148‑1 (MQ=255)
                                   cagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCg                                                                                                              >  2:254430/1‑148 (MQ=255)
                                     gcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCa                                                                                                            >  2:1243868/1‑148 (MQ=255)
                                     gcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCa                                                                                                            >  1:481204/1‑148 (MQ=255)
                                        cGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACgg                                                                                                         <  1:1878560/148‑1 (MQ=255)
                                          cTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                       <  1:2477929/148‑1 (MQ=255)
                                              ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGa                                                                                                   <  1:2360308/148‑1 (MQ=255)
                                                    ccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGAGTATTCTCTTCCTCATCCGCACGGGTCGGACCGaca                                                                                             >  2:2611262/1‑148 (MQ=255)
                                                        gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                         >  1:2692520/1‑148 (MQ=255)
                                                           gAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGt                                                                                      >  2:682533/1‑148 (MQ=255)
                                                           gAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGt                                                                                      >  1:1609796/1‑148 (MQ=255)
                                                             aGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAg                                                                                    <  1:1322260/148‑1 (MQ=255)
                                                                aGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTgg                                                                                 <  2:2431145/148‑1 (MQ=255)
                                                                 gAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGa                                                                                <  2:49810/148‑1 (MQ=255)
                                                                 gAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGa                                                                                <  2:1725733/148‑1 (MQ=255)
                                                                 gAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGa                                                                                <  1:1365497/148‑1 (MQ=255)
                                                                  aaGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGaa                                                                               >  2:143143/1‑148 (MQ=255)
                                                                   aGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAg                                                                              <  2:3097014/148‑1 (MQ=255)
                                                                    ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                             <  2:861290/148‑1 (MQ=255)
                                                                    ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                             >  2:1543021/1‑148 (MQ=255)
                                                                    ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                             >  1:2751441/1‑148 (MQ=255)
                                                                     gCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGAt                                                                            <  1:1239975/148‑1 (MQ=255)
                                                                     gCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGAt                                                                            <  1:2959760/148‑1 (MQ=255)
                                                                         tGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTg                                                                        >  2:817895/1‑148 (MQ=255)
                                                                           aTCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGc                                                                      >  2:2967426/1‑148 (MQ=255)
                                                                            tCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGcc                                                                     >  2:1601364/1‑148 (MQ=255)
                                                                              cAGTGTTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGc                                                                   <  2:545375/148‑1 (MQ=255)
                                                                                  ggTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                               <  1:269995/148‑1 (MQ=255)
                                                                                   gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                              <  2:1287477/148‑1 (MQ=255)
                                                                                   gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                              <  2:2721444/148‑1 (MQ=255)
                                                                                   gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                              <  2:1287456/148‑1 (MQ=255)
                                                                                       aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGc                                                          <  1:706801/148‑1 (MQ=255)
                                                                                       aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGc                                                          <  1:706720/148‑1 (MQ=255)
                                                                                                  cgcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCaa                                               <  1:2279193/148‑1 (MQ=255)
                                                                                                  cgcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCaa                                               <  1:2279246/148‑1 (MQ=255)
                                                                                                     ggTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAg                                            <  1:1423486/148‑1 (MQ=255)
                                                                                                        gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcg                                         >  1:2441072/1‑148 (MQ=255)
                                                                                                        gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGGATTTATATTAGTGCAGATTATTATATTGCTCATACGCACGGTTCGGACCGACACCCACGAAGTGGAAGATAAAGGCCGCAAAGGTATACGCAAAGCCAACAGCcg                                         >  1:1119213/1‑148 (MQ=255)
                                                                                                           gcgcACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGa                                      <  2:1810158/148‑1 (MQ=255)
                                                                                                           gcgcACGATGGTGTCACGTCGTGTCATTTTCGCTCCTCGTCGTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGa                                      >  1:928126/1‑148 (MQ=255)
                                                                                                                   tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTc                              <  1:2357433/148‑1 (MQ=255)
                                                                                                                   tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTACTCATCAGCACGGTTCGGACCCACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGACGCGAGCCGTTTc                              >  2:3224726/1‑148 (MQ=255)
                                                                                                                    ggTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                             <  1:381151/148‑1 (MQ=255)
                                                                                                                      tgtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTATCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAg                           >  2:161848/1‑148 (MQ=255)
                                                                                                                            gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                     <  1:3250322/148‑1 (MQ=255)
                                                                                                                            gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                     <  2:1076853/148‑1 (MQ=255)
                                                                                                                              cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCCACCGCCGCGAGCGGTTTCCAGTTGCCttt                   >  1:2828414/1‑148 (MQ=255)
                                                                                                                              cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                   >  1:1956123/1‑148 (MQ=255)
                                                                                                                              cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                   <  2:1285567/148‑1 (MQ=255)
                                                                                                                               gtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                  <  1:1689576/148‑1 (MQ=255)
                                                                                                                               gtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                  <  1:885213/148‑1 (MQ=255)
                                                                                                                               gtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                  <  2:1812163/148‑1 (MQ=255)
                                                                                                                               gtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                  <  1:653127/148‑1 (MQ=255)
                                                                                                                                tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                 >  2:2058341/1‑148 (MQ=255)
                                                                                                                                 ttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                >  2:720030/1‑148 (MQ=255)
                                                                                                                                     tttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAt            <  2:1030041/148‑1 (MQ=255)
                                                                                                                                      ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt           >  1:2411391/1‑148 (MQ=255)
                                                                                                                                        gcgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTc         >  1:1522884/1‑148 (MQ=255)
                                                                                                                                         cgcTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa        <  2:843930/148‑1 (MQ=255)
                                                                                                                                         cgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa        <  1:62889/148‑1 (MQ=255)
                                                                                                                                         cgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa        <  2:1652246/148‑1 (MQ=255)
                                                                                                                                             cctcGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGt    <  2:1846903/148‑1 (MQ=255)
                                                                                                                                             cctcGTATTTCTCTTCGTGCAGATTATTCTCGTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGt    <  2:2602348/148‑1 (MQ=255)
                                                                                                                                               tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt  <  2:1493618/148‑1 (MQ=255)
                                                                                                                                               tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt  <  2:290139/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TGCATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTT  >  CP009273/4075087‑4075377

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