Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,275,040 C→A 7.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,0400CA 7.0% 939.8 / 31.7 264intergenic (+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 (97/146);  new base A (12/9);  total (109/155)
Fisher's exact test for biased strand distribution p-value = 1.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAATTTC  >  CP009273/3274891‑3275189
                                                                                                                                                     |                                                                                                                                                     
gCTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTCAGGCCTTTCTTtttctttc                                                                                                                                                       >  2:117431/1‑150 (MQ=255)
 cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCg                                                                                                                                                      >  2:4021947/1‑150 (MQ=255)
 cTTAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCg                                                                                                                                                      <  2:1292846/150‑1 (MQ=255)
  ttAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATATTTTGGGCCGTTATTGGGGCTTCACTTCCCCCGGAAGGGCGTTCTTTTTCTTTCgg                                                                                                                                                     >  1:3516750/1‑149 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTGCTTTATTGGGGGGTCACTTACCCCGTAAGGCCGTTGTGTTTGTTTAGtt                                                                                                                                                    >  1:3796512/1‑150 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtt                                                                                                                                                    <  2:2608132/150‑1 (MQ=255)
   tAAACCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGATACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtt                                                                                                                                                    <  2:2011916/150‑1 (MQ=255)
    aaaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTGTTTCGttt                                                                                                                                                   >  2:393665/1‑150 (MQ=255)
     aaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGtttt                                                                                                                                                  <  1:2777930/150‑1 (MQ=255)
     aaCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGtttt                                                                                                                                                  <  2:2813157/150‑1 (MQ=255)
      aCCGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTg                                                                                                                                                 <  2:912093/150‑1 (MQ=255)
       ccGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGa                                                                                                                                                <  2:682472/150‑1 (MQ=255)
        cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                                                                                                                                               <  2:3423295/150‑1 (MQ=255)
        cGGGCAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGAt                                                                                                                                               <  1:368944/150‑1 (MQ=255)
            cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgt                                                                                                                                           <  2:2559431/150‑1 (MQ=255)
            cAAACGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgt                                                                                                                                           >  2:4403989/1‑150 (MQ=255)
             aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgtg                                                                                                                                          <  1:3897371/150‑1 (MQ=255)
             aaaCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCtgtg                                                                                                                                          >  1:284513/1‑150 (MQ=255)
               aCGCGCCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCa                                                                                                                                        <  2:2244972/150‑1 (MQ=255)
                 gcgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGc                                                                                                                                      >  2:368719/1‑150 (MQ=255)
                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCg                                                                                                                                     >  2:4469457/1‑150 (MQ=255)
                  cgcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCGTCGTTTTGTTAGTTTTCTACGTTATTGGGGGTTCACATACCCCGTCCGGGCTTTATTTTTGTTTAGTTTTGTTCTTTTCAGCg                                                                                                                                     >  1:3033324/1‑150 (MQ=255)
                   gcCAGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATGAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCgg                                                                                                                                    <  1:1346153/150‑1 (MQ=255)
                      aGCGTCGTTGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGgtgt                                                                                                                                 <  1:421284/150‑1 (MQ=255)
                             ttGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTCCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGc                                                                                                                          <  2:1285372/150‑1 (MQ=255)
                              tGCGCTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCg                                                                                                                         <  1:3970208/150‑1 (MQ=255)
                                cgcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGGTCACTTCCCCCGTACGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGAc                                                                                                                       >  1:1209774/1‑150 (MQ=255)
                                 gcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTACGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACg                                                                                                                      >  2:1526219/1‑150 (MQ=255)
                                 gcTGGATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACg                                                                                                                      >  2:1018329/1‑150 (MQ=255)
                                    ggATAGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGGGTCACTTCCCCCGTCAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGGGGGGTGCGACGCTa                                                                                                                   >  1:4343436/1‑150 (MQ=255)
                                        aGCGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTGTTGATAGTTTGCGCCGTTATTGGGCCTTCACTTCCCCCGTAAGGCGTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGTTGCGACGCtcaagc                                                                                                               >  1:3201302/1‑145 (MQ=255)
                                          cGGGCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGt                                                                                                             <  2:462719/150‑1 (MQ=255)
                                            ggCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCAAAAGCGTCt                                                                                                           >  1:2834183/1‑150 (MQ=255)
                                            ggCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCt                                                                                                           <  2:2072192/150‑1 (MQ=255)
                                            ggCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCt                                                                                                           >  1:4120550/1‑150 (MQ=255)
                                             gCTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCtt                                                                                                          <  2:3535003/150‑1 (MQ=255)
                                              cTACATGTGCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATATGTGCAGCGGTGTCGGATTCGAAGCTAAACCGTCGTa                                                                                                         >  1:4269409/1‑150 (MQ=255)
                                                     tgCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGAcc                                                                                                  <  2:1271605/150‑1 (MQ=255)
                                                     tgCAACAAATCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGAcc                                                                                                  >  2:2902890/1‑150 (MQ=255)
                                                     tgCAACAAATCTGGATTCAGGGTCAATGAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGAcc                                                                                                  <  2:3212795/150‑1 (MQ=255)
                                                           aaaTCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGt                                                                                            <  1:3926889/150‑1 (MQ=255)
                                                              tCTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTgg                                                                                         <  1:486766/150‑1 (MQ=255)
                                                               cTGGATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGt                                                                                        <  2:4443073/150‑1 (MQ=255)
                                                                 ggATTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGa                                                                                      <  1:7651/150‑1 (MQ=255)
                                                                   aTTCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGAcc                                                                                    >  2:1313835/1‑150 (MQ=255)
                                                                     tCAGGGTCAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGAAGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGc                                                                                  >  2:2163913/1‑150 (MQ=255)
                                                                        gggTCAATTAGCTTCGTTTTGATAGTTTTCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCaat                                                                               >  2:777916/1‑149 (MQ=255)
                                                                            cAATTAGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACATCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACCCAAAACCCGCTTATAACACCAACACATGGTGACCGCAAGGccc                                                                           >  1:1650065/1‑149 (MQ=255)
                                                                               ttAGCGTCGTTTTGATAGTTTGCTCCTTTAGTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGAt                                                                        <  2:838443/150‑1 (MQ=255)
                                                                               tgaGCTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTGTTCATTCGTATAGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGAt                                                                        <  2:2471234/148‑1 (MQ=255)
                                                                                   cTTCGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTTTGCAGCGGTGTCGGAAGCGACGCAAACGCGTCATATACGACCAAAACTTGGTTAACGCAAAGCCGGATAAAg                                                                    >  1:2445608/1‑150 (MQ=255)
                                                                                    ttCGTTTTGATAGTTGGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGc                                                                   <  2:4066949/150‑1 (MQ=255)
                                                                                      cGTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGt                                                                 >  2:1394540/1‑150 (MQ=255)
                                                                                       gTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGtt                                                                >  1:2178697/1‑150 (MQ=255)
                                                                                       gTTTTGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGtt                                                                <  2:1745852/150‑1 (MQ=255)
                                                                                          ttGATAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTgc                                                             <  2:4306590/150‑1 (MQ=255)
                                                                                           tGATAGTTTGCTCCTTTATTGGGGCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATACGACCTACAGTTGGTTTCCGCAAGGCCGGATAAAACGTTTgcg                                                            >  1:885644/1‑150 (MQ=255)
                                                                                             aTAGTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGcgcc                                                          <  2:2552720/150‑1 (MQ=255)
                                                                                              tAGTTTGCTCCTTTATTGCGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGcgccg                                                         <  1:3697383/150‑1 (MQ=255)
                                                                                                gTTTGCTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCa                                                       <  1:303780/150‑1 (MQ=255)
                                                                                                     cTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAg                                                  <  1:278990/150‑1 (MQ=255)
                                                                                                     cTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAg                                                  <  1:1572286/150‑1 (MQ=255)
                                                                                                     cTCCTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAg                                                  <  1:1456029/150‑1 (MQ=255)
                                                                                                        cTTTATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCaa                                               <  1:2739849/150‑1 (MQ=255)
                                                                                                         tttATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAAt                                              <  2:2253951/150‑1 (MQ=255)
                                                                                                           tATTGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGCTAAAGCGTTTGCGCCGCATCCAGCAATcc                                            >  2:2668575/1‑150 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          <  2:4181923/150‑1 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          <  2:3855775/150‑1 (MQ=255)
                                                                                                             ttGGGCCTTCACTTCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCt                                          <  1:595408/150‑1 (MQ=255)
                                                                                                                 gCCTTCACTTCCCCCGTAAGGACTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTg                                      <  1:292438/150‑1 (MQ=255)
                                                                                                                         ttCCCCCGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCttt                              <  1:1347238/150‑1 (MQ=255)
                                                                                                                           cccccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAt                            <  1:1598457/150‑1 (MQ=255)
                                                                                                                            ccacGTTAGGCCTTTCTTTTTCTTTCGATTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAAATCGTCTTAGCCGACCTTCTGTTGGTGACCGCTCGGGCGTATAAAGCGTTTGCGCTGCATCCAGCAATCCCTTTTGCTTCCTTTAtc                           <  1:2097440/147‑1 (MQ=255)
                                                                                                                             cccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtct                          <  2:2941121/150‑1 (MQ=255)
                                                                                                                              ccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtctt                         <  1:1139282/150‑1 (MQ=255)
                                                                                                                              ccGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATACGACCAACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCACTCCCATTTGCTTCCTGTAtctt                         >  1:4400858/1‑150 (MQ=255)
                                                                                                                               cGTAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTAAAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTTCGCCGCATCCAGCAATCCCATTTGCTTCCTTTAtcttt                        >  1:1409193/1‑150 (MQ=255)
                                                                                                                                 tAAGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttc                      <  1:307975/150‑1 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                     <  2:4366331/150‑1 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                     >  2:611584/1‑150 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttct                     <  2:3766822/150‑1 (MQ=255)
                                                                                                                                  aaGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACCCAAACCCGTCTTATACCACCAAAACTTGGTGACCGCCACGCCGGATAAAGCGTTTTCGCCCCAACCACCAAACCCAATTGCTTCCGTTAtcttttat                     >  1:928281/1‑148 (MQ=255)
                                                                                                                                   aGGCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttctt                    >  2:2179485/1‑150 (MQ=255)
                                                                                                                                    ggCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATTCGACGCAAACGCCTCTTATACGACCTACAGTTGGTGACCGCAAGGCCGGATAAATCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                   >  1:1751405/1‑150 (MQ=255)
                                                                                                                                    ggCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGAGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTAtcttttcttt                   <  1:1505754/150‑1 (MQ=255)
                                                                                                                                     gCCTTTCTTTTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTtttctttc                  <  2:926829/150‑1 (MQ=255)
                                                                                                                                       ctttctttTTCTTTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCGAACGCGTATTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTttgctttccc                <  1:659421/150‑3 (MQ=255)
                                                                                                                                        tttctttttcttTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAAc               >  2:3244145/1‑150 (MQ=255)
                                                                                                                                        tttctttttcttTCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAAc               <  1:3884971/150‑1 (MQ=255)
                                                                                                                                        tttctttttcttTAGTTTTGATCTGTGCAGCGGTGTCGGATGCGCCGCTAAAGCGTCTTATACGACCTACAGTTGGTTACCGCAAGGCCGGCTAAATCGTTTTCGCCGCAACCAGCAAACCCTTTTGCTTTCTTTATCTTTTATTTCaaa               >  1:3140549/1‑149 (MQ=255)
                                                                                                                                           ctttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGAt            <  1:2318297/150‑1 (MQ=255)
                                                                                                                                             ttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATca          <  2:2427096/150‑1 (MQ=255)
                                                                                                                                             ttttctttCGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATca          <  2:1860296/150‑1 (MQ=255)
                                                                                                                                               ttctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATcaca        <  1:1929389/150‑1 (MQ=255)
                                                                                                                                                tctttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACaa       <  1:207549/150‑1 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  2:277447/150‑1 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  1:1067851/150‑1 (MQ=255)
                                                                                                                                                  tttcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  1:1018329/150‑1 (MQ=255)
                                                                                                                                                  ttgcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAt     <  2:4326343/147‑1 (MQ=255)
                                                                                                                                                    tcGTTTTGATCTGTGCAGCGGTGTCGGATGCGACGCTAACGCGTCTTATCCGACCTACAGTTGGTGACCGCAAGGCCGGATAAAGCGTTTGCGCCGCATCCAGCAATCCCTTTTGCTTCCTTTATCTTTTCTTTCAACGATCACAAAttt   <  2:591663/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: