Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,976,568 A→T 5.1% Y84F (TAC→TTC)  aslB → putative AslA‑specific sulfatase‑maturating enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,976,5680AT 5.1% 1263.8 / 41.3 366Y84F (TAC→TTC) aslBputative AslA‑specific sulfatase‑maturating enzyme
Reads supporting (aligned to +/- strand):  ref base A (159/188);  new base T (10/9);  total (169/197)
Fisher's exact test for biased strand distribution p-value = 6.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 374 total aligned reads displayed.

ATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGC  >  CP009273/3976419‑3976717
                                                                                                                                                     |                                                                                                                                                     
aTCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTa                                                                                                                                                       >  1:3211629/1‑150 (MQ=255)
  cccTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAcc                                                                                                                                                     <  2:4943854/150‑1 (MQ=255)
   ccTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccg                                                                                                                                                    >  2:737311/1‑150 (MQ=255)
    cTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgc                                                                                                                                                   >  1:4772675/1‑150 (MQ=255)
     tttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgcc                                                                                                                                                  >  1:874390/1‑150 (MQ=255)
     tttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgcc                                                                                                                                                  <  2:2040038/150‑1 (MQ=255)
     tttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgcc                                                                                                                                                  <  1:3719582/150‑1 (MQ=255)
      ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgccg                                                                                                                                                 >  2:3921355/1‑150 (MQ=255)
      ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgccg                                                                                                                                                 >  2:2417814/1‑150 (MQ=255)
      ttACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTAccgccg                                                                                                                                                 <  2:4641505/150‑1 (MQ=255)
       tACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgt                                                                                                                                                >  2:1283991/1‑150 (MQ=255)
       tACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgt                                                                                                                                                >  2:4369087/1‑150 (MQ=255)
         ccGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgtgc                                                                                                                                              <  1:3059089/150‑1 (MQ=255)
           gcgAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCgccgtgccg                                                                                                                                            >  2:5098602/1‑150 (MQ=255)
             gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTTCCGCCGTGCCGTa                                                                                                                                          <  2:2772610/150‑1 (MQ=255)
             gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTa                                                                                                                                          <  1:4413531/150‑1 (MQ=255)
             gAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTa                                                                                                                                          <  1:3431659/150‑1 (MQ=255)
             gAAAAGACAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGACAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGAGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTa                                                                                                                                          >  1:598097/1‑150 (MQ=255)
                  gcgaGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCt                                                                                                                                     <  1:3455900/147‑1 (MQ=255)
                   ccAGTCACGCATATGGACGATGACTCGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCATCGAACCGCAAAACGATGTTGCTTTTACCTGGCAGGGCGGCGAACCTACGATACTCGGGCTGGCGTTTTTCCTCCGTGCAGTAGCGCTa                                                                                                                                    <  2:638676/150‑1 (MQ=255)
                    cAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTAc                                                                                                                                   >  2:2715521/1‑150 (MQ=255)
                     aGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACa                                                                                                                                  >  1:4687398/1‑150 (MQ=255)
                      gTCACGCATATGGACGATGACATGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAg                                                                                                                                 >  2:4740729/1‑150 (MQ=255)
                      gTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAg                                                                                                                                 >  2:2112486/1‑150 (MQ=255)
                        cACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGc                                                                                                                               >  1:3565381/1‑150 (MQ=255)
                          cGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGa                                                                                                                             >  1:2871579/1‑150 (MQ=255)
                          cGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGa                                                                                                                             >  2:4768546/1‑150 (MQ=255)
                           gCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGaa                                                                                                                            >  2:3850111/1‑150 (MQ=255)
                           gCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGaa                                                                                                                            >  1:1919496/1‑150 (MQ=255)
                            cATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGaaa                                                                                                                           <  2:1628386/150‑1 (MQ=255)
                            cATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGaaa                                                                                                                           <  1:2747972/150‑1 (MQ=255)
                             atatGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAat                                                                                                                          <  2:2472378/150‑1 (MQ=255)
                              tatGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAata                                                                                                                         <  1:3044817/150‑1 (MQ=255)
                                  gACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGt                                                                                                                     <  2:4442588/150‑1 (MQ=255)
                                   aCGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTg                                                                                                                    <  1:2255104/150‑1 (MQ=255)
                                    cGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGc                                                                                                                   <  2:3778933/150‑1 (MQ=255)
                                     gATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCt                                                                                                                  <  2:4251919/150‑1 (MQ=255)
                                      aTGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTg                                                                                                                 >  2:840165/1‑150 (MQ=255)
                                       tGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTgg                                                                                                                <  1:4154041/150‑1 (MQ=255)
                                        gACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGc                                                                                                               <  1:137103/150‑1 (MQ=255)
                                        gACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGc                                                                                                               >  2:3925319/1‑150 (MQ=255)
                                         acacGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCa                                                                                                              <  1:2033901/150‑1 (MQ=255)
                                          cacGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAg                                                                                                             >  2:1551813/1‑150 (MQ=255)
                                                ggAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGAt                                                                                                       >  2:1994897/1‑150 (MQ=255)
                                                 gAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATa                                                                                                      >  2:3520053/1‑150 (MQ=255)
                                                 gAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATa                                                                                                      >  2:905486/1‑150 (MQ=255)
                                                 gAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATa                                                                                                      >  2:1460383/1‑150 (MQ=255)
                                                  aaGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATaa                                                                                                     <  2:2892632/150‑1 (MQ=255)
                                                   aGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTTCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAg                                                                                                    <  1:1846227/150‑1 (MQ=255)
                                                    gCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGt                                                                                                   >  1:3466488/1‑150 (MQ=255)
                                                     cGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTa                                                                                                  <  1:4654581/150‑1 (MQ=255)
                                                       tATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAAc                                                                                                <  1:118892/150‑1 (MQ=255)
                                                         tGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAATTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAg                                                                                              >  1:752692/1‑150 (MQ=255)
                                                         tGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAg                                                                                              >  2:2559962/1‑150 (MQ=255)
                                                                cACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAg                                                                                       <  2:564693/150‑1 (MQ=255)
                                                                   tataTCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACt                                                                                    <  1:1011877/150‑1 (MQ=255)
                                                                   tataTCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTCTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACt                                                                                    <  2:1222399/150‑1 (MQ=255)
                                                                    atatCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTa                                                                                   <  2:1678121/150‑1 (MQ=255)
                                                                    atatCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTGACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTa                                                                                   <  2:5012118/150‑1 (MQ=255)
                                                                     tatCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTaa                                                                                  >  1:208587/1‑150 (MQ=255)
                                                                       tCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACg                                                                                <  2:317037/150‑1 (MQ=255)
                                                                              cAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGtgct                                                                         <  1:3745898/150‑1 (MQ=255)
                                                                               aGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGtgctg                                                                        >  2:775624/1‑150 (MQ=255)
                                                                               aGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGtgctg                                                                        <  2:4811019/150‑1 (MQ=255)
                                                                                 cGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTgctgct                                                                      <  1:803234/150‑1 (MQ=255)
                                                                                 cGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTgctgct                                                                      <  2:325014/150‑1 (MQ=255)
                                                                                 cGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTgctgct                                                                      <  1:1468128/150‑1 (MQ=255)
                                                                                        cAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCTTTTTTCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCCGACTAACGGCGTGCTGCTGGATGAc                                                               >  1:1977183/1‑150 (MQ=255)
                                                                                         aaaaCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACg                                                              >  2:4302655/1‑150 (MQ=255)
                                                                                          aaaCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTTCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGa                                                             >  1:4384809/1‑150 (MQ=255)
                                                                                          aaaCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGa                                                             >  2:1897865/1‑150 (MQ=255)
                                                                                           aaCGACGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGaa                                                            <  1:3921497/150‑1 (MQ=255)
                                                                                           aaCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGaa                                                            <  1:5124493/150‑1 (MQ=255)
                                                                                            aCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAAt                                                           <  1:4063983/150‑1 (MQ=255)
                                                                                            aCGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAAt                                                           <  1:1697709/150‑1 (MQ=255)
                                                                                              gAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATgg                                                         <  1:400941/150‑1 (MQ=255)
                                                                                               aaGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGt                                                        >  1:2413325/1‑150 (MQ=255)
                                                                                                   ggCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTgcgc                                                    >  2:3368789/1‑150 (MQ=255)
                                                                                                     cTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGt                                                  <  1:2189450/150‑1 (MQ=255)
                                                                                                       tttACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                >  1:4745836/1‑150 (MQ=255)
                                                                                                       tttACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGttt                                                <  1:2522722/150‑1 (MQ=255)
                                                                                                          aCCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                             >  2:2990635/1‑150 (MQ=255)
                                                                                                          aCCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                             <  1:4514491/150‑1 (MQ=255)
                                                                                                          aCCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTtctc                                             >  1:4966200/1‑150 (MQ=255)
                                                                                                             tGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTcgcg                                          <  2:4534968/150‑1 (MQ=255)
                                                                                                              ggCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCgg                                         <  2:3193151/150‑1 (MQ=255)
                                                                                                               gCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                        <  2:3426528/150‑1 (MQ=255)
                                                                                                               gCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGa                                        <  2:1108671/150‑1 (MQ=255)
                                                                                                                cAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAg                                       >  1:385759/1‑150 (MQ=255)
                                                                                                                     cggcgAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGcatca                                  <  1:2601884/150‑1 (MQ=255)
                                                                                                                           cccACCGCGACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTATCGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCt                            <  1:2303425/149‑1 (MQ=255)
                                                                                                                                  cTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGt                     <  2:3651513/150‑1 (MQ=255)
                                                                                                                                     cTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTa                  <  2:1840702/150‑1 (MQ=255)
                                                                                                                                      tCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTAt                 >  2:4239266/1‑150 (MQ=255)
                                                                                                                                       cGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCTCGGAGCATCATTTTCTTGTTGGTTTATc                >  2:1087490/1‑150 (MQ=255)
                                                                                                                                             ggCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGa          <  1:4126961/150‑1 (MQ=255)
                                                                                                                                                  ttttCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc     >  2:4677219/1‑150 (MQ=255)
                                                                                                                                                  ttttCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGcc     <  2:2322540/150‑1 (MQ=255)
                                                                                                                                                   ttACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGccc    >  1:2675252/1‑150 (MQ=255)
                                                                                                                                                     aCCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGc  <  2:3086969/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCCCTTTACCGCGAAAAGCCAGTCACGCATATGGACGATGACACGCTGGAAGCGTATGTCCGTCACTATATCGCTGCCAGCGAACCGCAAAACGAAGTGGCTTTTACCTGGCAGGGCGGCGAACCAACGCTACTCGGGCTGGCGTTTTACCGCCGTGCCGTAGCGCTACAGGCGAAATATGGTGCTGGCAGGAAGATAAGTAACAGCTTCCAGACTAACGGCGTGCTGCTGGATGACGAATGGTGCGCGTTTCTCGCGGAGCATCATTTTCTTGTTGGTTTATCGCTGGATGGCCCGC  >  CP009273/3976419‑3976717

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: