Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,575,339 G→T 92.5% 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,3390GT92.5% 777.9 / 47.7 253F481L (TTC→TTAydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (9/10);  new base T (118/116);  total (127/126)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 256 total aligned reads displayed.

ACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575193‑1575485
                                                                                                                                                  |                                                                                                                                                  
aCCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAAt                                                                                                                                                 <  2:1565321/150‑1 (MQ=255)
aCCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAAt                                                                                                                                                 >  1:1173512/1‑150 (MQ=255)
aCCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAAt                                                                                                                                                 <  2:4483600/150‑1 (MQ=255)
     gCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCt                                                                                                                                            >  2:3818427/1‑150 (MQ=255)
       cTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTa                                                                                                                                          <  1:4135692/150‑1 (MQ=255)
         tGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAg                                                                                                                                        >  2:4581150/1‑150 (MQ=255)
         tGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAg                                                                                                                                        >  2:2180074/1‑150 (MQ=255)
          gCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGt                                                                                                                                       <  1:2663954/150‑1 (MQ=255)
           cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                      >  2:1009806/1‑150 (MQ=255)
           cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                      >  2:4109916/1‑150 (MQ=255)
           cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                      >  1:241572/1‑150 (MQ=255)
           cATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTc                                                                                                                                      >  2:1757703/1‑150 (MQ=255)
              ctctTTAACGACCTGCGGATTGGCGGCGGCAATGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTc                                                                                                                                   <  2:1049492/150‑1 (MQ=255)
               tctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCa                                                                                                                                  >  2:2298368/1‑150 (MQ=255)
                ctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAg                                                                                                                                 >  2:3664303/1‑150 (MQ=255)
                 tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                >  1:2778332/1‑150 (MQ=255)
                 tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                >  2:1456406/1‑150 (MQ=255)
                   tAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                              <  1:1316545/150‑1 (MQ=255)
                      cGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgg                                                                                                                           >  1:4014595/1‑150 (MQ=255)
                      cGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgg                                                                                                                           >  1:1446072/1‑150 (MQ=255)
                       gACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTggg                                                                                                                          <  2:3924023/150‑1 (MQ=255)
                       gACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTggg                                                                                                                          <  2:3924291/150‑1 (MQ=255)
                       gACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTggg                                                                                                                          >  1:670639/1‑150 (MQ=255)
                       gACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTggg                                                                                                                          >  2:460745/1‑150 (MQ=255)
                        aCCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgggg                                                                                                                         <  1:4614667/150‑1 (MQ=255)
                         ccTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGt                                                                                                                        >  2:3251385/1‑150 (MQ=255)
                          cTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGtt                                                                                                                       >  2:1895165/1‑150 (MQ=255)
                              ggATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGc                                                                                                                   >  1:3145277/1‑150 (MQ=255)
                               gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                  <  2:4022621/150‑1 (MQ=255)
                               gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                  >  1:1063709/1‑150 (MQ=255)
                                aTTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCgg                                                                                                                 >  1:3355839/1‑150 (MQ=255)
                                aTTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCgg                                                                                                                 >  1:2483325/1‑150 (MQ=255)
                                 ttGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCggg                                                                                                                >  1:3410953/1‑150 (MQ=255)
                                     cggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAAt                                                                                                            >  2:2637005/1‑150 (MQ=255)
                                     cggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAAt                                                                                                            <  1:2119752/150‑1 (MQ=255)
                                      ggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATc                                                                                                           >  1:801081/1‑150 (MQ=255)
                                           cAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctg                                                                                                      <  1:4704629/150‑1 (MQ=255)
                                              ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   >  2:3337984/1‑150 (MQ=255)
                                              ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   <  1:3016481/150‑1 (MQ=255)
                                                ttATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactga                                                                                                 <  1:990365/150‑1 (MQ=255)
                                                  aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                               >  2:2730517/1‑150 (MQ=255)
                                                      ttttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGa                                                                                           <  2:4445717/150‑1 (MQ=255)
                                                       tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                          >  1:1823189/1‑150 (MQ=255)
                                                            tGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAt                                                                                     >  1:1862996/1‑150 (MQ=255)
                                                             gTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAAtt                                                                                    <  1:52073/150‑1 (MQ=255)
                                                               aGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtg                                                                                  <  2:4625379/150‑1 (MQ=255)
                                                                 aTCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgtg                                                                                >  1:4771041/1‑150 (MQ=255)
                                                                  tCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTgg                                                                               <  2:4098437/150‑1 (MQ=255)
                                                                        aGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTa                                                                         <  1:965633/150‑1 (MQ=255)
                                                                           ttGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATcc                                                                      <  1:3453322/150‑1 (MQ=255)
                                                                           ttGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATcc                                                                      <  1:3034498/150‑1 (MQ=255)
                                                                            tGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATccc                                                                     <  1:2289064/150‑1 (MQ=255)
                                                                             gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                    >  2:1228622/1‑150 (MQ=255)
                                                                             gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                    <  2:2362400/150‑1 (MQ=255)
                                                                               agagACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGa                                                                  <  2:3610600/150‑1 (MQ=255)
                                                                                gagaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGaa                                                                 >  2:1277267/1‑150 (MQ=255)
                                                                                 agaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAAt                                                                <  1:4453792/150‑1 (MQ=255)
                                                                                  gaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATg                                                               <  2:1397227/150‑1 (MQ=255)
                                                                                   aCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATgg                                                              >  1:1802377/1‑150 (MQ=255)
                                                                                    ccTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGa                                                             >  1:2753653/1‑150 (MQ=255)
                                                                                       aaCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAata                                                          >  1:4552480/1‑150 (MQ=255)
                                                                                          tGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAttt                                                       <  2:2125001/150‑1 (MQ=255)
                                                                                           gATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAtttt                                                      <  2:2778332/150‑1 (MQ=255)
                                                                                              tGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCt                                                   >  1:561420/1‑150 (MQ=255)
                                                                                               gtgttCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                  <  2:3705568/147‑1 (MQ=255)
                                                                                               gTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTc                                                  >  2:3433015/1‑150 (MQ=255)
                                                                                                  ttCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGACCAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTc                                               >  1:1002678/1‑150 (MQ=255)
                                                                                                    cTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCaa                                             <  2:1911958/150‑1 (MQ=255)
                                                                                                     tACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCaaa                                            <  2:566298/150‑1 (MQ=255)
                                                                                                      aCTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAc                                           >  1:3766137/1‑150 (MQ=255)
                                                                                                      aCTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAAc                                           <  2:3144235/150‑1 (MQ=255)
                                                                                                           ataGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTg                                      <  1:1741946/150‑1 (MQ=255)
                                                                                                             aGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgag                                    <  1:195131/150‑1 (MQ=255)
                                                                                                             aGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgag                                    <  2:4003175/150‑1 (MQ=255)
                                                                                                              gACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgaga                                   <  2:1832228/150‑1 (MQ=255)
                                                                                                                cAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAAt                                 >  2:858879/1‑150 (MQ=255)
                                                                                                                 aaGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                                >  1:2369626/1‑150 (MQ=255)
                                                                                                                    cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                             >  2:369355/1‑150 (MQ=255)
                                                                                                                    cGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAaga                             >  1:4219493/1‑150 (MQ=255)
                                                                                                                     gAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAagag                            >  2:3406223/1‑150 (MQ=255)
                                                                                                                      tttaatCTTTATTTCTCACCGTTTAAGATTTTTGGGTTAAGTGCTCATTGTTTGGGTTATGCTGGTAATCTTCTGACTGATTGCGAACTAATTTGTGGTAATTATCCCAGAATGGAATTTTTTCCTAGTCTAACCATTGTTAATAagatg                           <  2:2467833/148‑3 (MQ=255)
                                                                                                                      aataatCGTTATTTCTCACCGTATAAGATAATTGGCTTGAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAgg                           >  1:2772466/1‑150 (MQ=255)
                                                                                                                      aataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAgg                           <  2:924388/150‑1 (MQ=255)
                                                                                                                       ataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGt                          <  1:3664303/150‑1 (MQ=255)
                                                                                                                        taatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGtt                         >  2:3550711/1‑150 (MQ=255)
                                                                                                                         aatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTa                        <  1:1608486/150‑1 (MQ=255)
                                                                                                                                 tttCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTc                >  1:2758524/1‑150 (MQ=255)
                                                                                                                                 tttCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTc                <  2:1173512/150‑1 (MQ=255)
                                                                                                                                    ctcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGa             >  2:3394984/1‑150 (MQ=255)
                                                                                                                                     tcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAt            >  1:690980/1‑150 (MQ=255)
                                                                                                                                      cACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAtt           >  1:864934/1‑150 (MQ=255)
                                                                                                                                       aCCGTATAAGATTATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGGGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt          <  1:1757703/150‑1 (MQ=255)
                                                                                                                                       aCCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGACTGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt          >  2:4689390/1‑150 (MQ=255)
                                                                                                                                       aCCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt          >  1:4284154/1‑150 (MQ=255)
                                                                                                                                       aCCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAttt          <  1:4109916/150‑1 (MQ=255)
                                                                                                                                          gTATAAGATAATGGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTa       <  1:717276/150‑1 (MQ=255)
                                                                                                                                           tataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAt      <  1:3366050/150‑1 (MQ=255)
                                                                                                                                           tataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAt      >  2:984368/1‑150 (MQ=255)
                                                                                                                                             taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt    >  2:3736978/1‑150 (MQ=255)
                                                                                                                                             taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt    >  2:1637075/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACCACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575193‑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: