Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,023,757 T→A 42.9% K192I (AAA→ATA)  fadB ← fused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,023,7570TA42.9% 159.1 / inf 289K192I (AAA→ATA) fadBfused 3‑hydroxybutyryl‑CoA epimerase/delta(3)‑cis‑delta(2)‑trans‑enoyl‑CoA isomerase/enoyl‑CoA hydratase/3‑hydroxyacyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (72/93);  new base A (70/54);  total (142/147)
Fisher's exact test for biased strand distribution p-value = 3.31e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

Only 100 of 296 total aligned reads displayed.

GTGCGACCATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCATGATGCCC  >  CP009273/4023608‑4023906
                                                                                                                                                     |                                                                                                                                                     
gtgCGACCATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGtt                                                                                                                                                       <  1:2203868/150‑1 (MQ=255)
    gACCATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTc                                                                                                                                                   >  1:1241505/1‑150 (MQ=255)
       cATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGc                                                                                                                                                <  1:1635223/150‑1 (MQ=255)
       cATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGc                                                                                                                                                >  1:819422/1‑150 (MQ=255)
         tCCCTTTAGCGATGGTGAAGCTCAGGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCTGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAAGGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCtt                                                                                                                                              <  2:220920/150‑1 (MQ=255)
          cccTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCttt                                                                                                                                             <  2:2296359/150‑1 (MQ=255)
             tttAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGGTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGAc                                                                                                                                          <  1:824551/150‑1 (MQ=255)
              ttAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACt                                                                                                                                         <  2:486028/150‑1 (MQ=255)
              ttAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACt                                                                                                                                         >  1:2568358/1‑150 (MQ=255)
                aGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTAc                                                                                                                                       <  2:784246/150‑1 (MQ=255)
                 gCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACg                                                                                                                                      <  2:3915849/150‑1 (MQ=255)
                  cGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGc                                                                                                                                     >  1:432810/1‑150 (MQ=255)
                  cGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGc                                                                                                                                     >  1:771800/1‑150 (MQ=255)
                   gATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGcc                                                                                                                                    >  2:2713921/1‑150 (MQ=255)
                   gATGGTGAAGCTCATGGTGGCTTCAATCGTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGcc                                                                                                                                    <  2:1468579/150‑1 (MQ=255)
                       gTGAAGCTCTTTGTTGCTTCATTCTTTCTCAGTTTTAGTGGTTCCAGCTTCGGCTGAAGTTTTGCTTTCCAGTCGAGGTCGCCGTTATTGGCCTGGCGTAAAACCGCCTTTGCTCCTTCAACCAGTTTTTCTGCTTTGACTACTCCATcc                                                                                                                                <  1:3971010/150‑1 (MQ=255)
                         gAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATccac                                                                                                                              >  1:71794/1‑150 (MQ=255)
                          aaGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATccacc                                                                                                                             >  1:2912154/1‑150 (MQ=255)
                           aGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATccacca                                                                                                                            <  1:3066603/150‑1 (MQ=255)
                           aGCTCATGGTGGCTTCAATCTTGCTCAGATTTAGTGGTTCCAGCTTCGGCTGACGATTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATccacca                                                                                                                            <  2:4466419/150‑1 (MQ=255)
                              tCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCTAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGAc                                                                                                                         >  2:3474197/1‑150 (MQ=255)
                              tCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGGTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGGAATACCGCCTTGGCGCCGTCCACCAGTATTTCTGCTTTGACTACGCCATccaccggcc                                                                                                                         >  2:3472156/1‑146 (MQ=255)
                               cATGGTGGCTTCAATCTTGCTCAGTTTGAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGAcc                                                                                                                        <  2:3928208/150‑1 (MQ=255)
                                    tggCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGAttt                                                                                                                   <  1:3742315/150‑1 (MQ=255)
                                     ggCTTAAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGAtttt                                                                                                                  <  2:3882262/150‑1 (MQ=255)
                                      gCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTc                                                                                                                 <  1:2354060/150‑1 (MQ=255)
                                        ttCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAg                                                                                                               >  2:2390896/1‑150 (MQ=255)
                                        ttCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAg                                                                                                               >  2:2720495/1‑150 (MQ=255)
                                             tCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGTTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCt                                                                                                          <  2:770512/150‑1 (MQ=255)
                                               ttGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGa                                                                                                        >  2:503606/1‑150 (MQ=255)
                                                tGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGAt                                                                                                       >  1:2470855/1‑150 (MQ=255)
                                                tGCTCAGTTTTAGTCGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGAt                                                                                                       >  1:168134/1‑150 (MQ=255)
                                                   tCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCcg                                                                                                    <  2:1638823/150‑1 (MQ=255)
                                                     aGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCcgcg                                                                                                  >  1:1904309/1‑150 (MQ=255)
                                                      gTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCgcgc                                                                                                 >  2:2910738/1‑150 (MQ=255)
                                                         ttAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGa                                                                                              >  1:919491/1‑150 (MQ=255)
                                                         ttAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGa                                                                                              <  2:2881941/150‑1 (MQ=255)
                                                           aGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACa                                                                                            >  1:1938886/1‑150 (MQ=255)
                                                           aGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACa                                                                                            >  2:4190630/1‑150 (MQ=255)
                                                           aGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACa                                                                                            >  2:753759/1‑150 (MQ=255)
                                                             tGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATc                                                                                          <  1:2727909/150‑1 (MQ=255)
                                                              ggTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCt                                                                                         >  2:3095873/1‑150 (MQ=255)
                                                              ggTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCt                                                                                         >  2:3095738/1‑150 (MQ=255)
                                                              ggTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCt                                                                                         <  2:571655/150‑1 (MQ=255)
                                                                       ttCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACcggcg                                                                                <  1:3498608/150‑1 (MQ=255)
                                                                       ttCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACcggcg                                                                                <  1:1731841/150‑1 (MQ=255)
                                                                       ttCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACcggcg                                                                                <  2:455890/150‑1 (MQ=255)
                                                                            cTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAAt                                                                           >  1:4177291/1‑150 (MQ=255)
                                                                              gACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGa                                                                         <  1:1019288/150‑1 (MQ=255)
                                                                                    ttGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCa                                                                   >  1:1554138/1‑150 (MQ=255)
                                                                                     tGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAg                                                                  <  1:3214031/150‑1 (MQ=255)
                                                                                        tttCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAgcgc                                                               >  1:3020590/1‑150 (MQ=255)
                                                                                           ccAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACt                                                            <  2:1051715/150‑1 (MQ=255)
                                                                                             aGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGt                                                          >  2:4218646/1‑150 (MQ=255)
                                                                                              gTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTc                                                         <  2:889338/150‑1 (MQ=255)
                                                                                               tCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCa                                                        <  1:1833725/150‑1 (MQ=255)
                                                                                                cGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAg                                                       >  1:1003167/1‑150 (MQ=255)
                                                                                                 gAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAgc                                                      >  2:200746/1‑150 (MQ=255)
                                                                                                 gAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAgc                                                      >  1:257684/1‑150 (MQ=255)
                                                                                                 gAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAgc                                                      >  2:2450317/1‑150 (MQ=255)
                                                                                                        ccGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGc                                               >  2:4329875/1‑150 (MQ=255)
                                                                                                         cGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCa                                              <  1:1285093/150‑1 (MQ=255)
                                                                                                         cGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCa                                              <  2:831408/150‑1 (MQ=255)
                                                                                                          gTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCAt                                             <  1:4382631/150‑1 (MQ=255)
                                                                                                          gTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCAt                                             >  1:4276866/1‑150 (MQ=255)
                                                                                                          gTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCAt                                             <  2:819482/150‑1 (MQ=255)
                                                                                                           ttAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATa                                            <  2:2110355/150‑1 (MQ=255)
                                                                                                           ttAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGACTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATa                                            <  2:2109575/150‑1 (MQ=255)
                                                                                                             caTGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACg                                          <  2:2369975/149‑1 (MQ=255)
                                                                                                              aTGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGt                                         <  1:718335/150‑1 (MQ=255)
                                                                                                               tGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTg                                        <  1:1564189/150‑1 (MQ=255)
                                                                                                                ggCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTgg                                       >  1:539739/1‑150 (MQ=255)
                                                                                                                  ccTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCa                                     <  1:1364682/150‑1 (MQ=255)
                                                                                                                    tGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATa                                   >  1:3345512/1‑150 (MQ=255)
                                                                                                                     ggCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATAc                                  >  1:603783/1‑150 (MQ=255)
                                                                                                                       cGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGt                                <  2:394185/150‑1 (MQ=255)
                                                                                                                         tAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTAc                              >  1:3811517/1‑150 (MQ=255)
                                                                                                                           aaaCCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAg                            >  2:1671186/1‑150 (MQ=255)
                                                                                                                             aCCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGaa                          >  1:3608526/1‑150 (MQ=255)
                                                                                                                             aCCGCCTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGaa                          >  1:3985757/1‑150 (MQ=255)
                                                                                                                              ccGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCTCCCAGCATACGTGGCATACGTACAGAAc                         >  2:4237454/1‑150 (MQ=255)
                                                                                                                                gCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCg                       <  2:1708361/150‑1 (MQ=255)
                                                                                                                                 ccTTTGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGc                      >  2:3705630/1‑150 (MQ=255)
                                                                                                                                   tttGCGCCTTCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCa                    <  1:3472157/150‑1 (MQ=255)
                                                                                                                                    ttGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCaa                   <  2:3900775/150‑1 (MQ=255)
                                                                                                                                      gcgcCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAg                 <  2:944566/150‑1 (MQ=255)
                                                                                                                                        gcCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGcc               <  2:4438217/150‑1 (MQ=255)
                                                                                                                                        gcCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGcc               <  2:4529669/150‑1 (MQ=255)
                                                                                                                                        gcCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGcc               <  1:472499/150‑1 (MQ=255)
                                                                                                                                           ttCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAgg            <  1:1321797/150‑1 (MQ=255)
                                                                                                                                            tCAACCAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGc           >  1:2902127/1‑150 (MQ=255)
                                                                                                                                               aCCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatg        <  1:2629204/150‑1 (MQ=255)
                                                                                                                                                ccAGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatga       <  2:4149182/150‑1 (MQ=255)
                                                                                                                                                 cAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatgat      <  1:3140162/150‑1 (MQ=255)
                                                                                                                                                 cAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatgat      <  1:857482/150‑1 (MQ=255)
                                                                                                                                                  aGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatgatg     <  2:3885186/150‑1 (MQ=255)
                                                                                                                                                  aGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatgatg     <  2:2583555/150‑1 (MQ=255)
                                                                                                                                                  aGTATTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCatgatg     >  2:1684618/1‑150 (MQ=255)
                                                                                                                                                   gTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCATGATGc    >  1:4134222/1‑150 (MQ=255)
                                                                                                                                                    tttttCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCATGATGcc   <  2:1904309/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTGCGACCATCCCTTTAGCGATGGTGAAGCTCATGGTGGCTTCAATCTTGCTCAGTTTTAGTGGTTCCAGCTTCGGCTGACGTTTTGCTTTCCAGTCGAGGTCGCCGTTAATGGCCTGGCGTAAAACCGCCTTTGCGCCTTCAACCAGTTTTTCTGCTTTGACTACGCCATCCACCAGACCGATTTTCAGCGCCTGATCCGCGCCGACATCTTTACCGGCGGCAATGATTTCCAGCGCACTGTCAGCGCCCAGCATACGTGGCATACGTACAGAACCGCCAAAGCCAGGCATGATGCCC  >  CP009273/4023608‑4023906

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: