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,2340CA98.3% 576.8 / ‑0.4 192D296Y (GAC→TAC) fdoHformate dehydrogenase‑O, Fe‑S subunit
Reads supporting (aligned to +/- strand):  ref base C (4/0);  new base A (90/98);  total (94/98)
Fisher's exact test for biased strand distribution p-value = 5.56e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
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 192 total aligned reads displayed.

CATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075089‑4075378
                                                                                                                                                 |                                                                                                                                                
cATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTAtt                                                                                                                                                >  2:2519893/1‑148 (MQ=255)
cATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTAtt                                                                                                                                                >  2:3541696/1‑148 (MQ=255)
 aTCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTAttt                                                                                                                                               <  2:1230742/148‑1 (MQ=255)
  tCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtc                                                                                                                                              >  2:2212182/1‑148 (MQ=255)
    aaCCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTtctc                                                                                                                                            <  1:1279954/148‑1 (MQ=255)
      ccAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCtt                                                                                                                                          >  2:623270/1‑148 (MQ=255)
        aGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCg                                                                                                                                        <  2:1087583/148‑1 (MQ=255)
        aGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCg                                                                                                                                        <  2:3516338/148‑1 (MQ=255)
              aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  <  1:656383/148‑1 (MQ=255)
              aGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGa                                                                                                                                  <  2:2603456/148‑1 (MQ=255)
               ggACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGat                                                                                                                                 >  2:2791455/1‑148 (MQ=255)
                 aCGGTAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatta                                                                                                                               <  1:3099579/148‑1 (MQ=255)
                 aCGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatta                                                                                                                               <  2:28011/148‑1 (MQ=255)
                 aCGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGatta                                                                                                                               >  2:2885933/1‑148 (MQ=255)
                   gggAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGattatt                                                                                                                             <  2:9192/148‑1 (MQ=255)
                     gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           <  2:2947594/148‑1 (MQ=255)
                     gAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtct                                                                                                                           >  2:2984611/1‑148 (MQ=255)
                      aaCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATtctc                                                                                                                          >  2:720263/1‑148 (MQ=255)
                       aCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTctct                                                                                                                         >  1:1992723/1‑148 (MQ=255)
                       aCAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTctct                                                                                                                         >  2:1637267/1‑148 (MQ=255)
                        cAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCtt                                                                                                                        <  1:1925417/148‑1 (MQ=255)
                          aaaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTcc                                                                                                                      >  1:3137476/1‑148 (MQ=255)
                           aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTCGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCt                                                                                                                     >  2:2753888/1‑148 (MQ=255)
                           aaaGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTGTCTTCGTGGAGATTATTCTCTTCCt                                                                                                                     >  1:985113/1‑148 (MQ=255)
                             agcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCa                                                                                                                   <  1:1520094/148‑1 (MQ=255)
                              gcccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCAt                                                                                                                  <  1:1030100/148‑1 (MQ=255)
                               cccagcccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATc                                                                                                                 <  2:2318944/148‑1 (MQ=255)
                                    cccGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCAc                                                                                                            >  1:264280/1‑148 (MQ=255)
                                      cGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACgg                                                                                                          >  2:2273218/1‑148 (MQ=255)
                                      cGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACgg                                                                                                          >  1:397630/1‑148 (MQ=255)
                                        cTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGtt                                                                                                        <  2:3002523/148‑1 (MQ=255)
                                          cACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCGCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCg                                                                                                      >  1:2456001/1‑148 (MQ=255)
                                            ccgccgCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGa                                                                                                    <  1:25289/148‑1 (MQ=255)
                                             cgccgcCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAc                                                                                                   >  2:1576351/1‑148 (MQ=255)
                                             cgccgcCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGAc                                                                                                   >  2:944144/1‑148 (MQ=255)
                                               ccgccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCg                                                                                                 <  2:497049/148‑1 (MQ=255)
                                                  ccAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGaca                                                                                              >  2:2585220/1‑148 (MQ=255)
                                                   cAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGacac                                                                                             >  2:1284652/1‑148 (MQ=255)
                                                   cAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGacac                                                                                             >  2:3060042/1‑148 (MQ=255)
                                                     ggATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACAccc                                                                                           >  2:1004837/1‑148 (MQ=255)
                                                      gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTTCGCTGATCGTATTTATCTTCGTGGAGATTATTATCATCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          >  1:260177/1‑148 (MQ=255)
                                                      gATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCa                                                                                          <  1:3416846/148‑1 (MQ=255)
                                                       aTGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCAc                                                                                         >  2:2302627/1‑148 (MQ=255)
                                                        tGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACg                                                                                        <  1:623270/148‑1 (MQ=255)
                                                         gAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGt                                                                                       <  2:2440686/148‑1 (MQ=255)
                                                          aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCCCGGTTCGGACCGACACCCACGTa                                                                                      <  1:768396/148‑1 (MQ=255)
                                                          aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTa                                                                                      <  1:3475408/148‑1 (MQ=255)
                                                          aaGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTa                                                                                      <  1:1117795/148‑1 (MQ=255)
                                                           aGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAg                                                                                     <  1:3239542/148‑1 (MQ=255)
                                                            gCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGt                                                                                    <  1:1314177/148‑1 (MQ=255)
                                                             cAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCGCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTg                                                                                   <  2:1381641/148‑1 (MQ=255)
                                                                  ggCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTATTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGa                                                                              >  2:1218213/1‑148 (MQ=255)
                                                                     ggTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATAc                                                                           <  1:3106494/148‑1 (MQ=255)
                                                                      gTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACt                                                                          >  2:1142835/1‑148 (MQ=255)
                                                                      gTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACCTAGTGGCAGATACt                                                                          >  1:3048468/1‑148 (MQ=255)
                                                                           ccAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCg                                                                     <  1:3228893/148‑1 (MQ=255)
                                                                             aGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCa                                                                   >  2:1592901/1‑148 (MQ=255)
                                                                             aGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCa                                                                   >  2:401646/1‑148 (MQ=255)
                                                                                ggTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                                >  1:521509/1‑148 (MQ=255)
                                                                                ggTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAg                                                                <  1:2885933/148‑1 (MQ=255)
                                                                                 gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               <  1:1605501/148‑1 (MQ=255)
                                                                                 gTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAgg                                                               <  2:2149955/148‑1 (MQ=255)
                                                                                   tGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCGTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTa                                                             <  2:339655/148‑1 (MQ=255)
                                                                                    gATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAg                                                            <  1:3369768/148‑1 (MQ=255)
                                                                                     aTACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTGTCATTTGCGCTCCTCGTATTTCTCTTAGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTatc                                                           >  1:1910893/1‑146 (MQ=255)
                                                                                      tACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGcc                                                          >  2:491039/1‑148 (MQ=255)
                                                                                           tCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACCGGTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCaaa                                                     >  1:2847881/1‑148 (MQ=255)
                                                                                           tCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTGTCATTTGCGCGCCTCGTATTTCTCGTCGTGCAGATTATTATCTTCCTCCACCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCCAAGGTAGCCGCaaa                                                     >  1:1068584/1‑148 (MQ=255)
                                                                                             cGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGc                                                   <  1:328911/148‑1 (MQ=255)
                                                                                              ggCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGcc                                                  <  1:3055100/148‑1 (MQ=255)
                                                                                                cgcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTGCTGGGATTTCTCTTCGTGCAGATTTTTCTCTTGCTCATCCGCACGGTGCGGGCCGGCACCCCCGGAGTGGCAGAGACTGGGCGCAAAGGTAGCCGCAAAGccca                                                >  1:1009447/1‑146 (MQ=255)
                                                                                                cgcgGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGGAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGGAGCCGCAAAGCCaa                                                >  1:1009451/1‑148 (MQ=255)
                                                                                                    gtgtAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGc                                            <  1:843762/148‑1 (MQ=255)
                                                                                                         gcgcACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGa                                       >  2:2613554/1‑148 (MQ=255)
                                                                                                          cgcACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAg                                      >  2:2255418/1‑148 (MQ=255)
                                                                                                                aTGGTGTCACGTCGTTTCATTTGCGCACCTCGTAATTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGGAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGttt                                >  2:875641/1‑148 (MQ=255)
                                                                                                                  ggTGTCACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTcc                              <  2:1760779/148‑1 (MQ=255)
                                                                                                                    tgtCACGTCGTTTCATTTGGGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAg                            <  1:2116192/148‑1 (MQ=255)
                                                                                                                       cACGTCGTTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATg                         >  2:2664253/1‑148 (MQ=255)
                                                                                                                         cgtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGcc                       <  1:420584/148‑1 (MQ=255)
                                                                                                                          gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCACAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                      >  2:509503/1‑148 (MQ=255)
                                                                                                                          gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                      <  1:401646/148‑1 (MQ=255)
                                                                                                                          gtcgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCt                      <  2:1840945/148‑1 (MQ=255)
                                                                                                                            cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    <  1:2465065/148‑1 (MQ=255)
                                                                                                                            cgtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGGAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCttt                    <  2:384236/148‑1 (MQ=255)
                                                                                                                             gtTTCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTc                   <  1:507182/148‑1 (MQ=255)
                                                                                                                              tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  >  2:2074924/1‑148 (MQ=255)
                                                                                                                              tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  >  2:2075004/1‑148 (MQ=255)
                                                                                                                              tttCATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTcc                  >  2:3460762/1‑148 (MQ=255)
                                                                                                                                 cATTTGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGa               <  1:3515849/148‑1 (MQ=255)
                                                                                                                                    ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            >  2:862087/1‑148 (MQ=255)
                                                                                                                                    ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            >  1:2448840/1‑148 (MQ=255)
                                                                                                                                    ttGCGCTCCTCGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAAtt            >  2:2716536/1‑148 (MQ=255)
                                                                                                                                         ctcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACg       >  2:567802/1‑148 (MQ=255)
                                                                                                                                          tcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAAGAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACgg      <  1:2220959/148‑1 (MQ=255)
                                                                                                                                          tcctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAAGAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACgg      <  1:2307444/148‑1 (MQ=255)
                                                                                                                                           cctcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGt     <  1:731289/148‑1 (MQ=255)
                                                                                                                                             tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCTCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt   <  2:397630/148‑1 (MQ=255)
                                                                                                                                             tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt   >  2:320826/1‑148 (MQ=255)
                                                                                                                                             tcGTATTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGttt   <  1:2753888/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CATCAACCAGTTGAAGGACGGGAACAAAAAGCCCAGCCCGCTCACCGCCGCCAGGATGAAGCAGAAGGCGGTGATCCAGTGGTTGATACGTTCCGGCGCGGTGTAGCGCACGATGGTGTCACGTCGTTTCATTTGCGCTCCTCGTCTTTCTCTTCGTGCAGATTATTCTCTTCCTCATCCGCACGGTTCGGACCGACACCCACGTAGTGGAAGATACTGGCCGCAAAGGTAGCCGCAAAGCCAACAGCCGCGAGCGGTTTCCAGATGCCTTTCCAGAATTTCACGGTTTC  >  CP009273/4075089‑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: ‑