Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,892,304 C→A 100% P150P (CCG→CCT yieK ← putative 6‑phosphogluconolactonase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,892,3040CA98.1% 1016.6 / 7.4 327P150P (CCG→CCTyieKputative 6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (179/141);  total (181/146)
Fisher's exact test for biased strand distribution p-value = 2.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 327 total aligned reads displayed.

CGCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAGT  >  CP009273/3892157‑3892450
                                                                                                                                                   |                                                                                                                                                  
  cTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGa                                                                                                                                                >  2:3297572/1‑150 (MQ=255)
    ttACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCCTCTCTCCCTGGATAGGGAAc                                                                                                                                              >  2:147594/1‑150 (MQ=255)
       cccGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTcc                                                                                                                                           <  1:4324751/150‑1 (MQ=255)
         cgcgCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCAc                                                                                                                                         <  2:4277853/150‑1 (MQ=255)
         cgcgCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCAc                                                                                                                                         <  2:1987315/150‑1 (MQ=255)
          gcgcCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACg                                                                                                                                        <  2:1879340/150‑1 (MQ=255)
          gcgcCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACg                                                                                                                                        >  1:2102475/1‑150 (MQ=255)
            gccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGt                                                                                                                                      >  1:1565435/1‑150 (MQ=255)
             ccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGtt                                                                                                                                     <  2:461289/150‑1 (MQ=255)
             ccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGtt                                                                                                                                     <  1:1040328/150‑1 (MQ=255)
                cTGACAATTATCAGCTAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGc                                                                                                                                  <  2:3514847/150‑1 (MQ=255)
                 tGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCACCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                 >  1:4510086/1‑150 (MQ=255)
                 tGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                 <  2:2191648/150‑1 (MQ=255)
                  gACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTc                                                                                                                                >  1:1719721/1‑150 (MQ=255)
                    cAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCAt                                                                                                                              <  2:2657759/150‑1 (MQ=255)
                      aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATgg                                                                                                                            <  1:3855676/150‑1 (MQ=255)
                      aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATgg                                                                                                                            <  1:4616325/150‑1 (MQ=255)
                      aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATgg                                                                                                                            <  1:597214/150‑1 (MQ=255)
                      aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATgg                                                                                                                            >  1:1828481/1‑150 (MQ=255)
                      aTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGACACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCTTgg                                                                                                                            >  2:1792485/1‑150 (MQ=255)
                         aTCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGaaa                                                                                                                         >  1:1191777/1‑150 (MQ=255)
                           cAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                       <  2:5094187/150‑1 (MQ=255)
                           cAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                       >  1:775881/1‑150 (MQ=255)
                             gCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGACATggg                                                                                                                     >  2:4941744/1‑150 (MQ=255)
                              cAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGt                                                                                                                    <  2:3192599/150‑1 (MQ=255)
                               aaGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTc                                                                                                                   <  1:2077625/150‑1 (MQ=255)
                                aGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCg                                                                                                                  >  2:4197493/1‑150 (MQ=255)
                                aGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCg                                                                                                                  >  2:2797478/1‑150 (MQ=255)
                                aGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCg                                                                                                                  <  2:5178529/150‑1 (MQ=255)
                                  tttttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCACCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTa                                                                                                                <  1:3959595/150‑1 (MQ=255)
                                   ttttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTAt                                                                                                               <  2:58175/150‑1 (MQ=255)
                                    tttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTAtt                                                                                                              <  1:1132182/150‑1 (MQ=255)
                                     ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTc                                                                                                             >  2:1325152/1‑150 (MQ=255)
                                     ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTc                                                                                                             <  2:1248399/150‑1 (MQ=255)
                                     ttCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTc                                                                                                             <  1:1410847/150‑1 (MQ=255)
                                      tCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCg                                                                                                            >  2:5278388/1‑150 (MQ=255)
                                       cgccgcCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCgg                                                                                                           <  1:3404837/150‑1 (MQ=255)
                                       cgccgcCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCgg                                                                                                           <  1:2475066/150‑1 (MQ=255)
                                        gccgccATGATGCTTTTCGGCCCCATCGTCACGTAGCGGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACACTATCAACCATCTCTCCCGGGATAGGGAACTCCACGGTTTGCTCATGGTAATGGGTCGTATTCGGc                                                                                                          <  2:2594421/150‑1 (MQ=255)
                                           gccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCaaa                                                                                                       >  1:1866950/1‑150 (MQ=255)
                                             cATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAAtt                                                                                                     <  1:4389999/150‑1 (MQ=255)
                                              atgatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTg                                                                                                    >  1:5089933/1‑150 (MQ=255)
                                              atgatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTg                                                                                                    <  2:879548/150‑1 (MQ=255)
                                              atgatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTg                                                                                                    >  1:1011889/1‑150 (MQ=255)
                                                gatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGcc                                                                                                  >  1:3038697/1‑150 (MQ=255)
                                                 atgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCca                                                                                                 <  2:1057709/150‑1 (MQ=255)
                                                  tgCTTTTCGGCCCCATCGTCACGTAGCTGTTAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCcac                                                                                                >  2:2978206/1‑150 (MQ=255)
                                                  tgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCcac                                                                                                >  1:4196359/1‑150 (MQ=255)
                                                   gCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCcaca                                                                                               <  1:2281992/150‑1 (MQ=255)
                                                        tCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAAt                                                                                          >  1:2729449/1‑150 (MQ=255)
                                                            ccccATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGAcc                                                                                      >  1:4137775/1‑150 (MQ=255)
                                                              ccATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCAt                                                                                    <  1:1331021/150‑1 (MQ=255)
                                                               cATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATc                                                                                   >  2:2394803/1‑150 (MQ=255)
                                                               cATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATc                                                                                   <  1:1623630/150‑1 (MQ=255)
                                                                aTCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCt                                                                                  <  2:782761/150‑1 (MQ=255)
                                                                aTCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCACTCCCTGGATAGGGACCTCCACGGGTTGCTCATGGACAGGGGTCGTATTCGGCCACTTGCCACAAAAATGACCATCt                                                                                  >  2:4448874/1‑150 (MQ=255)
                                                                 tCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTg                                                                                 <  2:1584890/150‑1 (MQ=255)
                                                                   gTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCa                                                                               >  1:3988960/1‑150 (MQ=255)
                                                                    tCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCAc                                                                              <  1:955897/150‑1 (MQ=255)
                                                                      aCGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCATCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCt                                                                            <  2:524566/150‑1 (MQ=255)
                                                                       cGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTa                                                                           >  2:2038057/1‑150 (MQ=255)
                                                                         tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAAt                                                                         >  1:3600225/1‑150 (MQ=255)
                                                                         tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAAt                                                                         >  1:2059798/1‑150 (MQ=255)
                                                                          aGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATc                                                                        <  1:4645757/150‑1 (MQ=255)
                                                                           gCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATcc                                                                       >  1:909707/1‑150 (MQ=255)
                                                                            cTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATccc                                                                      >  1:3211522/1‑150 (MQ=255)
                                                                             tGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCa                                                                     >  2:4513663/1‑150 (MQ=255)
                                                                              gTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAg                                                                    >  2:5098202/1‑150 (MQ=255)
                                                                                 aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcac                                                                 >  2:4399112/1‑150 (MQ=255)
                                                                                  ggCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcacc                                                                <  2:4065635/150‑1 (MQ=255)
                                                                                    cACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcaccac                                                              >  1:2950626/1‑150 (MQ=255)
                                                                                      cTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCAccacca                                                            <  2:3671552/150‑1 (MQ=255)
                                                                                        aaCGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCaaa                                                          >  2:302104/1‑150 (MQ=255)
                                                                                           gAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATcc                                                       >  1:876819/1‑150 (MQ=255)
                                                                                            aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCa                                                      >  2:101379/1‑150 (MQ=255)
                                                                                            aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCa                                                      >  1:2680879/1‑150 (MQ=255)
                                                                                               aGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACATAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                   <  2:2373293/150‑1 (MQ=255)
                                                                                               aGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                   >  1:2180815/1‑150 (MQ=255)
                                                                                                gTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTcc                                                  >  2:1145954/1‑150 (MQ=255)
                                                                                                    ccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCt                                              <  2:3135861/150‑1 (MQ=255)
                                                                                                     cgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCtt                                             <  2:931442/150‑1 (MQ=255)
                                                                                                        ccAACTCGCCATGCTACACAATATCAACCATCTCTCCCTGGTTAGGGAACTCCACTGTTTGCTCATGGAAATGGTTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAc                                          <  1:2354566/150‑1 (MQ=255)
                                                                                                        ccAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAc                                          >  1:251572/1‑150 (MQ=255)
                                                                                                         cAACTCGCCATTCGCCACAATATCAACCTTCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCTTATTCGGCAAATAGCCACAAAAATGAACATCTGCACCTAATCACAGCACCACCAAAACCAGTCCTCCTTCAcg                                         >  2:1607723/1‑150 (MQ=255)
                                                                                                         cAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                                         >  2:4975815/1‑150 (MQ=255)
                                                                                                         cAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                                         >  2:4976349/1‑150 (MQ=255)
                                                                                                             tCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGcc                                     <  2:4686339/150‑1 (MQ=255)
                                                                                                                ccATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt                                  <  2:19047/150‑1 (MQ=255)
                                                                                                                  aTGCGCCACAATCTCAACCCTCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGttt                                <  2:5233128/150‑1 (MQ=255)
                                                                                                                       ccACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGAt                           <  1:60994/150‑1 (MQ=255)
                                                                                                                          caATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCAt                        >  2:3957602/1‑150 (MQ=255)
                                                                                                                                       ctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGtt           >  1:3335698/1‑150 (MQ=255)
                                                                                                                                       ctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGtt           >  2:4738986/1‑150 (MQ=255)
                                                                                                                                       ctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGtt           >  2:2084297/1‑150 (MQ=255)
                                                                                                                                         ctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTAt         >  2:2479054/1‑150 (MQ=255)
                                                                                                                                          tcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTGATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATc        >  1:1901994/1‑150 (MQ=255)
                                                                                                                                          tcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATc        <  2:630345/150‑1 (MQ=255)
                                                                                                                                             cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAAt     <  1:2231509/150‑1 (MQ=255)
                                                                                                                                               ggATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAg   >  1:921805/1‑150 (MQ=255)
                                                                                                                                                gATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAGt  <  1:2880629/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAGT  >  CP009273/3892157‑3892450

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: