Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,576,325 C→A 92.1% 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,3250CA92.1% 594.3 / 38.6 190E171D (GAG→GATnudKGDP‑mannose pyrophosphatase
Reads supporting (aligned to +/- strand):  ref base C (4/11);  new base A (89/86);  total (93/97)
Fisher's exact test for biased strand distribution p-value = 1.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.97e-01

Only 100 of 194 total aligned reads displayed.

CGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576177‑2576471
                                                                                                                                                    |                                                                                                                                                  
 gTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTAtctc                                                                                                                                                  <  1:4849316/150‑1 (MQ=255)
  tATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCg                                                                                                                                                 >  1:2167905/1‑150 (MQ=255)
    tCTTCGTTGCGCAGGTCGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGcc                                                                                                                                               >  1:3950389/1‑150 (MQ=255)
      ttCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCgg                                                                                                                                             <  1:1007655/150‑1 (MQ=255)
        cGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGtt                                                                                                                                           <  1:4309863/150‑1 (MQ=255)
         gTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGttt                                                                                                                                          >  2:2071379/1‑150 (MQ=255)
          ttGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGtttt                                                                                                                                         <  1:3733822/150‑1 (MQ=255)
          ttGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGtttt                                                                                                                                         <  1:4143237/150‑1 (MQ=255)
           tGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTg                                                                                                                                        <  1:4067096/150‑1 (MQ=255)
               cAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGatca                                                                                                                                    <  2:3239545/150‑1 (MQ=255)
                aGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGatcat                                                                                                                                   <  2:2449188/150‑1 (MQ=255)
                  gTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtct                                                                                                                                 <  1:159171/150‑1 (MQ=255)
                  gTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCAtct                                                                                                                                 <  1:159090/150‑1 (MQ=255)
                      cgcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAg                                                                                                                             <  1:422295/150‑1 (MQ=255)
                      cgcgGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAg                                                                                                                             <  1:3055902/150‑1 (MQ=255)
                         ggCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAgcgc                                                                                                                          >  2:180694/1‑150 (MQ=255)
                         ggCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAgcgc                                                                                                                          >  2:2794762/1‑150 (MQ=255)
                          gCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGcc                                                                                                                         >  1:1075610/1‑150 (MQ=255)
                          gCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGcc                                                                                                                         <  1:233725/150‑1 (MQ=255)
                               aTCGGATTTGTCGGATTGAGACAATAATAAAAAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCt                                                                                                                    <  2:2945873/150‑1 (MQ=255)
                                tCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTg                                                                                                                   <  2:2167905/150‑1 (MQ=255)
                                   gATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAAc                                                                                                                <  2:792829/150‑1 (MQ=255)
                                    aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                               <  1:3059291/150‑1 (MQ=255)
                                    aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                               <  1:1744861/150‑1 (MQ=255)
                                    aTTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACg                                                                                                               <  2:3272126/150‑1 (MQ=255)
                                        gTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAg                                                                                                           <  2:2538373/150‑1 (MQ=255)
                                         tCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGc                                                                                                          >  2:818646/1‑150 (MQ=255)
                                           ggATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTc                                                                                                        >  2:2426544/1‑150 (MQ=255)
                                            gATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCg                                                                                                       <  1:3264973/150‑1 (MQ=255)
                                                  gCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACt                                                                                                 >  2:4168403/1‑150 (MQ=255)
                                                    cAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTc                                                                                               >  1:2092062/1‑150 (MQ=255)
                                                     aataatAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCa                                                                                              >  1:2523943/1‑150 (MQ=255)
                                                       taataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAat                                                                                            <  2:431784/150‑1 (MQ=255)
                                                       taataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAat                                                                                            <  2:3622728/150‑1 (MQ=255)
                                                        aataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAata                                                                                           >  1:2862789/1‑150 (MQ=255)
                                                         ataaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAatat                                                                                          >  1:89151/1‑150 (MQ=255)
                                                          taaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatc                                                                                         >  1:4537455/1‑150 (MQ=255)
                                                          taaACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatc                                                                                         <  1:582998/150‑1 (MQ=255)
                                                            aacaacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatctt                                                                                       >  2:4930639/1‑150 (MQ=255)
                                                              caacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttca                                                                                     <  2:5030561/150‑1 (MQ=255)
                                                               aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATATTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTATTGATAATCTCCAGCGACAGTCTGAACTGCAGCTCGAGCACTTCAAAatcttcat                                                                                    >  1:3697466/1‑150 (MQ=255)
                                                               aacGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcat                                                                                    <  1:1819240/150‑1 (MQ=255)
                                                                acGAAATATTATTGTTATTTTTCAGTCCATTACATGTGACGTTTGCAAATAGTTAAGCAAGAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   <  2:855343/150‑1 (MQ=255)
                                                                acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   <  1:4298578/150‑1 (MQ=255)
                                                                acGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatc                                                                                   <  2:3742391/150‑1 (MQ=255)
                                                                 cGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatct                                                                                  >  2:2606032/1‑150 (MQ=255)
                                                                   aaaTATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATatcttcatcttc                                                                                >  2:274424/1‑150 (MQ=255)
                                                                         tattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACAcc                                                                          >  2:4088264/1‑150 (MQ=255)
                                                                         tattGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAAGAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACAcc                                                                          <  2:1424967/150‑1 (MQ=255)
                                                                          attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                         >  2:615133/1‑150 (MQ=255)
                                                                          attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                         >  2:3382369/1‑150 (MQ=255)
                                                                          attGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCg                                                                         <  2:3505448/150‑1 (MQ=255)
                                                                             gTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCa                                                                      >  1:313724/1‑150 (MQ=255)
                                                                              ttATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAc                                                                     <  2:1782199/150‑1 (MQ=255)
                                                                               tATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCAcc                                                                    >  1:2802294/1‑150 (MQ=255)
                                                                                  ttttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGc                                                                 >  2:3355489/1‑150 (MQ=255)
                                                                                  ttttCAGTCCATTAAATGTGAAGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGc                                                                 <  2:4170753/150‑1 (MQ=255)
                                                                                   tttCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCg                                                                <  1:91153/150‑1 (MQ=255)
                                                                                     tCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGtt                                                              >  2:4785430/1‑150 (MQ=255)
                                                                                     tCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGtt                                                              <  2:4968594/150‑1 (MQ=255)
                                                                                      cAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTa                                                             <  2:4513867/150‑1 (MQ=255)
                                                                                              tAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTg                                                     >  1:4754330/1‑150 (MQ=255)
                                                                                                aaTGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGAt                                                   >  2:106344/1‑150 (MQ=255)
                                                                                                    tgACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTc                                               >  2:1727396/1‑150 (MQ=255)
                                                                                                     gACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCa                                              >  2:87980/1‑150 (MQ=255)
                                                                                                       cGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACt                                            <  1:1314292/150‑1 (MQ=255)
                                                                                                       cGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACt                                            <  1:4040917/150‑1 (MQ=255)
                                                                                                       cGTTTGCAAAGAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACt                                            <  1:4956398/150‑1 (MQ=255)
                                                                                                          ttGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACCCCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTa                                         >  2:947830/1‑150 (MQ=255)
                                                                                                            gCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTAtt                                       >  2:2281394/1‑150 (MQ=255)
                                                                                                              aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                     >  2:717154/1‑150 (MQ=255)
                                                                                                              aaaTAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCg                                     >  2:4270288/1‑150 (MQ=255)
                                                                                                                 tAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCg                                  <  2:3895642/150‑1 (MQ=255)
                                                                                                                 tAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCg                                  <  1:1441537/150‑1 (MQ=255)
                                                                                                                 tAGTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCg                                  <  2:1271057/150‑1 (MQ=255)
                                                                                                                   gTTAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGAt                                <  2:4673222/150‑1 (MQ=255)
                                                                                                                     tAAGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaa                              >  1:904744/1‑150 (MQ=255)
                                                                                                                       aGCAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaa                            >  2:3668512/1‑150 (MQ=255)
                                                                                                                         cAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaaa                          >  1:3215935/1‑150 (MQ=255)
                                                                                                                         cAATAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATaaaaaa                          >  1:3215887/1‑150 (MQ=255)
                                                                                                                          aaTAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAg                         >  2:3556607/1‑150 (MQ=255)
                                                                                                                          aaTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAg                         >  2:2066450/1‑150 (MQ=255)
                                                                                                                           aTAACACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGt                        <  2:2373746/150‑1 (MQ=255)
                                                                                                                             aaCACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTgg                      <  2:3014943/150‑1 (MQ=255)
                                                                                                                             aaCACCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTgg                      <  1:2019188/150‑1 (MQ=255)
                                                                                                                               cacCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    >  1:2307876/1‑150 (MQ=255)
                                                                                                                               cacCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAt                    <  2:4300522/150‑1 (MQ=255)
                                                                                                                                acCGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGAtt                   <  1:1227919/150‑1 (MQ=255)
                                                                                                                                 ccGTCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCCGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTa                  >  2:3933142/1‑150 (MQ=255)
                                                                                                                                   gTCTTACCGTCACGTATATCGCAGGTTTTGAACATCTCCAGCGCCTGGCTGAACGGCCGCACGAGCACTTCAATATCTTAATCTTCGACACCGCCACCGGCGTTAGCGCGCAGATTGTAACAGTATTCGGCGATAAAAAAGTGGATTAga                >  1:973609/1‑149 (MQ=255)
                                                                                                                                    tCTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCt               <  2:4221033/150‑1 (MQ=255)
                                                                                                                                     cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              >  2:2845449/1‑150 (MQ=255)
                                                                                                                                     cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              >  1:2902207/1‑150 (MQ=255)
                                                                                                                                     cTTACCGTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTc              >  2:4162528/1‑150 (MQ=255)
                                                                                                                                           gTCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcac        <  1:2426544/150‑1 (MQ=255)
                                                                                                                                            tCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTACATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcaca       <  1:3350908/150‑1 (MQ=255)
                                                                                                                                            tCACGTATATCGCCGGTTTTGATCATCTCCAGCGCCGGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTcaca       <  1:603325/150‑1 (MQ=255)
                                                                                                                                             cACGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCacac      <  2:4466310/150‑1 (MQ=255)
                                                                                                                                               cGTATATCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCg    <  1:4374881/150‑1 (MQ=255)
                                                                                                                                                 tataTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGcc  <  2:2388624/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
CGTATCTTCGTTGCGCAGGTAGCGCGGCTCAATCGGATTTGTCGGATTGAGCCAATAATAAACAACGAAATATTATTGTTATTTTTCAGTCCATTAAATGTGACGTTTGCAAATAGTTAAGCAATAACACCGTCTTACCGTCACGTATCTCGCCGGTTTTGATCATCTCCAGCGCCTGGCTGAACGGCAGCTCGAGCACTTCAATATCTTCATCTTCGACACCGCCACCGGCGTTAGCGCGCTGATTGTCACTGTATTCGGCGATAAAAAAGTGGATTAGCTCAGTCACACCGCC  >  CP009273/2576177‑2576471

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: