Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,343,941 C→A 27.7% D184Y (GAT→TAT)  glpQ ← periplasmic glycerophosphodiester phosphodiesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,343,9410CA27.7% 402.6 / 153.4 188D184Y (GAT→TAT) glpQperiplasmic glycerophosphodiester phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base C (65/71);  new base A (17/35);  total (82/106)
Fisher's exact test for biased strand distribution p-value = 7.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.64e-01

Only 100 of 212 total aligned reads displayed.

ACCAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCC  >  CP009273/2343793‑2344090
                                                                                                                                                    |                                                                                                                                                     
aCCAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATcc                                                                                                                                                      <  2:2411050/150‑1 (MQ=255)
aCCAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATcc                                                                                                                                                      >  1:2859818/1‑150 (MQ=255)
  cAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATcctt                                                                                                                                                    <  2:3695520/150‑1 (MQ=255)
  cAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATactt                                                                                                                                                    <  2:4383293/150‑1 (MQ=255)
  cAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATactt                                                                                                                                                    >  2:4011563/1‑150 (MQ=255)
   aGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATccttc                                                                                                                                                   >  2:1471730/1‑150 (MQ=255)
   aGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATccttc                                                                                                                                                   <  1:3214469/150‑1 (MQ=255)
    gATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATacttcc                                                                                                                                                  <  2:1406868/150‑1 (MQ=255)
     aTTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATacttccc                                                                                                                                                 <  1:1059758/150‑1 (MQ=255)
        gAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATccttcccttc                                                                                                                                              <  1:4819908/150‑1 (MQ=255)
              cATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCtgatg                                                                                                                                        >  2:1130965/1‑150 (MQ=255)
               aTGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCtgatga                                                                                                                                       <  2:1009837/150‑1 (MQ=255)
                tGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTgatgat                                                                                                                                      <  1:2660667/150‑1 (MQ=255)
                tGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTgatgat                                                                                                                                      <  1:2659708/150‑1 (MQ=255)
                  cccATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATgg                                                                                                                                    <  2:767317/150‑1 (MQ=255)
                  cccATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATgg                                                                                                                                    >  1:426110/1‑150 (MQ=255)
                   ccATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGa                                                                                                                                   >  2:1322018/1‑150 (MQ=255)
                    cATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACAGGTGTAACCATATTTCTTCAGCACTACCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGaa                                                                                                                                  >  1:1842553/1‑150 (MQ=255)
                            gTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGgcg                                                                                                                          >  1:3651994/1‑150 (MQ=255)
                             ttCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGgcgc                                                                                                                         >  1:3542819/1‑150 (MQ=255)
                              tCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCt                                                                                                                        <  1:4690159/150‑1 (MQ=255)
                              tCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCt                                                                                                                        <  2:684162/150‑1 (MQ=255)
                               ccAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCtt                                                                                                                       <  1:4161875/150‑1 (MQ=255)
                                 aGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTg                                                                                                                     <  2:3701617/150‑1 (MQ=255)
                                   cTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCACTATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGAt                                                                                                                   >  1:1724715/1‑150 (MQ=255)
                                      aTTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTc                                                                                                                >  2:3097297/1‑150 (MQ=255)
                                      aTTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTc                                                                                                                <  2:386008/150‑1 (MQ=255)
                                       ttCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCt                                                                                                               >  1:2955126/1‑150 (MQ=255)
                                       ttCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCt                                                                                                               <  2:4072718/150‑1 (MQ=255)
                                        tCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTg                                                                                                              >  2:264211/1‑150 (MQ=255)
                                           tAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGAt                                                                                                           <  1:1996729/150‑1 (MQ=255)
                                           tAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGAt                                                                                                           <  1:1997072/150‑1 (MQ=255)
                                            aaTACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATa                                                                                                          <  1:880066/150‑1 (MQ=255)
                                            aaTACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATa                                                                                                          <  1:3176702/150‑1 (MQ=255)
                                               aCGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAAt                                                                                                       <  1:979831/150‑1 (MQ=255)
                                               aCGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAAt                                                                                                       <  1:2161930/150‑1 (MQ=255)
                                               aCGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAAt                                                                                                       <  1:219948/150‑1 (MQ=255)
                                                cGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTATACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACAACGGCGCTTTGATTTCTGGATAAATa                                                                                                      >  1:1312734/1‑150 (MQ=255)
                                                  cTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATAcc                                                                                                    <  1:66716/150‑1 (MQ=255)
                                                   ttCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCg                                                                                                   <  2:1813063/150‑1 (MQ=255)
                                                   ttCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCg                                                                                                   >  1:3784778/1‑150 (MQ=255)
                                                    tCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGa                                                                                                  <  1:3945183/150‑1 (MQ=255)
                                                      aGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGata                                                                                                >  1:4602101/1‑150 (MQ=255)
                                                      aGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGata                                                                                                <  1:1492967/150‑1 (MQ=255)
                                                      aGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGata                                                                                                <  1:4731640/150‑1 (MQ=255)
                                                        cTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATAtt                                                                                              >  1:3556877/1‑150 (MQ=255)
                                                         tcatcaGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCGGTGATTTATGGATAAATAACGATAttt                                                                                             >  1:4572038/1‑150 (MQ=255)
                                                          catcaGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATAtttt                                                                                            <  2:4861865/150‑1 (MQ=255)
                                                          catcaGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATAAAGATAtttt                                                                                            >  1:3291487/1‑150 (MQ=255)
                                                           atcaGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATAttttt                                                                                           <  1:2310064/150‑1 (MQ=255)
                                                             caGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTcc                                                                                         <  2:630097/150‑1 (MQ=255)
                                                              aGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTccc                                                                                        <  1:3243968/150‑1 (MQ=255)
                                                               gCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCg                                                                                       <  1:739745/150‑1 (MQ=255)
                                                                cATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCgg                                                                                      <  2:3235974/150‑1 (MQ=255)
                                                                  tCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTa                                                                                    <  1:3196738/150‑1 (MQ=255)
                                                                   cAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTag                                                                                   <  2:2734759/150‑1 (MQ=255)
                                                                    aaaaCATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTaga                                                                                  <  1:680936/150‑1 (MQ=255)
                                                                      aaCATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGt                                                                                <  1:365828/150‑1 (MQ=255)
                                                                      aaCATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGt                                                                                <  1:2596251/150‑1 (MQ=255)
                                                                      aaCATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGt                                                                                <  2:4457265/150‑1 (MQ=255)
                                                                           tGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGAttt                                                                           <  2:4043365/150‑1 (MQ=255)
                                                                            gCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTa                                                                          <  1:1054819/150‑1 (MQ=255)
                                                                             cAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTaa                                                                         <  1:2773469/150‑1 (MQ=255)
                                                                             cAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTaa                                                                         >  2:2971931/1‑150 (MQ=255)
                                                                                   aaCTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACAGATATTTTTCCCGGTAGAGTGATTTAACCCCTg                                                                   >  2:3476218/1‑150 (MQ=255)
                                                                                    aCTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGa                                                                  <  1:3080981/150‑1 (MQ=255)
                                                                                      tttATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAAc                                                                <  1:3634470/150‑1 (MQ=255)
                                                                                        tATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACaa                                                              >  2:384318/1‑150 (MQ=255)
                                                                                            gTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTc                                                          <  2:2498507/150‑1 (MQ=255)
                                                                                            gTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTc                                                          <  2:2934048/150‑1 (MQ=255)
                                                                                              cTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGAAATTTTTCCCGGTAGAGTGATTTAAACCCTGAACAAATTCaa                                                        >  2:285348/1‑150 (MQ=255)
                                                                                               tttACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAAt                                                       <  1:2330231/150‑1 (MQ=255)
                                                                                               tttACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAAt                                                       <  2:4595833/150‑1 (MQ=255)
                                                                                                 tACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAAtct                                                     <  2:4423666/150‑1 (MQ=255)
                                                                                                 tACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGTTTTCTGGATAAATAACGATATTTTTACCGGTAGAGTGATTTAAACCCAGAAAAAATTAAAtat                                                     >  1:190635/1‑148 (MQ=255)
                                                                                                  aCCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCGTTGATTTCTGGGTAAATACCCATATTTTTACCGGTAGAGTGATTTAACCCCAGAACAAATTCAAactc                                                    >  1:4145631/1‑150 (MQ=255)
                                                                                                   ccGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATctct                                                   <  1:761682/150‑1 (MQ=255)
                                                                                                     ggTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCtcttc                                                 >  2:456562/1‑150 (MQ=255)
                                                                                                      gtgtAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGTTTTCTGGGTTAAAAACGATATTTTTGCCGGTAGAGGGATTTAAACCCGGTACCAATTCAAtagcgtat                                                >  1:3362762/1‑143 (MQ=255)
                                                                                                      gtgtAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCtcttct                                                <  1:4423302/150‑1 (MQ=255)
                                                                                                        gtAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCttcttc                                              <  1:4331020/150‑1 (MQ=255)
                                                                                                        gtAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCttcttc                                              <  1:4404873/150‑1 (MQ=255)
                                                                                                              atatTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTTATGATGGAAGCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTGTAGTGATTTAACCCCAGAACAAAGTCAATCTATTATTCAAAgtt                                        >  2:3489339/1‑148 (MQ=255)
                                                                                                               tatTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGgtg                                       >  1:1112504/1‑150 (MQ=255)
                                                                                                                    cttcAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTAAAAGGTGTGCAc                                  >  1:2581252/1‑150 (MQ=255)
                                                                                                                     ttcAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCAcc                                 >  2:2802766/1‑150 (MQ=255)
                                                                                                                      tcAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCAccc                                >  2:2708929/1‑150 (MQ=255)
                                                                                                                       cAGCACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTTATTTATTGATAAATAACGATATTTTTACCCGTATATTGATTTAAACCCATAACAAAATAAAACTATTATTAAAAGGTGTGGAcccc                               >  1:4510377/1‑149 (MQ=255)
                                                                                                                        aGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAAAAACGATAATTTTCCCGGTAGAGTGATTTAACCCCAGAACAAAATCAAACACATCGTCAAAGGTGTGCACCCgg                              >  1:2522193/1‑150 (MQ=255)
                                                                                                                        aGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCGTTGATTTGTGGGTAAATACCCATATTTTTAACGGTAGAGTGATTTTACCCCAGAACAAAATCAATCTCTTCATAAAAGGTGTTCAccccg                              >  1:2839685/1‑148 (MQ=255)
                                                                                                                          cACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGaa                            <  2:3264645/150‑1 (MQ=255)
                                                                                                                          cACTTCCAGCGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGaa                            >  1:1547489/1‑150 (MQ=255)
                                                                                                                            cTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGTTTTCTGGATAAATACCGATATTTTTACCGGTAGAGTGTTTTAACCCCAGAACAAATTCAAACTCTTCTTCAAAGGTGTGCACCCGGAAGt                          >  1:3866915/1‑150 (MQ=255)
                                                                                                                             ttCCAGCGTTTGTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTc                         <  2:3382124/150‑1 (MQ=255)
                                                                                                                                   cGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACtt                   <  2:4602101/150‑1 (MQ=255)
                                                                                                                                   cGTTTTTGCCGCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACtt                   <  2:426110/150‑1 (MQ=255)
                                                                                                                                             gCAATATACTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTggg         >  2:192553/1‑150 (MQ=255)
                                                                                                                                              cAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGa        <  2:3556877/150‑1 (MQ=255)
                                                                                                                                               aaTATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGaa       <  2:2356661/150‑1 (MQ=255)
                                                                                                                                                  atCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACg    >  2:3465881/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
ACCAGATTGAGCTCCATGCCCATTTTGGGTTCCAGCTCATTCTTAATACGCTTCAGCTCATCAGCATCAAAACATTGCAAATAAACTTTATCGTCTTTACCGGTGTAACCATATTTCTTCAGCACTTCCAGCGTTTTTGCCGCAATATCCTTCCCTTCCTGATGATGGAACCACGGCGCTTTGATTTCTGGATAAATACCGATATTTTTCCCGGTAGAGTGATTTAACCCCTGAACAAATTCAATCTCTTCTTCAAAGGTGTGCACCCGGAAGTCGGACTTACCCATTGGGAAACGCC  >  CP009273/2343793‑2344090

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: