Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,274,047 T→A 5.2% *158R (TGA→AGA)  agaV → N‑acetylgalactosamine‑specific enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,274,0470TA 5.2% 1067.3 / 31.6 315*158R (TGA→AGA) agaVN‑acetylgalactosamine‑specific enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base T (172/125);  new base A (9/8);  total (182/133)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 331 total aligned reads displayed.

ACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATCT  >  CP009273/3273898‑3274196
                                                                                                                                                     |                                                                                                                                                     
aCTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACttt                                                                                                                                                       >  2:2296082/1‑150 (MQ=255)
aCTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACttt                                                                                                                                                       <  2:3169430/150‑1 (MQ=255)
aCTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACttt                                                                                                                                                       >  1:899751/1‑150 (MQ=255)
  tACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGa                                                                                                                                                     >  1:3239002/1‑150 (MQ=255)
    cGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAgg                                                                                                                                                   <  2:2967501/150‑1 (MQ=255)
       cAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGAt                                                                                                                                                <  2:3844611/150‑1 (MQ=255)
         aTGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGTAAACCGCTGGGGTGGTATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATtc                                                                                                                                              <  1:3394340/150‑1 (MQ=255)
         aTGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATtc                                                                                                                                              <  1:5122108/150‑1 (MQ=255)
           ggCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATtcat                                                                                                                                            >  1:424302/1‑150 (MQ=255)
             cAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATtcatca                                                                                                                                          >  1:2833341/1‑150 (MQ=255)
             cAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATtcatca                                                                                                                                          >  1:1470934/1‑150 (MQ=255)
               aacaacAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATg                                                                                                                                        >  1:444582/1‑150 (MQ=255)
               aacaacAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATg                                                                                                                                        <  1:850788/150‑1 (MQ=255)
                     aaaTCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATc                                                                                                                                  >  2:1747202/1‑150 (MQ=255)
                      aaTCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCa                                                                                                                                 >  2:3786723/1‑150 (MQ=255)
                       aTCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAg                                                                                                                                <  2:1421226/150‑1 (MQ=255)
                          gCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCt                                                                                                                             >  2:483455/1‑150 (MQ=255)
                             aaaaCGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGtt                                                                                                                          >  2:1030126/1‑150 (MQ=255)
                               aaCGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGc                                                                                                                        >  2:1760265/1‑150 (MQ=255)
                               aaCGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTGAGAGGGATTCATCTTGGAAATCAGCCTGTTGc                                                                                                                        <  2:4948630/150‑1 (MQ=255)
                                  ggTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAgg                                                                                                                     <  2:3590571/150‑1 (MQ=255)
                                       cTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCAttt                                                                                                                <  2:3087409/150‑1 (MQ=255)
                                         gtgGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGc                                                                                                              >  2:4504370/1‑150 (MQ=255)
                                          tgGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCg                                                                                                             <  1:2930160/150‑1 (MQ=255)
                                          tgGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTAGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCg                                                                                                             <  1:1314057/150‑1 (MQ=255)
                                            gATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGtt                                                                                                           <  2:714617/150‑1 (MQ=255)
                                            gATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGtt                                                                                                           >  2:4054050/1‑150 (MQ=255)
                                             aTGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTg                                                                                                          >  1:3733097/1‑150 (MQ=255)
                                              tGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTgg                                                                                                         <  1:1268820/150‑1 (MQ=255)
                                              tGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTgg                                                                                                         >  1:4054888/1‑150 (MQ=255)
                                                  ggCGATATCGCAGCATTTATCGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCAt                                                                                                     <  2:3210259/150‑1 (MQ=255)
                                                  ggCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCAt                                                                                                     <  1:725288/150‑1 (MQ=255)
                                                    cGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTa                                                                                                   >  2:2820309/1‑150 (MQ=255)
                                                    cGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTa                                                                                                   <  2:2097486/150‑1 (MQ=255)
                                                      atatCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATc                                                                                                 <  1:3932186/150‑1 (MQ=255)
                                                       tatCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCg                                                                                                >  2:4102363/1‑150 (MQ=255)
                                                         tCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGcc                                                                                              >  1:2311720/1‑150 (MQ=255)
                                                            cagcaTTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCttt                                                                                           <  2:1374910/150‑1 (MQ=255)
                                                             agcaTTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTa                                                                                          <  1:3453629/150‑1 (MQ=255)
                                                             agcaTTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTa                                                                                          >  2:1941382/1‑150 (MQ=255)
                                                              gcaTTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTAt                                                                                         <  2:3469439/150‑1 (MQ=255)
                                                                   tAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTg                                                                                    <  1:2279118/150‑1 (MQ=255)
                                                                       gACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCt                                                                                <  2:573363/150‑1 (MQ=255)
                                                                              aaaCCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATg                                                                         >  1:2328961/1‑150 (MQ=255)
                                                                              aaaCCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATg                                                                         <  2:4983003/150‑1 (MQ=255)
                                                                               aaCCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGt                                                                        <  2:1713113/150‑1 (MQ=255)
                                                                               aaCCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGt                                                                        >  2:4629886/1‑150 (MQ=255)
                                                                                 ccGCTGGGGTGGAATGCGTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGttt                                                                      <  2:4067783/150‑1 (MQ=255)
                                                                                   gCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTaa                                                                    >  2:1575127/1‑150 (MQ=255)
                                                                                     tGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACg                                                                  <  2:2461728/150‑1 (MQ=255)
                                                                                     tGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACg                                                                  >  1:3153486/1‑150 (MQ=255)
                                                                                      ggggTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACgg                                                                 >  1:1372065/1‑150 (MQ=255)
                                                                                         gTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCt                                                              <  1:4547756/150‑1 (MQ=255)
                                                                                          tGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTAGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTa                                                             >  2:2650693/1‑150 (MQ=255)
                                                                                           ggAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTaa                                                            <  1:3812098/150‑1 (MQ=255)
                                                                                            gAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAAc                                                           >  2:4364898/1‑150 (MQ=255)
                                                                                             aaTGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAAcc                                                          >  1:2753895/1‑150 (MQ=255)
                                                                                              aTGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAAccc                                                         >  2:4390960/1‑150 (MQ=255)
                                                                                              aTGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAAccc                                                         <  2:922693/150‑1 (MQ=255)
                                                                                               tGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCa                                                        >  2:1905031/1‑150 (MQ=255)
                                                                                               tGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCa                                                        >  2:4458961/1‑150 (MQ=255)
                                                                                               tGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTAGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCa                                                        >  1:1619417/1‑150 (MQ=255)
                                                                                                gCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCat                                                       <  2:3679653/150‑1 (MQ=255)
                                                                                                   tCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTAGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATg                                                    <  1:1020744/150‑1 (MQ=255)
                                                                                                     gTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATCTGTTTAACGGCCTAACCCCTATGCa                                                  >  2:2375836/1‑150 (MQ=255)
                                                                                                     gTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCa                                                  >  2:567032/1‑150 (MQ=255)
                                                                                                     gTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCAGTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCa                                                  >  2:5111273/1‑150 (MQ=255)
                                                                                                      ttCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCAc                                                 >  2:867277/1‑150 (MQ=255)
                                                                                                       tCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCTCCTTTATCGCTGGCATGGATATGTTTAACGGCCTAACCAATATGCAcc                                                >  2:593712/1‑150 (MQ=255)
                                                                                                           ggCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGccc                                            <  1:1641680/150‑1 (MQ=255)
                                                                                                                  cGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTc                                     <  1:3656897/150‑1 (MQ=255)
                                                                                                                     cAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTc                                  >  2:237792/1‑150 (MQ=255)
                                                                                                                     cAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTc                                  <  2:2942123/150‑1 (MQ=255)
                                                                                                                      agagCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCg                                 >  1:4720339/1‑150 (MQ=255)
                                                                                                                        agCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGc                               >  1:2646987/1‑150 (MQ=255)
                                                                                                                          ccTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGccc                             <  1:2172015/150‑1 (MQ=255)
                                                                                                                           ctgctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCg                            >  2:434541/1‑150 (MQ=255)
                                                                                                                            tgctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGt                           <  1:1044331/150‑1 (MQ=255)
                                                                                                                            tgctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGt                           >  2:1150714/1‑150 (MQ=255)
                                                                                                                             gctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGtt                          <  1:532628/150‑1 (MQ=255)
                                                                                                                              ctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTg                         >  2:1298771/1‑150 (MQ=255)
                                                                                                                              ctgTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTg                         >  2:4769756/1‑150 (MQ=255)
                                                                                                                               tgtgGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATAGGCACCGCCCGGTGGTCCTCGGCCCGTTgg                        <  2:2600317/150‑1 (MQ=255)
                                                                                                                                 tgGACCTCTTTAAATTACTTTGAGGGATTCATCACGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTc                      <  2:2554729/150‑1 (MQ=255)
                                                                                                                                  ggACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCg                     >  1:4183601/1‑150 (MQ=255)
                                                                                                                                  ggACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCAAATGCACCGCCCGGAGGTCCTCGGCCCGTTGGTCg                     >  2:4229868/1‑150 (MQ=255)
                                                                                                                                    aCCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCggg                   >  1:3910814/1‑150 (MQ=255)
                                                                                                                                    aCCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCggg                   <  2:5126896/150‑1 (MQ=255)
                                                                                                                                     ccTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGc                  >  1:432450/1‑150 (MQ=255)
                                                                                                                                       tctTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGAGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTg                >  1:3605879/1‑150 (MQ=255)
                                                                                                                                         tttAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGt              <  1:4899607/150‑1 (MQ=255)
                                                                                                                                            aaaTTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACt           >  2:3564815/1‑150 (MQ=255)
                                                                                                                                               ttACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTgg        >  1:1434617/1‑150 (MQ=255)
                                                                                                                                                tACTTAGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGc       <  1:3102362/150‑1 (MQ=255)
                                                                                                                                                  cTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGa     >  1:3511022/1‑150 (MQ=255)
                                                                                                                                                    ttGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATc   <  1:3525203/150‑1 (MQ=255)
                                                                                                                                                     tGAGGGATTCTTCATGGATATCAGCCTGTTGCAGGCTTTTGCGTTGGGCATTTTCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATCt  <  1:784470/150‑1 (MQ=255)
                                                                                                                                                     tGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATCt  <  1:2121711/150‑1 (MQ=255)
                                                                                                                                                     tGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATCt  <  1:2407588/150‑1 (MQ=255)
                                                                                                                                                     tGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGGTGGTCGGGCTGGTACTTGGCGATCt  >  2:3429976/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTACGCCAATGGCAAACAACAAATCGCCAAAACGGTTTCTGTGGATGCGGGCGATATCGCAGCATTTAACGACCTGAAAACCGCTGGGGTGGAATGCTTCGTTCAGGGCGTCCCGACAGAGCCTGCTGTGGACCTCTTTAAATTACTTTGAGGGATTCATCATGGAAATCAGCCTGTTGCAGGCATTTGCGTTGGGCATTATCGCCTTTATCGCTGGCCTGGATATGTTTAACGGCCTAACCCATATGCACCGCCCGGTGGTCCTCGGCCCGTTGGTCGGGCTGGTACTTGGCGATCT  >  CP009273/3273898‑3274196

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: