Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,576,325 C→A 100% E171D (GAG→GAT)  nudK ← GDP‑mannose pyrophosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,576,3250CA99.2% 857.9 / ‑1.7 244E171D (GAG→GAT) nudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/2);  new base A (124/116);  total (124/120)
Fisher's exact test for biased strand distribution p-value = 2.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-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 244 total aligned reads displayed.

ATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCG  >  CP009273/2576180‑2576469
                                                                                                                                                 |                                                                                                                                                
    tCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGt                                                                                                                                          >  1:3193956/1‑150 (MQ=255)
       ttGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATACTAACCAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGtttt                                                                                                                                       <  1:3865350/150‑1 (MQ=255)
          cgcAGGTATCGCGGCTCAATCTGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGat                                                                                                                                    <  1:2120920/150‑1 (MQ=255)
          cgcAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGat                                                                                                                                    <  1:3877016/150‑1 (MQ=255)
            cAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatca                                                                                                                                  >  2:210385/1‑150 (MQ=255)
             aGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatcat                                                                                                                                 <  1:4140429/150‑1 (MQ=255)
               gTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtct                                                                                                                               >  1:3022882/1‑150 (MQ=255)
                tAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtctc                                                                                                                              >  2:2842647/1‑150 (MQ=255)
                  gcgcGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCa                                                                                                                            <  2:3682528/150‑1 (MQ=255)
                  gcgcGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCa                                                                                                                            <  2:2285574/150‑1 (MQ=255)
                   cgcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAg                                                                                                                           <  2:189926/150‑1 (MQ=255)
                         tCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTg                                                                                                                     <  2:3513238/150‑1 (MQ=255)
                            aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGACATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                  <  1:1639431/150‑1 (MQ=255)
                            aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                  >  2:4098253/1‑150 (MQ=255)
                            aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                  <  1:3512190/150‑1 (MQ=255)
                            aTCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                  >  1:2410191/1‑150 (MQ=255)
                             tCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTg                                                                                                                 <  2:394987/150‑1 (MQ=255)
                               ggATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGaa                                                                                                               <  2:413475/150‑1 (MQ=255)
                               ggATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGaa                                                                                                               <  1:1800915/150‑1 (MQ=255)
                                 aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                             <  1:2943645/150‑1 (MQ=255)
                                 aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                             <  1:911729/150‑1 (MQ=255)
                                    tGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCa                                                                                                          <  2:4176047/150‑1 (MQ=255)
                                    tGTAGGATTGAGCCAATAATAAACCACGAAATATTATTGTTAATTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCTCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCa                                                                                                          <  2:1896850/150‑1 (MQ=255)
                                     gTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAg                                                                                                         >  2:351210/1‑150 (MQ=255)
                                       cGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCt                                                                                                       >  2:1035084/1‑150 (MQ=255)
                                       cGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCt                                                                                                       >  2:814700/1‑150 (MQ=255)
                                            tGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGc                                                                                                  >  1:1514915/1‑150 (MQ=255)
                                                ccAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACtt                                                                                              >  1:2555134/1‑150 (MQ=255)
                                                 cAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTc                                                                                             >  2:1736465/1‑150 (MQ=255)
                                                   ataataAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCaa                                                                                           >  2:4054558/1‑150 (MQ=255)
                                                    taataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAat                                                                                          >  1:136480/1‑150 (MQ=255)
                                                      ataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAatat                                                                                        <  1:1634787/150‑1 (MQ=255)
                                                          acaacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttc                                                                                    <  1:3219370/150‑1 (MQ=255)
                                                              cGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatct                                                                                >  2:2770371/1‑150 (MQ=255)
                                                              cGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatat                                                                                >  2:524021/1‑148 (MQ=255)
                                                               gAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTAACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCGGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatctt                                                                               >  2:3523882/1‑150 (MQ=255)
                                                                 aaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCg                                                                             >  1:714416/1‑150 (MQ=255)
                                                                   tattatTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGac                                                                           >  2:2386553/1‑150 (MQ=255)
                                                                   tattatTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGac                                                                           >  2:2218862/1‑150 (MQ=255)
                                                                   tattatTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGac                                                                           <  1:3355861/150‑1 (MQ=255)
                                                                    attattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGaca                                                                          >  1:2484278/1‑150 (MQ=255)
                                                                      tattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACAcc                                                                        <  1:4008450/150‑1 (MQ=255)
                                                                      tattGTTATTTTTCAGACCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACAGTCACGTATATCTCCGGTTATGATCATCTCCAGCGCATGTATGAACGTAAGCTCGAGCACTTAAATATCTTCATATTCGacaac                                                                        >  2:2355759/1‑148 (MQ=255)
                                                                       attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                       <  2:2692417/150‑1 (MQ=255)
                                                                       attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                       <  2:1937490/150‑1 (MQ=255)
                                                                        ttGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGc                                                                      >  1:2389/1‑150 (MQ=255)
                                                                        ttGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGc                                                                      >  1:1288169/1‑150 (MQ=255)
                                                                        ttGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGc                                                                      >  2:879453/1‑150 (MQ=255)
                                                                          gTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCa                                                                    >  2:1033580/1‑150 (MQ=255)
                                                                          gTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCa                                                                    >  2:3112900/1‑150 (MQ=255)
                                                                            tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                  >  2:642459/1‑150 (MQ=255)
                                                                             aTTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGGTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCg                                                                 <  1:2763375/150‑1 (MQ=255)
                                                                               ttttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGc                                                               >  1:3157185/1‑150 (MQ=255)
                                                                                 ttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGt                                                             <  1:3644627/150‑1 (MQ=255)
                                                                                  tCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGtt                                                            <  1:3775703/150‑1 (MQ=255)
                                                                                    aGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAg                                                          <  1:1145046/150‑1 (MQ=255)
                                                                                      tCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcg                                                        >  2:1761501/1‑150 (MQ=255)
                                                                                       ccATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgc                                                       >  1:2749065/1‑150 (MQ=255)
                                                                                         aTTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgcgc                                                     >  2:3991008/1‑150 (MQ=255)
                                                                                         aTTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAgcgcgc                                                     >  1:373604/1‑150 (MQ=255)
                                                                                             aaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAt                                                 <  2:454194/150‑1 (MQ=255)
                                                                                              aTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAtt                                                <  1:2123286/150‑1 (MQ=255)
                                                                                               tgtgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTg                                               >  1:4096680/1‑150 (MQ=255)
                                                                                                 tgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTc                                             >  2:2734678/1‑150 (MQ=255)
                                                                                                    cGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACt                                          <  2:3989075/150‑1 (MQ=255)
                                                                                                    cGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACt                                          <  1:3328464/150‑1 (MQ=255)
                                                                                                     gTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTg                                         >  1:2593725/1‑150 (MQ=255)
                                                                                                       ttGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTa                                       >  2:3368653/1‑150 (MQ=255)
                                                                                                        tGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTAt                                      >  2:1128617/1‑150 (MQ=255)
                                                                                                          cAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTc                                    <  1:2872642/150‑1 (MQ=255)
                                                                                                           aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                   <  2:1391714/150‑1 (MQ=255)
                                                                                                            aaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCgg                                  >  2:1829125/1‑150 (MQ=255)
                                                                                                             aTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGc                                 <  1:4147423/150‑1 (MQ=255)
                                                                                                             aTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGc                                 >  2:3789080/1‑150 (MQ=255)
                                                                                                              gaGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCg                                <  1:1334499/149‑1 (MQ=255)
                                                                                                               aGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGa                               <  2:2421789/150‑1 (MQ=255)
                                                                                                                gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                              >  1:977596/1‑150 (MQ=255)
                                                                                                                gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                              <  2:1501000/150‑1 (MQ=255)
                                                                                                                 ttAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATa                             >  2:1425272/1‑150 (MQ=255)
                                                                                                                  tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                            >  2:3038002/1‑150 (MQ=255)
                                                                                                                  tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                            >  2:1468045/1‑150 (MQ=255)
                                                                                                                   aaGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaa                           >  2:3604209/1‑150 (MQ=255)
                                                                                                                      cAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaaa                        >  1:3991661/1‑150 (MQ=255)
                                                                                                                           acacCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGa                   >  2:1579423/1‑150 (MQ=255)
                                                                                                                             acCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                 <  2:2163265/150‑1 (MQ=255)
                                                                                                                             acCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                 <  2:1154359/150‑1 (MQ=255)
                                                                                                                             acCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                 <  2:3391841/150‑1 (MQ=255)
                                                                                                                              ccGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGTGATAAAAAAGTGGATTa                <  1:1099404/150‑1 (MQ=255)
                                                                                                                              ccGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTa                <  1:246399/150‑1 (MQ=255)
                                                                                                                                gTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGc              >  2:106656/1‑150 (MQ=255)
                                                                                                                                 tCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCt             <  1:3646353/150‑1 (MQ=255)
                                                                                                                                  cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc            >  1:4114753/1‑150 (MQ=255)
                                                                                                                                  cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc            >  1:1215038/1‑150 (MQ=255)
                                                                                                                                  cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc            <  1:2080058/150‑1 (MQ=255)
                                                                                                                                  cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc            <  1:472737/150‑1 (MQ=255)
                                                                                                                                   ttACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCa           <  1:3162308/150‑1 (MQ=255)
                                                                                                                                      ccGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTc        <  1:225839/150‑1 (MQ=255)
                                                                                                                                        gTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac      <  1:204715/150‑1 (MQ=255)
                                                                                                                                          cACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac    <  2:2541672/150‑1 (MQ=255)
                                                                                                                                           aCGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCGGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACAcc   <  1:761139/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
ATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCG  >  CP009273/2576180‑2576469

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