Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,575,339 G→T 93.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,3390GT93.5% 765.5 / 40.5 245F481L (TTC→TTAydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (5/11);  new base T (104/125);  total (109/136)
Fisher's exact test for biased strand distribution p-value = 3.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.48e-01

Only 100 of 247 total aligned reads displayed.

ACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575196‑1575485
                                                                                                                                               |                                                                                                                                                  
aCGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTgg                                                                                                                                              <  1:55655/150‑1 (MQ=255)
  gCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCt                                                                                                                                            <  2:1361167/150‑1 (MQ=255)
    cTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTa                                                                                                                                          <  1:3004343/150‑1 (MQ=255)
    cTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTa                                                                                                                                          >  2:195404/1‑150 (MQ=255)
      tGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAg                                                                                                                                        >  2:2822963/1‑150 (MQ=255)
        taTCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATACTTGGCTTAAGTc                                                                                                                                      >  1:1378155/2‑150 (MQ=255)
         aTCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTcc                                                                                                                                     >  1:1060912/1‑150 (MQ=255)
         aTCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTcc                                                                                                                                     >  2:1835647/1‑150 (MQ=255)
          tctcTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCt                                                                                                                                    >  2:1336006/1‑150 (MQ=255)
            tctTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCa                                                                                                                                  >  2:954664/1‑150 (MQ=255)
              tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCTGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATCAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                >  2:864398/1‑150 (MQ=255)
              tttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgt                                                                                                                                >  2:3041983/1‑150 (MQ=255)
               ttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtg                                                                                                                               >  1:730637/1‑150 (MQ=255)
               ttAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtg                                                                                                                               <  2:3768809/150‑1 (MQ=255)
                tAACGACCTGCGGATTGGCGGCGTCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTTTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAgtgt                                                                                                                              <  2:4087181/150‑1 (MQ=255)
                 aaCGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGtt                                                                                                                             <  1:2795180/150‑1 (MQ=255)
                     aCCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTgggg                                                                                                                         <  2:3098300/150‑1 (MQ=255)
                       cTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGtt                                                                                                                       >  2:261649/1‑150 (MQ=255)
                        tGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTa                                                                                                                      <  2:3001889/150‑1 (MQ=255)
                        tGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTa                                                                                                                      <  1:4326007/150‑1 (MQ=255)
                        tGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTa                                                                                                                      >  2:725071/1‑150 (MQ=255)
                        tGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTa                                                                                                                      >  2:2896129/1‑150 (MQ=255)
                         gCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTAt                                                                                                                     <  1:2857108/150‑1 (MQ=255)
                          cGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATg                                                                                                                    <  2:3614736/150‑1 (MQ=255)
                            gATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCg                                                                                                                  >  2:2472773/1‑150 (MQ=255)
                              ttGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCggg                                                                                                                >  2:745917/1‑150 (MQ=255)
                                  cggcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGAGTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAAt                                                                                                            <  2:393722/150‑1 (MQ=255)
                                    gcggcAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCg                                                                                                          <  1:3076035/150‑1 (MQ=255)
                                        cAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctg                                                                                                      <  2:1037538/150‑1 (MQ=255)
                                         aaGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctga                                                                                                     >  2:1566372/1‑150 (MQ=255)
                                          aGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgac                                                                                                    <  1:1738426/150‑1 (MQ=255)
                                          aGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgaa                                                                                                    <  2:773302/150‑2 (MQ=255)
                                           ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   <  2:134850/150‑1 (MQ=255)
                                           ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   >  2:3650073/1‑150 (MQ=255)
                                           ggTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgact                                                                                                   <  2:2271119/150‑1 (MQ=255)
                                            gTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTctgactg                                                                                                  >  2:3160991/1‑150 (MQ=255)
                                               aTCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATg                                                                                               <  2:792340/150‑1 (MQ=255)
                                                tCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATgg                                                                                              <  1:111349/150‑1 (MQ=255)
                                                 cTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGc                                                                                             <  2:3429391/150‑1 (MQ=255)
                                                  tttttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCg                                                                                            <  2:3958173/150‑1 (MQ=255)
                                                   ttttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGa                                                                                           <  2:455886/150‑1 (MQ=255)
                                                   ttttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGa                                                                                           <  2:819944/150‑1 (MQ=255)
                                                   ttttGCTGTAGATCCGACAGTTTGTAGAGACCTAACTGATTGTCTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGa                                                                                           <  2:2339523/150‑1 (MQ=255)
                                                    tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                          <  1:1628288/150‑1 (MQ=255)
                                                    tttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGaa                                                                                          <  1:1608927/150‑1 (MQ=255)
                                                     ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                         <  2:4223418/150‑1 (MQ=255)
                                                     ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                         <  2:1479314/150‑1 (MQ=255)
                                                     ttGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAAc                                                                                         >  1:2100264/1‑150 (MQ=255)
                                                       gCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACaa                                                                                       >  1:2387289/1‑150 (MQ=255)
                                                       gCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACaa                                                                                       <  1:3653701/150‑1 (MQ=255)
                                                             gATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgt                                                                                 <  1:3285253/150‑1 (MQ=255)
                                                             gATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgt                                                                                 <  2:661044/150‑1 (MQ=255)
                                                              aTCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTtgtg                                                                                <  2:4001091/150‑1 (MQ=255)
                                                               tCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTgg                                                                               <  1:3243949/150‑1 (MQ=255)
                                                                  gTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTaa                                                                            >  2:1175568/1‑150 (MQ=255)
                                                                         tGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATccc                                                                     >  2:3727775/1‑150 (MQ=255)
                                                                          gTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCa                                                                    >  2:4281074/1‑150 (MQ=255)
                                                                           tAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAg                                                                   >  2:124013/1‑150 (MQ=255)
                                                                           tAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAg                                                                   <  2:4358493/150‑1 (MQ=255)
                                                                           tAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAg                                                                   <  1:2455522/150‑1 (MQ=255)
                                                                            agagACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAATTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGa                                                                  >  2:3912586/1‑150 (MQ=255)
                                                                             gagaCCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGaa                                                                 >  1:3557417/1‑150 (MQ=255)
                                                                                 ccTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGa                                                             <  2:1450906/150‑1 (MQ=255)
                                                                                  cTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGaa                                                            <  1:4158063/150‑1 (MQ=255)
                                                                                   tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                           <  2:2100814/150‑1 (MQ=255)
                                                                                   tAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAat                                                           <  2:239570/150‑1 (MQ=255)
                                                                                    aaCTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAata                                                          <  2:2067125/150‑1 (MQ=255)
                                                                                      cTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAtt                                                        <  2:2789026/150‑1 (MQ=255)
                                                                                       tGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATAttt                                                       >  1:1946961/1‑150 (MQ=255)
                                                                                          ttGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTcc                                                    <  2:3575241/150‑1 (MQ=255)
                                                                                          ttGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTcc                                                    >  2:1733958/1‑150 (MQ=255)
                                                                                               ttCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGACCAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTc                                               >  1:3978411/1‑150 (MQ=255)
                                                                                               ttCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTc                                               >  1:3976927/1‑150 (MQ=255)
                                                                                                tCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCa                                              <  2:1195423/150‑1 (MQ=255)
                                                                                                      gTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAg                                        <  2:3331957/150‑1 (MQ=255)
                                                                                                      gTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAg                                        >  1:1679987/1‑150 (MQ=255)
                                                                                                      gTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAg                                        >  1:1601359/1‑150 (MQ=255)
                                                                                                          aGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgag                                    <  1:3063030/150‑1 (MQ=255)
                                                                                                          aGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgag                                    >  2:575405/1‑150 (MQ=255)
                                                                                                           gACAAGCGAATACTCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTgaga                                   <  2:3407337/150‑1 (MQ=255)
                                                                                                             cAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAAt                                 <  2:351489/150‑1 (MQ=255)
                                                                                                             cAAGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAAt                                 <  2:350744/150‑1 (MQ=255)
                                                                                                              aaGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                                >  1:1035207/1‑150 (MQ=255)
                                                                                                              aaGCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATa                                <  2:1081707/150‑1 (MQ=255)
                                                                                                                gCGAATAATCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAag                              <  2:1414447/150‑1 (MQ=255)
                                                                                                                  gAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAagag                            >  1:1024749/1‑150 (MQ=255)
                                                                                                                   aataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAgg                           >  1:1185654/1‑150 (MQ=255)
                                                                                                                   aataatCGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAgg                           <  1:1307901/150‑1 (MQ=255)
                                                                                                                         cGTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGTGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATcc                     <  1:519572/150‑1 (MQ=255)
                                                                                                                          gTTATTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCa                    <  1:1835647/150‑1 (MQ=255)
                                                                                                                             aTTTCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGt                 >  1:768180/1‑150 (MQ=255)
                                                                                                                               ttCTCACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCa               >  1:2024395/1‑150 (MQ=255)
                                                                                                                                 ctcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGa             >  2:1539392/1‑150 (MQ=255)
                                                                                                                                  tcACCGTATAAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGAt            <  1:4055301/150‑1 (MQ=255)
                                                                                                                                         ataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtg     >  2:3953841/1‑150 (MQ=255)
                                                                                                                                         ataAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtg     >  2:3784123/1‑150 (MQ=255)
                                                                                                                                          taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt    <  1:837441/150‑1 (MQ=255)
                                                                                                                                          taAGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgt    <  1:3041808/150‑1 (MQ=255)
                                                                                                                                           aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg   <  1:2617115/150‑1 (MQ=255)
                                                                                                                                           aaGATAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTAtgtg   <  1:2569657/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ACGCCTTGCATCTCTTTAACGACCTGCGGATTGGCGGCGGCAAGGTTATCTTTTTGCTGTAGATCCGTCAGTTTGTAGAGACCTAACTGATTGTTTTCTACTGTATAGACAAGCGAATAATCGTTATTTCTCACCGTATAAGAGAATTGGCTTAAGTCCTCAGTGTTGGGGTTATGCGGGTAATCGTCTGACTGATGGCGAACAAATTTGTGGTAATTATCCCAGAATGGAATATTTTCCTCGTCAAACCAGTGAGAATAAGAGGTTATCCAGGTCAGATTTTTATGTGG  >  CP009273/1575196‑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: