Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,575,339 G→T 100% F481L (TTC→TTA ydeN ← putative Ser‑type periplasmic non‑aryl sulfatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,575,3390GT95.4% 726.4 / 23.6 233F481L (TTC→TTAydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (6/5);  new base T (100/122);  total (106/127)
Fisher's exact test for biased strand distribution p-value = 5.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 233 total aligned reads displayed.

TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGT  >  CP009273/1575192‑1575483
                                                                                                                                                   |                                                                                                                                                
 aCCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAt                                                                                                                                               <  2:2452026/150‑1 (MQ=255)
  ccACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAtt                                                                                                                                              >  2:4486715/1‑150 (MQ=255)
  ccACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAtt                                                                                                                                              >  1:3189898/1‑150 (MQ=255)
  ccACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAtt                                                                                                                                              >  2:2153874/1‑150 (MQ=255)
   cACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTg                                                                                                                                             >  2:4572220/1‑150 (MQ=255)
    aCGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTgg                                                                                                                                            <  1:2825483/150‑1 (MQ=255)
     cGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATCAGATAATTGGc                                                                                                                                           >  2:741359/1‑150 (MQ=255)
       ccTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCtt                                                                                                                                         >  2:876006/1‑150 (MQ=255)
       ccTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCtt                                                                                                                                         <  1:3046841/150‑1 (MQ=255)
          tGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAg                                                                                                                                      >  1:2703416/1‑150 (MQ=255)
           gCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGt                                                                                                                                     <  1:3373868/150‑1 (MQ=255)
           gCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGt                                                                                                                                     <  1:2389597/150‑1 (MQ=255)
            cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                    >  1:1213192/1‑150 (MQ=255)
              tctcTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCt                                                                                                                                  <  1:4393958/150‑1 (MQ=255)
                 ctTTAACGACCTGCGGATTGGTGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAg                                                                                                                               <  1:4268275/150‑1 (MQ=255)
                 ctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAg                                                                                                                               >  1:3763082/1‑150 (MQ=255)
                  tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                              <  1:389373/150‑1 (MQ=255)
                  tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                              >  2:3143342/1‑150 (MQ=255)
                   ttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtg                                                                                                                             <  2:4294117/150‑1 (MQ=255)
                   ttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAgtg                                                                                                                             >  2:1745066/1‑150 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                            >  2:2667394/1‑150 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgg                                                                                                                            >  2:593687/1‑149 (MQ=255)
                           cTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGtt                                                                                                                     >  1:1174555/1‑150 (MQ=255)
                            tGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTa                                                                                                                    <  1:3490089/150‑1 (MQ=255)
                              cGGATTGGCGGCGTCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGATACCTAACTGATTGTTTTCTACTGTATAGCCAAGCTAATAATCGTTATTTCTCACCGTATAAGTTAATTTTCTTAAGTCCTCAGTGTTGGGGTTTTg                                                                                                                  <  2:606010/150‑1 (MQ=255)
                                gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                >  2:1298336/1‑150 (MQ=255)
                                gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                <  1:4139757/150‑1 (MQ=255)
                                 aTTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCgg                                                                                                               <  2:2056334/150‑1 (MQ=255)
                                  ttGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCggg                                                                                                              >  1:2930889/1‑150 (MQ=255)
                                  ttGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCggg                                                                                                              <  2:3327315/150‑1 (MQ=255)
                                  ttGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCggg                                                                                                              >  1:3808451/1‑150 (MQ=255)
                                    ggcggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTa                                                                                                            >  1:2569938/1‑150 (MQ=255)
                                    ggcggcggcAAGGTCATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTa                                                                                                            >  1:3052683/1‑150 (MQ=255)
                                     gcggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTaa                                                                                                           <  2:2446637/150‑1 (MQ=255)
                                       ggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATc                                                                                                         >  2:429579/1‑150 (MQ=255)
                                       ggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATc                                                                                                         <  1:2187157/150‑1 (MQ=255)
                                            cAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctg                                                                                                    <  2:2542629/150‑1 (MQ=255)
                                            cAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctg                                                                                                    <  2:434453/150‑1 (MQ=255)
                                             aaGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTATACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctta                                                                                                   <  1:1183967/150‑3 (MQ=255)
                                              aGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgac                                                                                                  >  1:2807259/1‑150 (MQ=255)
                                               ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                 <  2:2496267/150‑1 (MQ=255)
                                               ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                 >  2:1669672/1‑150 (MQ=255)
                                                gTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactg                                                                                                >  2:3748332/1‑150 (MQ=255)
                                                  tATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGAt                                                                                              <  1:3621803/150‑1 (MQ=255)
                                                     cTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGc                                                                                           >  2:3512594/1‑150 (MQ=255)
                                                         ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                       <  2:2797941/150‑1 (MQ=255)
                                                               tAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAttt                                                                                 <  2:1371950/150‑1 (MQ=255)
                                                               tAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAttt                                                                                 <  2:4141187/150‑1 (MQ=255)
                                                                  aTCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgtg                                                                              >  2:4476920/1‑150 (MQ=255)
                                                                  aTCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgtg                                                                              <  1:3001757/150‑1 (MQ=255)
                                                                   tCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTgg                                                                             <  2:2000026/150‑1 (MQ=255)
                                                                      gTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTaa                                                                          <  1:2367789/150‑1 (MQ=255)
                                                                      gTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTaa                                                                          >  1:212931/1‑150 (MQ=255)
                                                                        cAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAtt                                                                        >  1:2623248/1‑150 (MQ=255)
                                                                        cAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAtt                                                                        >  2:3961562/1‑150 (MQ=255)
                                                                            ttGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATcc                                                                    <  2:1845436/150‑1 (MQ=255)
                                                                              gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                  <  2:117779/150‑1 (MQ=255)
                                                                                     ccTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGACCAAATTTGTGGTAATTATCCCAGAATGGa                                                           >  1:4457356/1‑150 (MQ=255)
                                                                                     ccTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGa                                                           <  2:530830/150‑1 (MQ=255)
                                                                                       tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                         >  2:486035/1‑150 (MQ=255)
                                                                                       tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                         >  2:1736855/1‑150 (MQ=255)
                                                                                          cTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAtt                                                      <  1:3675191/150‑1 (MQ=255)
                                                                                          cTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCAAGAATGGAATAtt                                                      >  2:1953115/1‑150 (MQ=255)
                                                                                             aTTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTc                                                   >  2:1353070/1‑150 (MQ=255)
                                                                                              ttGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTcc                                                  <  1:654605/150‑1 (MQ=255)
                                                                                               tGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCt                                                 <  1:1039137/150‑1 (MQ=255)
                                                                                               tGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCt                                                 >  2:460554/1‑150 (MQ=255)
                                                                                                        cTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAcc                                        <  2:932453/150‑1 (MQ=255)
                                                                                                         tGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCa                                       <  1:473061/150‑1 (MQ=255)
                                                                                                           tataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGt                                     <  1:693906/150‑1 (MQ=255)
                                                                                                           tataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGt                                     <  1:1544271/150‑1 (MQ=255)
                                                                                                             taGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTga                                   >  1:1406204/1‑150 (MQ=255)
                                                                                                                aCAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGaa                                <  1:843429/150‑1 (MQ=255)
                                                                                                                  aaGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                              >  2:1038685/1‑150 (MQ=255)
                                                                                                                  aaGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                              >  1:1796620/1‑150 (MQ=255)
                                                                                                                    gCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAag                            <  2:3580216/150‑1 (MQ=255)
                                                                                                                        ataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGt                        >  1:2426061/1‑150 (MQ=255)
                                                                                                                          aatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTa                      >  2:1129722/1‑150 (MQ=255)
                                                                                                                           atCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTAt                     <  1:1876796/150‑1 (MQ=255)
                                                                                                                           atCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTAt                     >  2:894128/1‑150 (MQ=255)
                                                                                                                           atCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTAt                     >  2:2853831/1‑150 (MQ=255)
                                                                                                                             cGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATcc                   >  2:415705/1‑150 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                  <  2:2876918/150‑1 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                  <  1:2583819/150‑1 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                  >  2:4161358/1‑150 (MQ=255)
                                                                                                                               ttATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAg                 <  1:2227010/150‑1 (MQ=255)
                                                                                                                               ttATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAg                 <  1:3330923/150‑1 (MQ=255)
                                                                                                                                tATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAgg                >  2:2180558/1‑150 (MQ=255)
                                                                                                                                   ttCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCa             <  2:3198738/150‑1 (MQ=255)
                                                                                                                                        aCCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt        <  1:876006/150‑1 (MQ=255)
                                                                                                                                         ccGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtttt       >  1:3615476/1‑150 (MQ=255)
                                                                                                                                         ccGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtttt       >  1:145599/1‑150 (MQ=255)
                                                                                                                                          cGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttttt      >  2:2701973/1‑150 (MQ=255)
                                                                                                                                           gTATAATATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCATAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa     <  1:3974772/150‑1 (MQ=255)
                                                                                                                                           gTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa     <  1:1604542/150‑1 (MQ=255)
                                                                                                                                           gTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa     <  1:1081306/150‑1 (MQ=255)
                                                                                                                                           gTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa     >  1:907428/1‑150 (MQ=255)
                                                                                                                                            tataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAt    <  2:3584421/150‑1 (MQ=255)
                                                                                                                                             ataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtg   >  2:3805214/1‑150 (MQ=255)
                                                                                                                                              taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt  <  1:4075236/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGT  >  CP009273/1575192‑1575483

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: