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,3040CA96.2% 796.6 / 19.7 260P150P (CCG→CCT) yieKputative 6‑phosphogluconolactonase
Reads supporting (aligned to +/- strand):  ref base C (3/7);  new base A (133/117);  total (136/124)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 260 total aligned reads displayed.

GCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAG  >  CP009273/3892158‑3892449
                                                                                                                                                  |                                                                                                                                                 
gcTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAgg                                                                                                                                                 >  1:281355/1‑149 (MQ=255)
  tttACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGa                                                                                                                                               >  1:1217643/1‑149 (MQ=255)
    tACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAAc                                                                                                                                             >  2:3614408/1‑149 (MQ=255)
    tACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCACCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAAc                                                                                                                                             >  2:2039252/1‑149 (MQ=255)
      cccGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTc                                                                                                                                           <  2:2997236/149‑1 (MQ=255)
      cccGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTc                                                                                                                                           <  2:1681117/149‑1 (MQ=255)
      cccGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTc                                                                                                                                           <  2:3709108/149‑1 (MQ=255)
        cgcgCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCa                                                                                                                                         <  2:2466211/149‑1 (MQ=255)
         gcgcCGCTGACTATTATCAGCTAGTTTTTCGCCGCCATGTTGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCAc                                                                                                                                        <  1:2730803/149‑1 (MQ=255)
          cgccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACg                                                                                                                                       <  1:3571087/149‑1 (MQ=255)
           gccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACgg                                                                                                                                      >  2:743062/1‑149 (MQ=255)
           gccgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACgg                                                                                                                                      <  2:2285207/149‑1 (MQ=255)
             cgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGtt                                                                                                                                    <  1:1070309/149‑1 (MQ=255)
             cgcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGtt                                                                                                                                    <  1:42977/149‑1 (MQ=255)
              gcTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGttt                                                                                                                                   <  1:195015/149‑1 (MQ=255)
                 gACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                >  2:2584245/1‑149 (MQ=255)
                 gACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCt                                                                                                                                >  1:1382694/1‑149 (MQ=255)
                  aCAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTc                                                                                                                               >  1:2982650/1‑149 (MQ=255)
                    aaTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCAt                                                                                                                             <  1:705064/149‑1 (MQ=255)
                    aaTTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCAt                                                                                                                             >  1:76801/1‑149 (MQ=255)
                      ttATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATgg                                                                                                                           >  2:3974941/1‑149 (MQ=255)
                       tATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGa                                                                                                                          >  1:3673116/1‑149 (MQ=255)
                         tCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGaaa                                                                                                                        <  2:1455131/149‑1 (MQ=255)
                           aGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATg                                                                                                                      <  1:3099861/149‑1 (MQ=255)
                             cAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATggg                                                                                                                    >  1:2925272/1‑149 (MQ=255)
                                        ccgccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTc                                                                                                            <  1:862819/146‑1 (MQ=255)
                                           ccATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCaaa                                                                                                      >  2:1002225/1‑149 (MQ=255)
                                            cATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAAt                                                                                                     <  1:2242992/149‑1 (MQ=255)
                                             atgatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAAtt                                                                                                    >  2:1294722/1‑149 (MQ=255)
                                               gatgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGc                                                                                                  >  2:136806/1‑149 (MQ=255)
                                                 tgCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCca                                                                                                <  1:1054681/149‑1 (MQ=255)
                                                       tCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACaaaaa                                                                                          <  2:826321/149‑1 (MQ=255)
                                                       tCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACaaaaa                                                                                          >  1:3614927/1‑149 (MQ=255)
                                                       tCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACaaaaa                                                                                          <  2:1092512/149‑1 (MQ=255)
                                                         ggCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATg                                                                                        >  2:862046/1‑149 (MQ=255)
                                                          gCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGa                                                                                       <  1:67512/149‑1 (MQ=255)
                                                              cATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCAt                                                                                   >  2:2891180/1‑149 (MQ=255)
                                                              cATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCAt                                                                                   >  1:1270954/1‑149 (MQ=255)
                                                               aTCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATc                                                                                  <  2:1827684/149‑1 (MQ=255)
                                                                        tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTaa                                                                         >  1:3711955/1‑149 (MQ=255)
                                                                        tAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTaa                                                                         >  1:2789468/1‑149 (MQ=255)
                                                                         aGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAAt                                                                        >  2:2375033/1‑149 (MQ=255)
                                                                         aGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAAt                                                                        >  1:725772/1‑149 (MQ=255)
                                                                           cTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATcc                                                                      <  2:3591374/149‑1 (MQ=255)
                                                                           cTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATcc                                                                      >  1:2167755/1‑149 (MQ=255)
                                                                            tGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATccc                                                                     <  2:633178/149‑1 (MQ=255)
                                                                            tGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATccc                                                                     <  2:2696250/149‑1 (MQ=255)
                                                                              tCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAg                                                                   >  2:2720474/1‑149 (MQ=255)
                                                                                aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGca                                                                 >  2:102875/1‑149 (MQ=255)
                                                                                aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGca                                                                 >  1:3582980/1‑149 (MQ=255)
                                                                                aGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGca                                                                 >  1:3650973/1‑149 (MQ=255)
                                                                                 ggCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcac                                                                >  2:989697/1‑149 (MQ=255)
                                                                                  gCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcacc                                                               >  2:415414/1‑149 (MQ=255)
                                                                                   cACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcacca                                                              >  1:3515531/1‑149 (MQ=255)
                                                                                   cACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATAACAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcacca                                                              >  2:3909180/1‑149 (MQ=255)
                                                                                    aCTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGcaccac                                                             <  1:3682691/149‑1 (MQ=255)
                                                                                       aaCGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCaa                                                          >  1:643049/1‑149 (MQ=255)
                                                                                          gAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATc                                                       <  2:254912/149‑1 (MQ=255)
                                                                                          gAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATc                                                       <  2:2218393/149‑1 (MQ=255)
                                                                                          gAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCGGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATc                                                       <  2:2330110/149‑1 (MQ=255)
                                                                                           aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATcc                                                      <  2:2166647/149‑1 (MQ=255)
                                                                                           aaaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATcc                                                      >  2:2353415/1‑149 (MQ=255)
                                                                                            aaaGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCa                                                     <  1:1269428/149‑1 (MQ=255)
                                                                                               gTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                  >  1:1817643/1‑149 (MQ=255)
                                                                                               gTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTc                                                  <  1:3142819/149‑1 (MQ=255)
                                                                                                   ccgccCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGcc                                              <  1:2543505/149‑1 (MQ=255)
                                                                                                         aaCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcg                                        >  1:2851041/1‑149 (MQ=255)
                                                                                                          aCTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCAcgc                                       >  2:3726310/1‑149 (MQ=255)
                                                                                                            tCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACgcgc                                     >  2:2506617/1‑149 (MQ=255)
                                                                                                            tCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACgcgc                                     <  1:1453469/149‑1 (MQ=255)
                                                                                                               ccATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAg                                  >  2:97896/1‑149 (MQ=255)
                                                                                                                cATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt                                 >  2:935567/1‑149 (MQ=255)
                                                                                                                cATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt                                 <  2:2038819/149‑1 (MQ=255)
                                                                                                                cATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt                                 <  2:4073807/149‑1 (MQ=255)
                                                                                                                cATGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGt                                 >  1:855589/1‑149 (MQ=255)
                                                                                                                  tGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGttt                               >  1:1620803/1‑149 (MQ=255)
                                                                                                                  tGCGCCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGttt                               >  2:333660/1‑149 (MQ=255)
                                                                                                                   gcgcCACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTc                              <  2:1598768/149‑1 (MQ=255)
                                                                                                                      ccACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGa                           >  2:276063/1‑149 (MQ=255)
                                                                                                                      ccACAATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGa                           <  1:4025155/149‑1 (MQ=255)
                                                                                                                       cacaATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGAt                          <  2:2981932/149‑1 (MQ=255)
                                                                                                                         caATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa                        <  2:281355/149‑1 (MQ=255)
                                                                                                                         caATATCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCa                        <  1:910167/149‑1 (MQ=255)
                                                                                                                            tatCAACCATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGc                     <  2:308285/149‑1 (MQ=255)
                                                                                                                             atCACCCATCTCTCCCTGGATAGGGAACTCCACGGTGTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCt                    <  2:3118069/149‑1 (MQ=255)
                                                                                                                                   cATCTCTCCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGt              >  2:127480/1‑149 (MQ=255)
                                                                                                                                     tctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAg            <  1:2584245/149‑1 (MQ=255)
                                                                                                                                     tctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAg            <  1:515221/149‑1 (MQ=255)
                                                                                                                                      ctctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGt           <  1:1290627/149‑1 (MQ=255)
                                                                                                                                        ctcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTa         >  1:1005421/1‑149 (MQ=255)
                                                                                                                                         tcCCTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTAt        <  1:831410/149‑1 (MQ=255)
                                                                                                                                          cccTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATc       <  1:3974956/149‑1 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     <  2:548590/149‑1 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     <  1:3928954/149‑1 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     <  2:3891631/149‑1 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     >  1:1644914/1‑149 (MQ=255)
                                                                                                                                            cTGGATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCaa     <  2:355144/149‑1 (MQ=255)
                                                                                                                                              ggATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATa   >  1:1923562/1‑149 (MQ=255)
                                                                                                                                               gATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAg  >  2:1778211/1‑149 (MQ=255)
                                                                                                                                               gATAGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAg  >  1:1255840/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCTTTACCCGCGCCGCTGACAATTATCAGCAAGTTTTTCGCCGCCATGATGCTTTTCGGCCCCATCGTCACGTAGCTGTCAGGCACTAACGAAAAGTCGCCGCCCAACTCGCCATGCGCCACAATATCAACCATCTCTCCCTGGATCGGGAACTCCACGGTTTGCTCATGGAAATGGGTCGTATTCGGCAAATTGCCACAAAAATGACCATCTGCACCTAATCCCAGCACCACCAAATCCAGTCCGCCTTCACGCGCCAGTTTCTGATCATGCTCGCGGTAGTTATCAATAG  >  CP009273/3892158‑3892449

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