Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,168,795 C→A 73.9% D336Y (GAT→TAT)  gatZ ← D‑tagatose 1,6‑bisphosphate aldolase 2, subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,168,7950CA73.9% 315.3 / 167.1 188D336Y (GAT→TAT) gatZD‑tagatose 1,6‑bisphosphate aldolase 2, subunit
Reads supporting (aligned to +/- strand):  ref base C (17/32);  new base A (76/63);  total (93/95)
Fisher's exact test for biased strand distribution p-value = 1.99e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 191 total aligned reads displayed.

TTTGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAGC  >  CP009273/2168646‑2168942
                                                                                                                                                     |                                                                                                                                                   
tttGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGCATATCCAGTAATGAATc                                                                                                                                                     >  2:2994700/1‑150 (MQ=255)
tttGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATc                                                                                                                                                     >  2:2994703/1‑150 (MQ=255)
  tGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGt                                                                                                                                                   >  1:1786069/1‑150 (MQ=255)
   gggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGtt                                                                                                                                                  <  2:1102708/150‑1 (MQ=255)
   gggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGtt                                                                                                                                                  >  1:1873531/1‑150 (MQ=255)
    ggAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTa                                                                                                                                                 <  1:2557907/150‑1 (MQ=255)
       aGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaaa                                                                                                                                              <  2:3886599/150‑1 (MQ=255)
       aGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGACAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTaaca                                                                                                                                              >  1:3460904/1‑148 (MQ=255)
         aTACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAAcc                                                                                                                                            >  2:3187737/1‑150 (MQ=255)
          tACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAAccc                                                                                                                                           <  1:93573/150‑1 (MQ=255)
           aCTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCg                                                                                                                                          <  1:2831898/150‑1 (MQ=255)
             tGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTa                                                                                                                                        <  2:4485576/150‑1 (MQ=255)
               aCTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACg                                                                                                                                      <  1:4640085/150‑1 (MQ=255)
                cTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGa                                                                                                                                     >  1:849924/1‑150 (MQ=255)
                    tCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataat                                                                                                                                 >  1:4268715/1‑150 (MQ=255)
                     cATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGataata                                                                                                                                <  2:305281/150‑1 (MQ=255)
                     cATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGataata                                                                                                                                <  2:4577867/150‑1 (MQ=255)
                      aTGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAatat                                                                                                                               <  2:2320824/150‑1 (MQ=255)
                      aTGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAatat                                                                                                                               <  2:3702924/150‑1 (MQ=255)
                       tGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtt                                                                                                                              <  1:2022043/150‑1 (MQ=255)
                        gCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATAttt                                                                                                                             <  1:2380875/150‑1 (MQ=255)
                        gCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATAttt                                                                                                                             <  1:4593536/150‑1 (MQ=255)
                        gCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAttt                                                                                                                             <  2:1051681/150‑1 (MQ=255)
                          ccAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATAttttt                                                                                                                           <  2:245383/150‑1 (MQ=255)
                           cAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtttttt                                                                                                                          <  1:1069450/150‑1 (MQ=255)
                           cAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATAtttttt                                                                                                                          <  1:1069385/150‑1 (MQ=255)
                              tGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCa                                                                                                                       <  1:2336036/150‑1 (MQ=255)
                                 gATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTa                                                                                                                    <  1:3428736/150‑1 (MQ=255)
                                    gTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCAGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTATTTCCAGTATTg                                                                                                                 >  2:4687221/1‑150 (MQ=255)
                                      ccACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCg                                                                                                               <  2:2061075/150‑1 (MQ=255)
                                          gCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtc                                                                                                           >  2:257459/1‑150 (MQ=255)
                                          gCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtc                                                                                                           >  2:1142443/1‑150 (MQ=255)
                                           ccTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtcg                                                                                                          <  2:2948952/150‑1 (MQ=255)
                                            cTTCAAGATTCACCCTCATCGTTTCGACGCTATTTTTAATCCGACTCTGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATCTTTTTTCCAGTATTGCGGTtcgt                                                                                                         <  2:2046337/150‑1 (MQ=255)
                                            cTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTtcgt                                                                                                         <  2:1202643/150‑1 (MQ=255)
                                               cAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATGATATTTTTTCCAGTATTGCGGTTCGTCGa                                                                                                      >  1:4287753/1‑150 (MQ=255)
                                                 aGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGc                                                                                                    >  1:1394152/1‑150 (MQ=255)
                                                   aTTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCAt                                                                                                  >  2:3681269/1‑150 (MQ=255)
                                                               cGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAAt                                                                                      <  1:206574/150‑1 (MQ=255)
                                                               cGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAAt                                                                                      <  2:247282/150‑1 (MQ=255)
                                                                  ttCGACGCTATTATTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTAAAAACCCGTACGATAATATATTTTCCAGTATTGCGGGTCGTCGAGCATCACTTCTTCAATTAc                                                                                   >  2:4691295/1‑150 (MQ=255)
                                                                        gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                             >  1:2072344/1‑150 (MQ=255)
                                                                        gCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAg                                                                             >  1:1663664/1‑150 (MQ=255)
                                                                          tATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACTGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGc                                                                           >  2:617868/1‑150 (MQ=255)
                                                                           aTTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCa                                                                          <  2:479200/150‑1 (MQ=255)
                                                                           aTTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCa                                                                          <  2:2102012/150‑1 (MQ=255)
                                                                           aTTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCa                                                                          <  2:1398467/150‑1 (MQ=255)
                                                                              tttAATCCGACTATGCGGCCAGTAATAACGAACACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAAcc                                                                       >  2:2767380/1‑150 (MQ=255)
                                                                                tAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGc                                                                     >  1:4488465/1‑150 (MQ=255)
                                                                                 aaTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCt                                                                    <  1:3083032/150‑1 (MQ=255)
                                                                                 aaTCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCt                                                                    <  1:654426/150‑1 (MQ=255)
                                                                                   tCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGc                                                                  >  2:1431881/1‑150 (MQ=255)
                                                                                     cGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGa                                                                >  2:2853376/1‑150 (MQ=255)
                                                                                      gACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAt                                                               >  1:1710150/1‑150 (MQ=255)
                                                                                        cTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGAttt                                                             <  1:2223153/150‑1 (MQ=255)
                                                                                           tGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCa                                                          <  1:154951/150‑1 (MQ=255)
                                                                                            gCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAg                                                         <  1:407863/150‑1 (MQ=255)
                                                                                             cGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACAGCAAGGCAACCGCTGCGATTTTCAgg                                                        >  1:3657356/1‑150 (MQ=255)
                                                                                               gCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAgggg                                                      <  1:3206278/150‑1 (MQ=255)
                                                                                                  aGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGa                                                   <  1:3416477/150‑1 (MQ=255)
                                                                                                   gTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGAt                                                  <  2:3285214/150‑1 (MQ=255)
                                                                                                    taataaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATa                                                 <  1:3932543/150‑1 (MQ=255)
                                                                                                     aataaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATaa                                                >  1:2467932/1‑150 (MQ=255)
                                                                                                       taaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGt                                              >  1:1199660/1‑150 (MQ=255)
                                                                                                       taaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGt                                              <  2:1105741/150‑1 (MQ=255)
                                                                                                       taaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGt                                              <  2:1105315/150‑1 (MQ=255)
                                                                                                        aaCGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGtt                                             <  2:373622/150‑1 (MQ=255)
                                                                                                          cGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTcc                                           >  1:4831357/1‑150 (MQ=255)
                                                                                                           gAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCt                                          <  2:690641/150‑1 (MQ=255)
                                                                                                           gAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCt                                          <  1:816759/150‑1 (MQ=255)
                                                                                                           gAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCt                                          >  1:3778791/1‑150 (MQ=255)
                                                                                                           gAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCt                                          <  2:2505512/150‑1 (MQ=255)
                                                                                                            aaTACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTg                                         <  1:4142606/150‑1 (MQ=255)
                                                                                                              tACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCt                                       <  2:1853151/150‑1 (MQ=255)
                                                                                                              tACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCt                                       >  2:102308/1‑150 (MQ=255)
                                                                                                               aCGATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTTAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTc                                      <  1:69251/150‑1 (MQ=255)
                                                                                                                 gATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCaa                                    >  1:4062717/1‑150 (MQ=255)
                                                                                                                 gATCCGACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCaa                                    >  1:3959341/1‑150 (MQ=255)
                                                                                                                      gACAGGCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtgt                               <  1:726700/150‑1 (MQ=255)
                                                                                                                       aCAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTtgtg                              >  2:1208552/1‑150 (MQ=255)
                                                                                                                        cAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGc                             <  1:167673/150‑1 (MQ=255)
                                                                                                                          ggCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCa                           >  2:3573873/1‑150 (MQ=255)
                                                                                                                          ggCTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCa                           >  2:4360595/1‑150 (MQ=255)
                                                                                                                            cTGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGt                         >  1:53942/1‑150 (MQ=255)
                                                                                                                             tGTAACGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTg                        >  1:4326817/1‑150 (MQ=255)
                                                                                                                                aaCGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaa                     <  2:783602/150‑1 (MQ=255)
                                                                                                                                aaCGAATATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCaa                     >  1:847904/1‑150 (MQ=255)
                                                                                                                                   gAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAata                  >  1:4032591/1‑150 (MQ=255)
                                                                                                                                    aaTATCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAatat                 >  2:2242368/1‑150 (MQ=255)
                                                                                                                                     atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc                >  2:2255196/1‑150 (MQ=255)
                                                                                                                                     atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc                <  1:1126203/150‑1 (MQ=255)
                                                                                                                                     atatCCAGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATc                >  1:3850713/1‑150 (MQ=255)
                                                                                                                                           aGTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcg          <  2:4326817/150‑1 (MQ=255)
                                                                                                                                            gTAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcgc         >  1:740777/1‑150 (MQ=255)
                                                                                                                                             tAATGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTcgcg        >  1:4001784/1‑150 (MQ=255)
                                                                                                                                                tGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTaa     >  2:4186236/1‑150 (MQ=255)
                                                                                                                                                tGAATAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTaa     >  2:1317475/1‑150 (MQ=255)
                                                                                                                                                  aaTAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAg   >  2:1778816/1‑150 (MQ=255)
                                                                                                                                                  aaTAGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAg   <  2:4773347/150‑1 (MQ=255)
                                                                                                                                                   aTCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAGc  <  2:3909407/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTTGGGAAGATACTGACTAATCATGCCCAGTGGGATGTCCACGCCTTCAAGATTCACCATCATCGTTTCGACGCTATTTTTAATCCGACTATGCGGCCAGTAATAACGAATACGATCCGACAGGCTGTAACGAATATCCAGTAATGAATCGTTAAAACCCGTACGATAATATTTTTTCCAGTATTGCGGTTCGTCGAGCATCACTTCTTCAATTACCGCCAGGCAACCGCTGCGATTTTCAGGGGCGATAAGTTCCTGCTCAATTTGTGCCAGTGCAAATATCGCCTCGCGTAAAGC  >  CP009273/2168646‑2168942

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: