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,3390GT96.8% 822.6 / 14.5 246F481L (TTC→TTAydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (5/3);  new base T (100/138);  total (105/141)
Fisher's exact test for biased strand distribution p-value = 2.92e-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 250 total aligned reads displayed.

TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575192‑1575485
                                                                                                                                                   |                                                                                                                                                  
 aCCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAt                                                                                                                                                 >  2:911030/1‑150 (MQ=255)
  ccACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAtt                                                                                                                                                <  1:2505366/150‑1 (MQ=255)
  ccACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAtt                                                                                                                                                <  1:1745374/150‑1 (MQ=255)
    aCGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTgg                                                                                                                                              <  2:830065/150‑1 (MQ=255)
    aCGACTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGTTTATCTTTTTGCTGTAGATCCGTCTGTTTGTAGAGACCTAACTGATTGTTTTCTTCTGTTTAGTCTAGCGATTAATCGTTATTTCTCACCGTATAAGATAATTgg                                                                                                                                              <  2:550358/150‑1 (MQ=255)
      gCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCt                                                                                                                                            <  2:2848619/150‑1 (MQ=255)
        cTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTa                                                                                                                                          >  1:3598820/1‑150 (MQ=255)
             aTCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTcc                                                                                                                                     >  2:427804/1‑150 (MQ=255)
                 ctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCTTATAAGATAATTGGCTTAAGTCCTCAg                                                                                                                                 <  1:195731/150‑1 (MQ=255)
                  tttCACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                <  1:3910840/150‑1 (MQ=255)
                  tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                <  1:3046529/150‑1 (MQ=255)
                   ttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtg                                                                                                                               >  1:3083276/1‑150 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                              >  1:4338597/1‑150 (MQ=255)
                     aaCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGtt                                                                                                                             <  1:3741410/150‑1 (MQ=255)
                     aaCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGtt                                                                                                                             >  1:3170754/1‑150 (MQ=255)
                     aaCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGtt                                                                                                                             >  1:2487988/1‑150 (MQ=255)
                      aCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTg                                                                                                                            <  2:3299420/150‑1 (MQ=255)
                       cGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgg                                                                                                                           <  1:1936224/150‑1 (MQ=255)
                        gACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGACTAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTggg                                                                                                                          >  2:3509752/1‑150 (MQ=255)
                         aCCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgggg                                                                                                                         <  2:1025190/150‑1 (MQ=255)
                         aCCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgggg                                                                                                                         >  2:2271879/1‑150 (MQ=255)
                          ccTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAGTGGCTTAAGTCCTCAGTGTTGGGGt                                                                                                                        >  1:757296/1‑150 (MQ=255)
                           cTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGtt                                                                                                                       >  2:3887409/1‑150 (MQ=255)
                                    ggcggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTa                                                                                                              <  1:320393/150‑1 (MQ=255)
                                      cggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAAt                                                                                                            <  1:531329/150‑1 (MQ=255)
                                        gcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCg                                                                                                          <  2:2061004/150‑1 (MQ=255)
                                           gcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTct                                                                                                       <  2:2780560/150‑1 (MQ=255)
                                               ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   <  2:2585128/150‑1 (MQ=255)
                                               ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   >  1:3993979/1‑150 (MQ=255)
                                                gTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactg                                                                                                  >  1:1019673/1‑150 (MQ=255)
                                                 ttATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactga                                                                                                 <  1:4513487/150‑1 (MQ=255)
                                                  tATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGAt                                                                                                <  1:2088834/150‑1 (MQ=255)
                                                     cTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGc                                                                                             >  1:2900026/1‑150 (MQ=255)
                                                       ttttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGa                                                                                           <  1:71952/150‑1 (MQ=255)
                                                        tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                          <  2:1029909/150‑1 (MQ=255)
                                                         ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                         <  2:2519366/150‑1 (MQ=255)
                                                         ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                         <  1:164252/150‑1 (MQ=255)
                                                              gTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAtt                                                                                    <  2:2173777/150‑1 (MQ=255)
                                                              gTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAtt                                                                                    <  2:2820832/150‑1 (MQ=255)
                                                               tAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAttt                                                                                   <  2:3182317/150‑1 (MQ=255)
                                                                    ccGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGt                                                                              >  2:2145443/1‑150 (MQ=255)
                                                                    ccGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGt                                                                              >  2:2145483/1‑150 (MQ=255)
                                                                       tCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAt                                                                           <  2:1347607/150‑1 (MQ=255)
                                                                       tCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAt                                                                           <  2:1652090/150‑1 (MQ=255)
                                                                        cAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAtt                                                                          >  2:3964205/1‑150 (MQ=255)
                                                                        cAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAtt                                                                          <  2:116158/150‑1 (MQ=255)
                                                                             tGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATccc                                                                     <  2:499185/150‑1 (MQ=255)
                                                                              gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                    <  1:2638837/150‑1 (MQ=255)
                                                                               tAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAg                                                                   >  1:2942772/1‑150 (MQ=255)
                                                                                agagACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGa                                                                  >  1:104198/1‑150 (MQ=255)
                                                                                agagACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGa                                                                  <  1:1216989/150‑1 (MQ=255)
                                                                                  agaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAAt                                                                >  2:1785576/1‑150 (MQ=255)
                                                                                  agaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAAt                                                                <  1:2102693/150‑1 (MQ=255)
                                                                                   gaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATg                                                               >  2:4569006/1‑150 (MQ=255)
                                                                                       tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                           >  2:4335283/1‑150 (MQ=255)
                                                                                       tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                           <  2:4239573/150‑1 (MQ=255)
                                                                                          cTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAtt                                                        <  1:3359841/150‑1 (MQ=255)
                                                                                           tGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAttt                                                       <  1:2278281/150‑1 (MQ=255)
                                                                                              ttGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTcc                                                    <  1:2848412/150‑1 (MQ=255)
                                                                                                gTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                  <  1:182992/150‑1 (MQ=255)
                                                                                                gTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                  <  2:2370134/150‑1 (MQ=255)
                                                                                                  tttCTACAGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGt                                                >  2:193156/1‑150 (MQ=255)
                                                                                                        cTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAcc                                          >  2:3053348/1‑150 (MQ=255)
                                                                                                         tGTATAGACTAGCGAATAATCTTTATTTCTCACCTTATAAGATAATTGGGTTAAGTCCTCAGTTTTGGGGTTTTGCGGTTAATCGTCTGACTGATGGAGAACAAATTTGTGTTAATTATCCCAGAATGTAATATTTTCCTCGTCAAACCa                                         <  2:3712204/150‑1 (MQ=255)
                                                                                                         tGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCa                                         >  2:4107682/1‑150 (MQ=255)
                                                                                                         tGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCa                                         >  1:4510180/1‑150 (MQ=255)
                                                                                                          gTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAg                                        >  1:2990656/1‑150 (MQ=255)
                                                                                                           tataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGt                                       >  1:1561662/1‑150 (MQ=255)
                                                                                                              aGACAATCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCATCGTCAAACCAGTgag                                    >  2:3000817/1‑150 (MQ=255)
                                                                                                               gACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgaga                                   >  1:2759544/1‑150 (MQ=255)
                                                                                                                aCAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGaa                                  >  2:1825009/1‑150 (MQ=255)
                                                                                                                   aGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATaa                               <  1:1462455/150‑1 (MQ=255)
                                                                                                                     cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                             >  1:186914/1‑150 (MQ=255)
                                                                                                                     cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                             >  2:3998601/1‑150 (MQ=255)
                                                                                                                      gAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAagag                            >  1:740687/1‑150 (MQ=255)
                                                                                                                      gAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAagag                            >  1:424180/1‑150 (MQ=255)
                                                                                                                        ataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGt                          <  1:1639774/150‑1 (MQ=255)
                                                                                                                         taatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGtt                         <  1:1621662/150‑1 (MQ=255)
                                                                                                                             cGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATcc                     >  1:4210805/1‑150 (MQ=255)
                                                                                                                             cGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATcc                     <  2:4531121/150‑1 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                    >  2:622735/1‑150 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                    >  1:122227/1‑150 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                    <  2:4324610/150‑1 (MQ=255)
                                                                                                                               ttATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAg                   <  2:2662483/150‑1 (MQ=255)
                                                                                                                               ttATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAg                   >  2:466322/1‑150 (MQ=255)
                                                                                                                                tATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAgg                  >  2:2430396/1‑150 (MQ=255)
                                                                                                                                tATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGTTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAgg                  >  2:1634476/1‑150 (MQ=255)
                                                                                                                                 aTTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGt                 >  2:3819415/1‑150 (MQ=255)
                                                                                                                                 aTTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGt                 >  2:3464786/1‑150 (MQ=255)
                                                                                                                                    tctcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTGATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAg              >  2:3713205/1‑150 (MQ=255)
                                                                                                                                     ctcACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGa             <  2:783186/150‑1 (MQ=255)
                                                                                                                                       cACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtt           <  2:4496891/150‑1 (MQ=255)
                                                                                                                                        aCCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt          >  1:2327102/1‑150 (MQ=255)
                                                                                                                                          cGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttttt        <  2:4416676/150‑1 (MQ=255)
                                                                                                                                          cGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttttt        <  2:889918/150‑1 (MQ=255)
                                                                                                                                           gTAGAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa       <  2:4043954/150‑1 (MQ=255)
                                                                                                                                               aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg   <  1:3678107/150‑1 (MQ=255)
                                                                                                                                               aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg   >  1:2054286/1‑150 (MQ=255)
                                                                                                                                               aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg   <  1:4407036/150‑1 (MQ=255)
                                                                                                                                                agatAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTgg  >  1:1509830/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575192‑1575485

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: