Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,575,339 G→T 69.1% 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,3390GT69.1% 266.9 / 229.3 210F481L (TTC→TTAydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (35/30);  new base T (75/70);  total (110/100)
Fisher's exact test for biased strand distribution p-value = 8.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.16e-01

Only 100 of 211 total aligned reads displayed.

TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGGC  >  CP009273/1575192‑1575486
                                                                                                                                                   |                                                                                                                                                   
    aCGCCTTGCATCTCTTTACCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTgg                                                                                                                                               <  1:2697658/150‑1 (MQ=255)
            cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                       >  2:4555881/1‑150 (MQ=255)
              tctcTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCt                                                                                                                                     <  1:118926/150‑1 (MQ=255)
               ctctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTc                                                                                                                                    <  2:1534103/150‑1 (MQ=255)
                tctTTCACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCa                                                                                                                                   <  1:1619042/150‑1 (MQ=255)
                  tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                 <  1:2612785/150‑1 (MQ=255)
                  tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                 <  1:185061/150‑1 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                               <  1:4211883/150‑1 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                               <  1:3742020/150‑1 (MQ=255)
                    tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                               >  1:1709141/1‑150 (MQ=255)
                           cTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGtt                                                                                                                        >  2:922245/1‑150 (MQ=255)
                             gCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTAt                                                                                                                      <  2:3513833/150‑1 (MQ=255)
                             gCGGAGTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTAt                                                                                                                      <  2:2555213/150‑1 (MQ=255)
                               ggATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAATCGAATAATCGTTATTTCTCACCGAATAAGATAATTGTCTTAAGTCCTCAGTGTTGGGGTTATGc                                                                                                                    >  1:4621522/1‑150 (MQ=255)
                                gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                   >  2:1134886/1‑150 (MQ=255)
                                       ggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATc                                                                                                            <  2:1909837/150‑1 (MQ=255)
                                        gcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCg                                                                                                           <  2:1395650/150‑1 (MQ=255)
                                            cAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctg                                                                                                       >  1:1783560/1‑150 (MQ=255)
                                             aaGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctga                                                                                                      <  2:3964820/150‑1 (MQ=255)
                                                gTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactg                                                                                                   >  2:4464710/1‑150 (MQ=255)
                                                 ttATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactga                                                                                                  >  1:2147086/1‑150 (MQ=255)
                                                 ttATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactga                                                                                                  >  1:3325907/1‑150 (MQ=255)
                                                   aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                                <  1:1223234/150‑1 (MQ=255)
                                                   aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                                >  2:2100641/1‑150 (MQ=255)
                                                   aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                                <  1:3038204/150‑1 (MQ=255)
                                                   aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                                >  1:2796107/1‑150 (MQ=255)
                                                    tCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATgg                                                                                               <  1:4558290/150‑1 (MQ=255)
                                                        tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                           >  2:3122103/1‑150 (MQ=255)
                                                        tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                           >  2:159632/1‑150 (MQ=255)
                                                              gTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAtt                                                                                     >  2:2691171/1‑150 (MQ=255)
                                                               tAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAttt                                                                                    <  1:3706228/150‑1 (MQ=255)
                                                               tAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAttt                                                                                    >  1:1207394/1‑150 (MQ=255)
                                                                    ccGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGt                                                                               >  2:4794177/1‑150 (MQ=255)
                                                                    ccGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGt                                                                               >  2:742057/1‑150 (MQ=255)
                                                                     cGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTa                                                                              >  2:1865127/1‑150 (MQ=255)
                                                                      gTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTaa                                                                             >  1:1237363/1‑150 (MQ=255)
                                                                      gTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCATATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTaa                                                                             <  1:4872767/150‑1 (MQ=255)
                                                                       tCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTTTGGTAAt                                                                            >  2:7642/1‑150 (MQ=255)
                                                                       tCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAt                                                                            <  2:3287454/150‑1 (MQ=255)
                                                                       tCAGTGTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAAt                                                                            <  1:2723778/150‑1 (MQ=255)
                                                                         aGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTa                                                                          <  2:4335636/150‑1 (MQ=255)
                                                                          gTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTAt                                                                         <  1:4057772/150‑1 (MQ=255)
                                                                           tttGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATc                                                                        <  1:3569453/150‑1 (MQ=255)
                                                                             tGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATccc                                                                      >  1:768851/1‑150 (MQ=255)
                                                                             tGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATccc                                                                      <  1:2365020/150‑1 (MQ=255)
                                                                              gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                     <  1:1424068/150‑1 (MQ=255)
                                                                              gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                     <  1:1422079/150‑1 (MQ=255)
                                                                                agagACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGa                                                                   >  1:2073168/1‑150 (MQ=255)
                                                                                  agaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAAt                                                                 <  1:348241/150‑1 (MQ=255)
                                                                                   gaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATg                                                                >  1:3345377/1‑150 (MQ=255)
                                                                                   gaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATg                                                                >  2:284041/1‑150 (MQ=255)
                                                                                   gaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATg                                                                <  1:3107511/150‑1 (MQ=255)
                                                                                     ccTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGa                                                              <  1:1607717/150‑1 (MQ=255)
                                                                                        aaCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAata                                                           <  1:2414558/150‑1 (MQ=255)
                                                                                        aaCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAata                                                           <  1:1205937/150‑1 (MQ=255)
                                                                                         aCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAatat                                                          >  1:4836266/1‑150 (MQ=255)
                                                                                         aCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAatat                                                          >  1:4705123/1‑150 (MQ=255)
                                                                                           tGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAttt                                                        <  2:4332502/150‑1 (MQ=255)
                                                                                            gATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAtttt                                                       >  2:2643279/1‑150 (MQ=255)
                                                                                                gTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                   <  1:4168029/150‑1 (MQ=255)
                                                                                                gTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                   >  2:2388487/1‑150 (MQ=255)
                                                                                                 ttttCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCg                                                  <  2:4247524/150‑1 (MQ=255)
                                                                                                    tCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCa                                               <  2:4185755/150‑1 (MQ=255)
                                                                                                     cTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCaa                                              <  1:3968851/150‑1 (MQ=255)
                                                                                                     cTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCaa                                              >  1:3908058/1‑150 (MQ=255)
                                                                                                       aCTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAc                                            <  1:1211409/150‑1 (MQ=255)
                                                                                                       aCTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAc                                            >  1:3413045/1‑150 (MQ=255)
                                                                                                       aCTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAc                                            >  1:3411344/1‑150 (MQ=255)
                                                                                                        cTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAcc                                           <  2:1499950/150‑1 (MQ=255)
                                                                                                           tataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGt                                        <  2:4284394/150‑1 (MQ=255)
                                                                                                            ataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTg                                       <  1:928894/150‑1 (MQ=255)
                                                                                                            ataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTg                                       >  2:3734938/1‑150 (MQ=255)
                                                                                                                 cAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAAt                                  <  1:4353703/150‑1 (MQ=255)
                                                                                                                 cAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAAt                                  <  1:1028141/150‑1 (MQ=255)
                                                                                                                  aaGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGACCAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                                 >  1:1758040/1‑150 (MQ=255)
                                                                                                                    gCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAag                               <  1:4191380/150‑1 (MQ=255)
                                                                                                                     cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                              >  2:3210092/1‑150 (MQ=255)
                                                                                                                     cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                              >  2:779508/1‑150 (MQ=255)
                                                                                                                       aataatCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAgg                            >  2:3307727/1‑150 (MQ=255)
                                                                                                                         taatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGtt                          >  2:481639/1‑150 (MQ=255)
                                                                                                                         taatCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGtt                          <  1:3899723/150‑1 (MQ=255)
                                                                                                                          aatCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTa                         <  2:4745239/150‑1 (MQ=255)
                                                                                                                             cGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATcc                      >  2:3359038/1‑150 (MQ=255)
                                                                                                                              gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                     >  1:1326729/1‑150 (MQ=255)
                                                                                                                               ttATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAg                    >  1:2376651/1‑150 (MQ=255)
                                                                                                                                tATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAgg                   <  2:426143/150‑1 (MQ=255)
                                                                                                                                  tttCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTc                 >  2:2248521/1‑150 (MQ=255)
                                                                                                                                   ttCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCa                >  1:256852/1‑150 (MQ=255)
                                                                                                                                   ttCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCa                >  1:3438596/1‑150 (MQ=255)
                                                                                                                                      tcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGGGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAt             <  2:1748550/150‑1 (MQ=255)
                                                                                                                                       cACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtt            >  2:1162958/1‑150 (MQ=255)
                                                                                                                                         ccGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtttt          <  2:3785685/150‑1 (MQ=255)
                                                                                                                                          cGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttttt         >  1:17575/1‑150 (MQ=255)
                                                                                                                                            tataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAt       >  2:243700/1‑150 (MQ=255)
                                                                                                                                             ataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtg      <  1:2626350/150‑1 (MQ=255)
                                                                                                                                             ataAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtg      >  2:1277401/1‑150 (MQ=255)
                                                                                                                                              taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt     <  1:605877/150‑1 (MQ=255)
                                                                                                                                              taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt     >  1:4722903/1‑150 (MQ=255)
                                                                                                                                               aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg    >  1:743059/1‑150 (MQ=255)
                                                                                                                                                 gagaATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGACCAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTATTATGTGGc  >  2:2125344/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGGC  >  CP009273/1575192‑1575486

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: