Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,074,579 G→T 100% R4S (CGC→AGC)  fdhE ← formate dehydrogenase formation protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,074,5790GT98.6% 670.3 / 0.6 209R4S (CGC→AGC) fdhEformate dehydrogenase formation protein
Reads supporting (aligned to +/- strand):  ref base G (1/2);  new base T (109/97);  total (110/99)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 209 total aligned reads displayed.

TCACCCAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATA  >  CP009273/4074434‑4074720
                                                                                                                                                 |                                                                                                                                             
tCACCCAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTaa                                                                                                                                             >  2:1519903/1‑148 (MQ=255)
 cACCCAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAAt                                                                                                                                            >  1:2578904/1‑148 (MQ=255)
 cACCCAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAAt                                                                                                                                            >  1:2708083/1‑148 (MQ=255)
     cAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTc                                                                                                                                        <  2:1244611/148‑1 (MQ=255)
       gagGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCAt                                                                                                                                      <  2:3537203/148‑1 (MQ=255)
       gagGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCAt                                                                                                                                      <  2:1732727/148‑1 (MQ=255)
       gagGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCAt                                                                                                                                      <  2:1710328/148‑1 (MQ=255)
         ggATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTc                                                                                                                                    <  1:3003517/148‑1 (MQ=255)
           attattTTCTGCCAGCTCGCGCAGACGCTCGGCGCGCCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTCATACTCATTCAg                                                                                                                                  >  2:2033728/1‑148 (MQ=255)
            ttattTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGc                                                                                                                                 >  1:3284689/1‑148 (MQ=255)
            ttattTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGc                                                                                                                                 <  2:1798720/148‑1 (MQ=255)
             tattTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGct                                                                                                                                <  1:2264083/148‑1 (MQ=255)
                tttCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGctttc                                                                                                                             >  1:990015/1‑148 (MQ=255)
                tttCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGctttc                                                                                                                             >  1:1775232/1‑148 (MQ=255)
                  tCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGctttctt                                                                                                                           <  2:1657220/148‑1 (MQ=255)
                    tGCCAGCTCGCGCAGACGCTCGGCGCGGCGTTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCtttt                                                                                                                         <  1:300910/148‑1 (MQ=255)
                     gCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTc                                                                                                                        <  1:2726721/148‑1 (MQ=255)
                     gCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTc                                                                                                                        <  1:2791318/148‑1 (MQ=255)
                      ccAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTcc                                                                                                                       >  2:324951/1‑148 (MQ=255)
                         gCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGtt                                                                                                                    <  2:1665447/148‑1 (MQ=255)
                         gCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGtt                                                                                                                    <  2:1538753/148‑1 (MQ=255)
                          cTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTg                                                                                                                   >  1:1581588/1‑148 (MQ=255)
                          cTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTg                                                                                                                   >  1:1581313/1‑148 (MQ=255)
                            cgcgCAGACGCTCGGCGCGGCGGTTTTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTc                                                                                                                 <  1:737503/148‑1 (MQ=255)
                             gcgcAGACTCTCGGCGCGGCTGTTGTATTAATTCTTGAGCCTAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCATCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCt                                                                                                                <  2:825551/148‑1 (MQ=255)
                                  gACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCgg                                                                                                           >  2:1787193/1‑148 (MQ=255)
                                    cGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGAc                                                                                                         >  1:3469170/1‑148 (MQ=255)
                                    cGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGAc                                                                                                         <  2:192022/148‑1 (MQ=255)
                                     gCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGAcc                                                                                                        >  2:2493476/1‑148 (MQ=255)
                                           cgcggcgGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACgg                                                                                                  >  2:2945324/1‑148 (MQ=255)
                                                   ttGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAgcg                                                                                          >  1:1223944/1‑148 (MQ=255)
                                                   ttGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAgcg                                                                                          >  1:1163834/1‑148 (MQ=255)
                                                   ttGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAgcg                                                                                          >  1:1001019/1‑148 (MQ=255)
                                                    tGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAgcgc                                                                                         >  2:745560/1‑148 (MQ=255)
                                                    tGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAgcgc                                                                                         >  1:362423/1‑148 (MQ=255)
                                                     gTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGcgcg                                                                                        >  1:3356964/1‑148 (MQ=255)
                                                     gTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGcgcg                                                                                        >  2:658269/1‑148 (MQ=255)
                                                       ataAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCggg                                                                                      <  1:1940915/148‑1 (MQ=255)
                                                       ataAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCggg                                                                                      <  2:2708083/148‑1 (MQ=255)
                                                       ataAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCggg                                                                                      <  2:2578904/148‑1 (MQ=255)
                                                         aaaTTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCGTGCGGACCTCACGGTACCAGCGCGGGTg                                                                                    <  1:3067984/148‑1 (MQ=255)
                                                            ttCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATg                                                                                 >  2:1624342/1‑148 (MQ=255)
                                                             tCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGt                                                                                >  1:1711654/1‑148 (MQ=255)
                                                                 gAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTttct                                                                            <  1:2929965/148‑1 (MQ=255)
                                                                   gCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTttcttc                                                                          <  2:1239234/148‑1 (MQ=255)
                                                                      gAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGcc                                                                       <  2:2814361/148‑1 (MQ=255)
                                                                      gAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGcc                                                                       <  2:2814372/148‑1 (MQ=255)
                                                                       aGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGccc                                                                      <  1:999640/148‑1 (MQ=255)
                                                                         ggAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAg                                                                    >  1:622946/1‑148 (MQ=255)
                                                                          gAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAgg                                                                   >  2:773663/1‑148 (MQ=255)
                                                                            aCAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcg                                                                 >  2:1580855/1‑148 (MQ=255)
                                                                             cAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGgcgc                                                                <  1:1981392/148‑1 (MQ=255)
                                                                               aTAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTg                                                              <  2:1891867/148‑1 (MQ=255)
                                                                               aTAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTg                                                              >  2:1600651/1‑148 (MQ=255)
                                                                                 aaCGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGt                                                            >  2:3155381/1‑148 (MQ=255)
                                                                                 aaCGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGt                                                            <  2:228215/148‑1 (MQ=255)
                                                                                  aCGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTa                                                           >  2:3285605/1‑148 (MQ=255)
                                                                                    ggcggAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAAc                                                         <  1:1546454/148‑1 (MQ=255)
                                                                                      cggAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAAccc                                                       <  1:3508466/148‑1 (MQ=255)
                                                                                      cggAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAAccc                                                       >  2:1005514/1‑148 (MQ=255)
                                                                                      cggAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAAccc                                                       >  2:1005482/1‑148 (MQ=255)
                                                                                        gAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCAt                                                     >  2:2836600/1‑148 (MQ=255)
                                                                                        gAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCAt                                                     <  2:1168216/148‑1 (MQ=255)
                                                                                         aaTCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATc                                                    >  2:3040909/1‑148 (MQ=255)
                                                                                         aaTCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATc                                                    <  2:3020415/148‑1 (MQ=255)
                                                                                         aaTCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATc                                                    >  2:849302/1‑148 (MQ=255)
                                                                                           tCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCt                                                  >  1:929363/1‑148 (MQ=255)
                                                                                           tCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCt                                                  >  1:1099904/1‑148 (MQ=255)
                                                                                           tCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCt                                                  >  1:3537679/1‑148 (MQ=255)
                                                                                            cATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCtt                                                 >  2:217733/1‑148 (MQ=255)
                                                                                            cATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCtt                                                 <  2:3356964/148‑1 (MQ=255)
                                                                                              tatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTcc                                               >  1:2121582/1‑148 (MQ=255)
                                                                                              tatCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTcc                                               <  2:2901213/148‑1 (MQ=255)
                                                                                              tatCCGCCGTACGTTTCTCGCTTGAACACAGCTAATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGCACCTCACGGTACAAGCGCGGGTGATGTTTCTTCGCCCATGCGCTGGTAACCCATCATTcc                                               >  1:1506803/1‑148 (MQ=255)
                                                                                                    ccgTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATg                                         <  2:390789/148‑1 (MQ=255)
                                                                                                    ccgTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATg                                         >  2:1278832/1‑148 (MQ=255)
                                                                                                      gTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGc                                       >  1:877136/1‑148 (MQ=255)
                                                                                                        aCGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCgg                                     >  1:2140303/1‑148 (MQ=255)
                                                                                                          gTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTa                                   <  2:794535/148‑1 (MQ=255)
                                                                                                           tttCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTaa                                  >  2:1413495/1‑148 (MQ=255)
                                                                                                                 gCTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGc                            <  1:1624342/148‑1 (MQ=255)
                                                                                                                  cTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGcc                           <  2:1831297/148‑1 (MQ=255)
                                                                                                                  cTTGAACCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGcc                           <  1:1149156/148‑1 (MQ=255)
                                                                                                                      aaCCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTc                       >  2:1856475/1‑148 (MQ=255)
                                                                                                                      aaCCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTc                       >  2:635015/1‑148 (MQ=255)
                                                                                                                       aCCCAGCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCa                      >  2:2654971/1‑148 (MQ=255)
                                                                                                                            gCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCaa                 <  1:3115041/148‑1 (MQ=255)
                                                                                                                            gCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCaa                 <  2:3519562/148‑1 (MQ=255)
                                                                                                                            gCTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCaa                 <  2:2124723/148‑1 (MQ=255)
                                                                                                                             cTCATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCaaa                >  1:881882/1‑148 (MQ=255)
                                                                                                                               cATCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAgg              >  2:2574405/1‑148 (MQ=255)
                                                                                                                                 tCTTGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGggc            <  1:1843038/148‑1 (MQ=255)
                                                                                                                                    tGCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGggcggc         >  1:2460988/1‑148 (MQ=255)
                                                                                                                                     gCGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGgcggcg        <  1:657988/148‑1 (MQ=255)
                                                                                                                                      cGGGATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCAACAATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGt       >  1:2999716/1‑148 (MQ=255)
                                                                                                                                        ggATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAg     >  1:1533705/1‑148 (MQ=255)
                                                                                                                                        ggATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAg     >  1:10790/1‑148 (MQ=255)
                                                                                                                                         gATTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGa    >  1:3409731/1‑148 (MQ=255)
                                                                                                                                          aTTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGat   >  2:417534/1‑148 (MQ=255)
                                                                                                                                          aTTATGCTAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGat   <  1:133217/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
TCACCCAGAGGATTATTTTCTGCCAGCTCGCGCAGACGCTCGGCGCGGCGGTTGTATAAATTCTTGAGCCGAGGGAACAATAACGGCGGAATCATATCCGCCGTACGTTTCTCGCTTGAACCCAGCTCATCTTGCGGGATTATGCGAATACTCATTCAGCTTTCTTTTCCGTTGTCTTGCGGACCTCACGGTACCAGCGCGGGTGATGTTTCTTCGCCCAGGCGCTGGTAACCCATCCTTCCACCATGGCGGTAATCGTGCCTTTCACCCAAAGGGCGGCGTAGATA  >  CP009273/4074434‑4074720

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