Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,275,040 C→A 19.0% intergenic (+56/‑295) agaA → / → agaS pseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,275,0400CA19.0% 657.4 / 123.2 262intergenic (+56/‑295)agaA/agaSpseudogene, N‑acetylgalactosamine‑6‑phosphate deacetylase fragment; putative N‑acetylgalactosamine‑6‑phosphate deacetylase/tagatose‑6‑phosphate ketose/aldose isomerase
Reads supporting (aligned to +/- strand):  ref base C (90/121);  new base A (20/29);  total (112/150)
Fisher's exact test for biased strand distribution p-value = 8.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

Only 100 of 282 total aligned reads displayed.

GCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTC  >  CP009273/3274891‑3275189
                                                                                                                                                     |                                                                                                                                                     
gCTTAAACCGGGCAAACGCGCCAGTGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTtttctttc                                                                                                                                                       <  1:2333281/150‑1 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTGTTTCGtt                                                                                                                                                    >  2:3714716/1‑150 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtt                                                                                                                                                    >  2:3204915/1‑150 (MQ=255)
    aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGttt                                                                                                                                                   >  2:1465056/1‑150 (MQ=255)
        cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                                                                                                                                               <  1:2689034/150‑1 (MQ=255)
        cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                                                                                                                                               <  1:3960441/150‑1 (MQ=255)
          ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCt                                                                                                                                             <  2:412079/150‑1 (MQ=255)
          ggCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCt                                                                                                                                             <  1:3989227/150‑1 (MQ=255)
             aaaCGCGCCAGCGTCGTTGTGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgtg                                                                                                                                          <  2:2517630/150‑1 (MQ=255)
             aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTATTTCGTTTTGTTAtgtg                                                                                                                                          >  2:1813448/1‑150 (MQ=255)
               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGGCTTTCTTTTTCTTTCGTTGTGTTATGTGCa                                                                                                                                        >  1:2622996/1‑150 (MQ=255)
               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCa                                                                                                                                        <  2:1278067/150‑1 (MQ=255)
                cgcgCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAg                                                                                                                                       <  1:4066826/150‑1 (MQ=255)
                cgcgCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAg                                                                                                                                       <  2:2477109/150‑1 (MQ=255)
                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGTTAGTTTGGTCCTTTATTGGGGGGTCACGTACCCCCGAAGGCCTTTAGTTTTCTTTGGTTTTGTTCTGTGCAGCg                                                                                                                                     >  1:783149/1‑150 (MQ=255)
                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCg                                                                                                                                     <  1:3592885/150‑1 (MQ=255)
                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCg                                                                                                                                     <  2:3798486/150‑1 (MQ=255)
                    ccAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGgt                                                                                                                                   <  2:1038193/150‑1 (MQ=255)
                     cAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCGTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtg                                                                                                                                  >  1:1908539/1‑150 (MQ=255)
                           cgtTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGAt                                                                                                                            >  2:1261995/1‑150 (MQ=255)
                               gcgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGa                                                                                                                        >  2:3548699/1‑150 (MQ=255)
                               gcgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCATTCTTTTTCTTTCGTTTTGATCTGTTCAGATGTGTCGGATgaga                                                                                                                        >  2:2863052/1‑150 (MQ=255)
                                cgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGAc                                                                                                                       <  2:2944532/150‑1 (MQ=255)
                                    ggATAGCGGGCTACATGTTCAACAAATCTGGATTCAGGGTCAATAAGAGTCGTTATGATAGTTTGATCCTTTATTGGGCCTTCACTTCACCCGAAAGGCCTATATTTTTCTTTAGTTTTGATCTGTTAAGCGGTGTCTGAAGCGaaccta                                                                                                                   >  2:1965583/1‑145 (MQ=255)
                                       tAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGGTAGTTTGGTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGGTTTGATCTGTGCAGCGGTGTCGGGTTCGACGCTAAcg                                                                                                                >  1:2577356/1‑150 (MQ=255)
                                       tAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcg                                                                                                                <  1:3518761/150‑1 (MQ=255)
                                        aGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcgc                                                                                                               <  2:930967/150‑1 (MQ=255)
                                         gCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAcgcg                                                                                                              >  1:1378197/1‑150 (MQ=255)
                                              cTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTa                                                                                                         <  2:3988618/150‑1 (MQ=255)
                                              cTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTa                                                                                                         <  2:2478682/150‑1 (MQ=255)
                                               tACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGGTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTATTTTTATTTCGTTTTGATCTGTGCAGCGGTGTGGGGTGGGACGCTAAAGCGGCTTAt                                                                                                        >  1:767106/1‑150 (MQ=255)
                                                  aTGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCg                                                                                                     >  2:1638738/1‑150 (MQ=255)
                                                  aTGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCg                                                                                                     >  2:3826860/1‑150 (MQ=255)
                                                          caaATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGGTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAACGCGTCTTATACGAccaaaag                                                                                             >  1:1314714/1‑145 (MQ=255)
                                                             aTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTCTTGGGCCTTCACTTCCCCCGTCAGGCCGTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGCTGCGACGCTAACGCCACTTATCCCACCAACCGTTg                                                                                          >  2:3356897/1‑150 (MQ=255)
                                                             aTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAACGCGTCTTATCCGACCAACACTTg                                                                                          >  1:1589982/1‑150 (MQ=255)
                                                             aTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTg                                                                                          <  2:3236579/150‑1 (MQ=255)
                                                                tGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTg                                                                                       >  1:3892725/1‑150 (MQ=255)
                                                                  gATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAc                                                                                     >  2:821000/1‑150 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    <  1:681575/150‑1 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTATTTTTCTTTCGGTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    >  2:653031/1‑150 (MQ=255)
                                                                    ttCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCGTTATTGGGGCTTCACGTCCCCCCGAAGGCCTTTATTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGGATCGCCCCAAACCCCTCTTAAACGACCAACACAAGGTGAccc                                                                                   >  1:437859/1‑149 (MQ=255)
                                                                            cAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCg                                                                           <  1:3730421/150‑1 (MQ=255)
                                                                            cAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCg                                                                           <  1:2034236/150‑1 (MQ=255)
                                                                             aaTTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCgg                                                                          <  2:4160580/150‑1 (MQ=255)
                                                                              aTTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGa                                                                         <  1:3065791/150‑1 (MQ=255)
                                                                              aTTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAACGCGTCTTATCCGAACAACAGTTGGTGAACGCAAGGCCGGa                                                                         >  1:2496814/1‑150 (MQ=255)
                                                                              aTTAGCTTCGTTTTGATAGTGTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGa                                                                         <  1:2717942/150‑1 (MQ=255)
                                                                                tAGCTTCGTTTTGATAGTTTGCTCCTTTATTTGGGCTTCACTTCCCCCGTAAGGCCTTTATTTTTCTTTCGTTTTGATCTGTGCATCGGTGTCGGATTCGACGCTAACGCGTCATATCCGACCAACAGTTTGTGAACGCAAAGCCGGATa                                                                       >  1:3215459/1‑150 (MQ=255)
                                                                                tAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATa                                                                       <  2:2769246/150‑1 (MQ=255)
                                                                                  gCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATaaa                                                                     >  2:3319007/1‑150 (MQ=255)
                                                                                   cTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAg                                                                    <  2:970792/150‑1 (MQ=255)
                                                                                    ttCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGc                                                                   >  2:1357677/1‑150 (MQ=255)
                                                                                     tCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCg                                                                  <  1:3279134/150‑1 (MQ=255)
                                                                                     tCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCg                                                                  <  1:1082895/150‑1 (MQ=255)
                                                                                     tCGTTTTGATAGTGTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCg                                                                  <  2:2291258/150‑1 (MQ=255)
                                                                                       gTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTTGTGACCGCAAGGCCGGATAAAGCGtt                                                                >  1:2795828/1‑150 (MQ=255)
                                                                                        ttttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGttt                                                               <  2:4215650/150‑1 (MQ=255)
                                                                                        ttttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGttt                                                               <  1:3690913/150‑1 (MQ=255)
                                                                                        ttttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGttt                                                               <  1:1507904/150‑1 (MQ=255)
                                                                                         tttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTg                                                              <  1:2241090/150‑1 (MQ=255)
                                                                                           tGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgcg                                                            >  2:3265836/1‑150 (MQ=255)
                                                                                            gATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgcgc                                                           >  2:1427277/1‑150 (MQ=255)
                                                                                            gATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgcgc                                                           >  1:2134138/1‑150 (MQ=255)
                                                                                                      tCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGc                                                 <  1:3356897/150‑1 (MQ=255)
                                                                                                      tCCTTTATTGGGCCTTCACGTCCCCCGTAAGGCCTTTCTTGTTCTGTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGc                                                 <  1:1159857/150‑1 (MQ=255)
                                                                                                       ccTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGAGACGCTAAACCCTCATATACCACCAAAACATTGTGAACCCAAAGCCCGATAAAGCGTTTTCGCCGCaaacacca                                                >  1:2953603/1‑143 (MQ=255)
                                                                                                         tttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              <  1:3204772/150‑1 (MQ=255)
                                                                                                         tttATTGGGCCGTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              >  2:3093044/1‑150 (MQ=255)
                                                                                                          ttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATc                                             <  2:2289677/150‑1 (MQ=255)
                                                                                                          ttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATc                                             <  2:3075945/150‑1 (MQ=255)
                                                                                                            aTTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATccc                                           >  2:1734355/1‑150 (MQ=255)
                                                                                                                 gCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTg                                      <  2:4260675/150‑1 (MQ=255)
                                                                                                                  ccTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGc                                     <  2:1106158/150‑1 (MQ=255)
                                                                                                                   cTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCt                                    <  1:2445139/150‑1 (MQ=255)
                                                                                                                    ttCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCtt                                   >  1:2384309/1‑150 (MQ=255)
                                                                                                                        cTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCtt                               <  1:3684101/150‑1 (MQ=255)
                                                                                                                        cTTCCCCCGTAAGGCCTTTCTTTTTCTGTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCtt                               <  1:871589/150‑1 (MQ=255)
                                                                                                                           cccccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAt                            <  1:1174782/150‑1 (MQ=255)
                                                                                                                            ccccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAAAGCGTCTTATCCGACCAAAAGTTGGTGACCCCAAGGCCGGATAAAACGTTTTCGCCGCAACCACCAAAACCATTTGCTTCCTTTAtc                           >  1:1464250/1‑150 (MQ=255)
                                                                                                                               cGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttt                        <  2:57561/150‑1 (MQ=255)
                                                                                                                                gTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                       <  1:707042/150‑1 (MQ=255)
                                                                                                                                gTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                       <  2:539930/150‑1 (MQ=255)
                                                                                                                                gTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctttt                       <  2:3029542/150‑1 (MQ=255)
                                                                                                                                 tAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttc                      <  2:1086339/150‑1 (MQ=255)
                                                                                                                                 tAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtattttc                      >  2:3583652/1‑150 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                     >  2:4148011/1‑150 (MQ=255)
                                                                                                                                    ggCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                   <  1:960693/150‑1 (MQ=255)
                                                                                                                                     gCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTtttctttc                  <  2:1030071/150‑1 (MQ=255)
                                                                                                                                      ccTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCa                 <  2:797659/150‑1 (MQ=255)
                                                                                                                                       ctttctttTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCaa                <  1:2356832/150‑1 (MQ=255)
                                                                                                                                       ctttctttTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCaa                <  1:180760/150‑1 (MQ=255)
                                                                                                                                           ctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt            >  1:2439640/1‑150 (MQ=255)
                                                                                                                                              tttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcac         >  2:905160/1‑150 (MQ=255)
                                                                                                                                              tttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCAACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcac         >  2:1838772/1‑150 (MQ=255)
                                                                                                                                                tctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaa       <  1:1664573/150‑1 (MQ=255)
                                                                                                                                                tctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaa       <  1:595031/150‑1 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  2:3591944/150‑1 (MQ=255)
                                                                                                                                                     cGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTc  <  2:2097561/150‑1 (MQ=255)
                                                                                                                                                     cGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTc  <  2:1593818/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTC  >  CP009273/3274891‑3275189

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