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,2340CA97.7% 509.2 / 3.6 171D296Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base A (70/97);  total (71/100)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 171 total aligned reads displayed.

ATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075090‑4075378
                                                                                                                                                |                                                                                                                                                
aaaaaCCCGTTGAAGGACGGGAACAAAAAGACCAGCACGCTCACCGCAGCCAGGATGAAGCAGAAGGCTGTGATCCATTTTTTGATACGTTCCGGCGCGGTGTAGCGCACTATGGTGTCACTTAGTTTCATTTTCGCTACTCGTAttt                                                                                                                                               <  2:2089461/145‑1 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             <  1:344828/148‑1 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             >  2:180142/1‑148 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             <  2:3102143/148‑1 (MQ=255)
   aaCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtctc                                                                                                                                            <  2:317650/148‑1 (MQ=255)
    aCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTctct                                                                                                                                           <  1:523570/148‑1 (MQ=255)
        gTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGt                                                                                                                                       >  1:148248/1‑148 (MQ=255)
        gTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGt                                                                                                                                       >  2:1437241/1‑148 (MQ=255)
          tGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGc                                                                                                                                     >  2:2678939/1‑148 (MQ=255)
          tGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCGCGTATTTCTCTTCGTGc                                                                                                                                     >  2:2678930/1‑148 (MQ=255)
            aaGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                   <  1:2438396/148‑1 (MQ=255)
            aaGGACGGGAAAAAAAAGCCCAGCCCTCTCACCGCCGCCAGGATGAATCAGATGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAg                                                                                                                                   <  1:1822306/148‑1 (MQ=255)
              ggACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                 <  2:1067354/148‑1 (MQ=255)
               gACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatt                                                                                                                                <  1:1743296/148‑1 (MQ=255)
                aCGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatta                                                                                                                               <  1:306712/148‑1 (MQ=255)
                    gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           <  1:3399140/148‑1 (MQ=255)
                    gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           <  1:3399336/148‑1 (MQ=255)
                      aCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTctct                                                                                                                         <  1:376765/148‑1 (MQ=255)
                       cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCtt                                                                                                                        <  2:2248257/148‑1 (MQ=255)
                       cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCtt                                                                                                                        <  2:2248629/148‑1 (MQ=255)
                        aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTc                                                                                                                       >  2:247527/1‑148 (MQ=255)
                        aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTc                                                                                                                       <  2:694466/148‑1 (MQ=255)
                         aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCAATTGCGCACCGCGTATTTCTCTTCGTGCAGTTTATTCTCTTcc                                                                                                                      >  1:2322190/1‑148 (MQ=255)
                          aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     >  2:3438061/1‑148 (MQ=255)
                          aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     >  1:1454766/1‑148 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    <  1:3250906/148‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                  <  1:1745890/148‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                  <  1:1746637/148‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                  <  1:2958/148‑1 (MQ=255)
                             gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                  <  1:2966/148‑1 (MQ=255)
                              cccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGGATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 >  2:2420256/1‑148 (MQ=255)
                                  gcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCa                                                                                                             <  2:2273833/148‑1 (MQ=255)
                                  gcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCa                                                                                                             <  1:2739065/148‑1 (MQ=255)
                                   cccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCAc                                                                                                            >  2:2436548/1‑148 (MQ=255)
                                   cccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCAc                                                                                                            >  2:2436552/1‑148 (MQ=255)
                                     cGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACgg                                                                                                          >  1:2505620/1‑148 (MQ=255)
                                      gCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGt                                                                                                         <  2:722424/148‑1 (MQ=255)
                                      gCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGt                                                                                                         <  1:2738104/148‑1 (MQ=255)
                                       ctctCCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                        <  1:1415521/148‑1 (MQ=255)
                                          aCCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCgg                                                                                                     <  1:1092171/148‑1 (MQ=255)
                                               cgccAGGAGGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGa                                                                                                <  2:2102019/148‑1 (MQ=255)
                                                     taTGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          >  1:788423/2‑148 (MQ=255)
                                                     gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          >  2:2212681/1‑148 (MQ=255)
                                                     gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          >  2:2025708/1‑148 (MQ=255)
                                                          aGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAg                                                                                     <  2:1094933/148‑1 (MQ=255)
                                                          aGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTATCTTACTCATCCGCACGGTTCGGACCGACACCCACGTAg                                                                                     >  1:3330130/1‑148 (MQ=255)
                                                              gAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGa                                                                                 <  1:1999199/148‑1 (MQ=255)
                                                                     gTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTCCTCGTATTTCTCGTAGTGCAGATTATTCTCTTCCTCATCCGCACGGGTCGGACCGACACCCACCGAGTGGAAGATAca                                                                          >  1:1621961/1‑147 (MQ=255)
                                                                     gTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACt                                                                          >  1:1621938/1‑148 (MQ=255)
                                                                     gTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACt                                                                          >  2:1672322/1‑148 (MQ=255)
                                                                       gATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTgg                                                                        <  2:2864030/148‑1 (MQ=255)
                                                                       gATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTgg                                                                        <  2:125937/148‑1 (MQ=255)
                                                                       gATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACCTAGTGGAAGATACTgg                                                                        >  1:921658/1‑148 (MQ=255)
                                                                        aTCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTCCTCGTATTTATATTCGTGCAGATTATTATCTTACTCATCCGCACGGTTCGGACCGACACCCACGGAGTGGAAGATaaaggc                                                                       >  1:2664396/1‑143 (MQ=255)
                                                                        aTCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGc                                                                       <  2:447151/148‑1 (MQ=255)
                                                                        aTCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGc                                                                       <  2:3059728/148‑1 (MQ=255)
                                                                        aTCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGc                                                                       <  2:2062981/148‑1 (MQ=255)
                                                                               ggTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                                <  2:2595287/148‑1 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               <  1:3469455/148‑1 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               >  1:1882872/1‑148 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               <  1:3468821/148‑1 (MQ=255)
                                                                                gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               <  1:75955/148‑1 (MQ=255)
                                                                                 ttGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGt                                                              <  2:2029074/148‑1 (MQ=255)
                                                                                    aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGc                                                           >  2:3116750/1‑148 (MQ=255)
                                                                                               cgcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCaa                                                >  1:1083767/1‑148 (MQ=255)
                                                                                                  ggTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAg                                             <  2:3448593/148‑1 (MQ=255)
                                                                                                     gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcg                                          <  2:907072/148‑1 (MQ=255)
                                                                                                     gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcg                                          <  1:284356/148‑1 (MQ=255)
                                                                                                         cgcACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAg                                      <  2:2821353/148‑1 (MQ=255)
                                                                                                             cGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGt                                  <  2:2964616/148‑1 (MQ=255)
                                                                                                             cGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGt                                  <  2:2745288/148‑1 (MQ=255)
                                                                                                             cGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGt                                  <  2:1169434/148‑1 (MQ=255)
                                                                                                               aTGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                >  2:3126288/1‑148 (MQ=255)
                                                                                                               aTGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                >  2:3125991/1‑148 (MQ=255)
                                                                                                               aTGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                >  2:2722271/1‑148 (MQ=255)
                                                                                                               aTGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCGTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                >  2:2722328/1‑148 (MQ=255)
                                                                                                                 ggTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              >  2:259992/1‑148 (MQ=255)
                                                                                                                  gtgtCACGTCGTTTCATTTGCGCTCCGCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCa                             >  2:1018819/1‑148 (MQ=255)
                                                                                                                   tgtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAg                            <  2:1787034/148‑1 (MQ=255)
                                                                                                                   tgtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAg                            <  1:2492031/148‑1 (MQ=255)
                                                                                                                    gtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGa                           <  1:659237/148‑1 (MQ=255)
                                                                                                                      cACGTCGTTTCATTTGCGCTACACGTATTTCACTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAATGAAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCACATg                         >  2:875526/1‑148 (MQ=255)
                                                                                                                        cgtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGcc                       <  1:2887610/148‑1 (MQ=255)
                                                                                                                           cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    >  2:1482473/1‑148 (MQ=255)
                                                                                                                           cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGGAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    >  2:3467151/1‑148 (MQ=255)
                                                                                                                            gtTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                   >  1:268982/1‑148 (MQ=255)
                                                                                                                              ttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  1:847995/1‑148 (MQ=255)
                                                                                                                               tCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAg                <  2:2977795/148‑1 (MQ=255)
                                                                                                                                 aTTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGaa              >  2:1513860/1‑148 (MQ=255)
                                                                                                                                   ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            <  1:792376/148‑1 (MQ=255)
                                                                                                                                   ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            <  2:209858/148‑1 (MQ=255)
                                                                                                                                    tGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAttt           <  2:2270234/148‑1 (MQ=255)
                                                                                                                                    tGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAttt           <  1:3245027/148‑1 (MQ=255)
                                                                                                                                      cgcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCa         >  2:1879115/1‑148 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCAc        <  1:18098/148‑1 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCAc        <  1:1797297/148‑1 (MQ=255)
                                                                                                                                         tcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACgg      >  2:2350476/1‑148 (MQ=255)
                                                                                                                                         tcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACgg      <  1:132871/148‑1 (MQ=255)
                                                                                                                                          cctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGt     <  1:23889/148‑1 (MQ=255)
                                                                                                                                           ctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGtt    <  1:2682260/148‑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: