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% 649.9 / 5.8 216D296Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (0/5);  new base A (88/123);  total (88/128)
Fisher's exact test for biased strand distribution p-value = 8.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.87e-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 216 total aligned reads displayed.

ATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075090‑4075378
                                                                                                                                                |                                                                                                                                                
aTCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCttt                                                                                                                                               <  2:897616/148‑1 (MQ=255)
aTCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTAttt                                                                                                                                               <  1:3735435/148‑1 (MQ=255)
 tCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTAGTATTtc                                                                                                                                              >  2:691802/1‑148 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTtct                                                                                                                                             <  1:2706864/148‑1 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             <  2:1109664/148‑1 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             >  1:2738870/1‑148 (MQ=255)
  cAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtct                                                                                                                                             <  1:1327340/148‑1 (MQ=255)
    aCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTctct                                                                                                                                           <  2:3404591/148‑1 (MQ=255)
     ccAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCtt                                                                                                                                          <  2:1685013/148‑1 (MQ=255)
     ccAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCtt                                                                                                                                          <  2:3581807/148‑1 (MQ=255)
        gTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGt                                                                                                                                       <  1:1395164/148‑1 (MQ=255)
         ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                      >  2:229472/1‑148 (MQ=255)
         ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                      >  2:132443/1‑148 (MQ=255)
         ttGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTg                                                                                                                                      <  1:1590588/148‑1 (MQ=255)
             aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  <  2:3786134/148‑1 (MQ=255)
              ggACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                 >  2:1766125/1‑148 (MQ=255)
                aCGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatta                                                                                                                               <  1:3164645/148‑1 (MQ=255)
                 cgtgAACAAAAAGCCCAGCCCGCTCACCGCCGCAAGGAAGAAGCAGAAGGCGGTTATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCTCTCCTATTATTTCTCTTCGTGCAGattat                                                                                                                              <  2:2905703/145‑1 (MQ=255)
                 cGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattat                                                                                                                              <  2:2906799/148‑1 (MQ=255)
                  gggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattatt                                                                                                                             <  1:1283802/148‑1 (MQ=255)
                  gggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattatt                                                                                                                             <  2:1929394/148‑1 (MQ=255)
                   ggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtc                                                                                                                            <  2:1913735/148‑1 (MQ=255)
                    gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           <  1:3406384/148‑1 (MQ=255)
                    gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           <  1:3405417/148‑1 (MQ=255)
                     aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtctc                                                                                                                          <  2:432026/148‑1 (MQ=255)
                     aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtctc                                                                                                                          >  1:1184283/1‑148 (MQ=255)
                     aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtctc                                                                                                                          <  1:3473299/148‑1 (MQ=255)
                      aCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTatct                                                                                                                         >  1:832408/1‑148 (MQ=255)
                        aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTCCTCGTATTTATCTTCGTGCAGATTATTGTATTc                                                                                                                       >  1:1713651/1‑148 (MQ=255)
                        aaaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTc                                                                                                                       <  1:527372/148‑1 (MQ=255)
                          aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     >  1:3378255/1‑148 (MQ=255)
                          aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     <  2:1772542/148‑1 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    <  1:3584648/148‑1 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    >  1:1468107/1‑148 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    <  1:1063553/148‑1 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    <  1:291721/148‑1 (MQ=255)
                           aaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTc                                                                                                                    <  2:1861647/148‑1 (MQ=255)
                            agcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                   >  1:1273057/1‑148 (MQ=255)
                              cccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 <  2:1417651/148‑1 (MQ=255)
                                 agcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGc                                                                                                              <  1:408325/148‑1 (MQ=255)
                                   cccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCAc                                                                                                            >  2:2296939/1‑148 (MQ=255)
                                      gCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTCGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGt                                                                                                         <  1:207047/148‑1 (MQ=255)
                                       cTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                        <  2:3038771/148‑1 (MQ=255)
                                           ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGa                                                                                                    <  1:707978/148‑1 (MQ=255)
                                               cgccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGa                                                                                                <  2:1237908/148‑1 (MQ=255)
                                                    ggATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACAccc                                                                                           >  1:3192579/1‑148 (MQ=255)
                                                     gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          >  1:2995404/1‑148 (MQ=255)
                                                     gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTATCTTCGTGCAGATTATTATCTTCCTCATACGCACGGTTCGGACCGACACCCa                                                                                          >  1:3671950/1‑148 (MQ=255)
                                                      aTGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGGTCCTCGTATTTCTCTTCGTGCAGATTATTGTCTTCCTCATCCGCACGGTTCGGACCGGCACCCAc                                                                                         >  1:2447578/1‑148 (MQ=255)
                                                       tGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACg                                                                                        >  1:2290105/1‑148 (MQ=255)
                                                         aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTa                                                                                      <  1:2328871/148‑1 (MQ=255)
                                                         aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGGATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTa                                                                                      >  1:1148166/1‑148 (MQ=255)
                                                                aGGCGGTGATCCAGTGGTTGATACGTTGCGGCGCGGTGTAGCGCTCGATGGTGTGACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAg                                                                               <  2:3394991/148‑1 (MQ=255)
                                                                 ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                              >  2:2082244/1‑148 (MQ=255)
                                                                 ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                              >  2:3763062/1‑148 (MQ=255)
                                                                    ggTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATAc                                                                           <  2:1009533/148‑1 (MQ=255)
                                                                      tGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTg                                                                         >  2:2429527/1‑148 (MQ=255)
                                                                            aGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCa                                                                   >  2:1542623/1‑148 (MQ=255)
                                                                             gTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCaa                                                                  <  1:862684/148‑1 (MQ=255)
                                                                              tGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCaaa                                                                 <  2:2738870/148‑1 (MQ=255)
                                                                              tGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCaaa                                                                 <  2:824859/148‑1 (MQ=255)
                                                                               ggTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                                >  2:578579/1‑148 (MQ=255)
                                                                                 ttGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTATATTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGt                                                              >  1:424648/1‑148 (MQ=255)
                                                                                 ttGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATCCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGt                                                              <  1:2765721/148‑1 (MQ=255)
                                                                                  tGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTa                                                             >  2:203106/1‑148 (MQ=255)
                                                                                   gATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAg                                                            <  1:1104417/148‑1 (MQ=255)
                                                                                    aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGc                                                           <  2:1783291/148‑1 (MQ=255)
                                                                                             ggCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAg                                                            >  2:2198655/1‑138 (MQ=255)
                                                                                                gcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAAc                                               <  2:1539059/148‑1 (MQ=255)
                                                                                                gcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAAc                                               <  2:561887/148‑1 (MQ=255)
                                                                                                  ggTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAg                                             <  2:3460239/148‑1 (MQ=255)
                                                                                                   gtgtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGc                                            >  2:1599584/1‑148 (MQ=255)
                                                                                                   gtgtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGc                                            <  1:3622027/148‑1 (MQ=255)
                                                                                                    tgtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGcc                                           >  2:3779884/1‑148 (MQ=255)
                                                                                                     gtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcg                                          <  1:1253295/148‑1 (MQ=255)
                                                                                                      tAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCcgc                                         >  2:3390767/1‑148 (MQ=255)
                                                                                                         cgcACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAg                                      <  2:3320007/148‑1 (MQ=255)
                                                                                                             cGATGGTGTCACGTCGTTTAATTTTCGCTGCTCGGATTTATAGTCGTGGAGATTATTCTCTTACTAATCCGCACCGTTAGGACCGACACCCACCCAGTGGAAGATAAAGGCCGCAAAGGTACCCCCAAAACCACAACCCCCGAGCGGt                                  >  1:3232401/1‑148 (MQ=255)
                                                                                                               aTGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                <  1:2620373/148‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTc                               <  1:2793516/148‑1 (MQ=255)
                                                                                                                tGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTc                               <  1:1812336/148‑1 (MQ=255)
                                                                                                                   tgtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAg                            <  2:3063142/148‑1 (MQ=255)
                                                                                                                    gtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGa                           <  2:1273057/148‑1 (MQ=255)
                                                                                                                    gtCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGa                           <  1:3282158/148‑1 (MQ=255)
                                                                                                                        cgtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGcc                       <  2:2950328/148‑1 (MQ=255)
                                                                                                                         gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                      <  2:2492894/148‑1 (MQ=255)
                                                                                                                             tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  <  2:76088/148‑1 (MQ=255)
                                                                                                                             tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  <  2:76069/148‑1 (MQ=255)
                                                                                                                              ttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  2:108710/1‑148 (MQ=255)
                                                                                                                              ttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCa                 >  1:1962068/1‑148 (MQ=255)
                                                                                                                                cATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGa               <  2:1715385/148‑1 (MQ=255)
                                                                                                                                  tttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAt             <  1:3091050/148‑1 (MQ=255)
                                                                                                                                  tttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAt             <  1:1120037/148‑1 (MQ=255)
                                                                                                                                       gcTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCAc        <  2:1467763/148‑1 (MQ=255)
                                                                                                                                        ctcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       >  2:69987/1‑148 (MQ=255)
                                                                                                                                        ctcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       >  2:2000147/1‑148 (MQ=255)
                                                                                                                                        ctcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       >  2:2000910/1‑148 (MQ=255)
                                                                                                                                        ctcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       >  2:1002710/1‑148 (MQ=255)
                                                                                                                                         tcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACgg      <  2:3114593/148‑1 (MQ=255)
                                                                                                                                             cGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTc  >  1:3514771/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                                
ATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075090‑4075378

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: ‑