Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,576,325 T→A 100% N153Y (AAT→TAT)  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,576,3250TA98.1% 871.5 / 3.8 261N153Y (AAT→TAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base A (146/110);  total (149/112)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 264 total aligned reads displayed.

CGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCT  >  CP009273/1576180‑1576473
                                                                                                                                                 |                                                                                                                                                    
   gtttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                               >  2:2304713/1‑150 (MQ=255)
    tttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGa                                                                                                                                              >  1:3393701/1‑150 (MQ=255)
    tttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGACTAGa                                                                                                                                              >  2:2273651/1‑150 (MQ=255)
    tgtgtttATCTTCCGGTACCGGCACATGACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGa                                                                                                                                              <  2:5100948/148‑1 (MQ=255)
       gtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGacac                                                                                                                                           >  1:3681211/1‑150 (MQ=255)
             cTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgc                                                                                                                                     >  2:2490550/1‑150 (MQ=255)
              ttCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcg                                                                                                                                    >  1:2892371/1‑150 (MQ=255)
              ttCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcg                                                                                                                                    <  1:4826573/150‑1 (MQ=255)
               tCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgg                                                                                                                                   <  2:1427495/150‑1 (MQ=255)
                ccGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcggg                                                                                                                                  <  1:1686197/150‑1 (MQ=255)
                 cGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggc                                                                                                                                 >  2:1598647/1‑150 (MQ=255)
                 cGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggc                                                                                                                                 >  1:4619100/1‑150 (MQ=255)
                  ggTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcg                                                                                                                                <  2:3557944/150‑1 (MQ=255)
                   gTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcgg                                                                                                                               >  1:1964672/1‑150 (MQ=255)
                   gTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcgg                                                                                                                               <  2:3522255/150‑1 (MQ=255)
                    tACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGcgggcggg                                                                                                                              <  2:1287550/150‑1 (MQ=255)
                     aCCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGa                                                                                                                             <  1:1471660/150‑1 (MQ=255)
                      ccGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAg                                                                                                                            <  2:4324602/150‑1 (MQ=255)
                         gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTc                                                                                                                         <  1:4561267/150‑1 (MQ=255)
                             aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAcc                                                                                                                     >  1:913037/1‑150 (MQ=255)
                             aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAcc                                                                                                                     >  1:2584025/1‑150 (MQ=255)
                                  tGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAACGCGGGCGGGAGCTCGACCGGTCa                                                                                                                >  2:4385678/1‑150 (MQ=255)
                                  tGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCa                                                                                                                >  2:2544950/1‑150 (MQ=255)
                                    aTTTTGGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCatt                                                                                                              <  2:221524/150‑1 (MQ=255)
                                       tttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCattatt                                                                                                           <  2:1831914/150‑1 (MQ=255)
                                        ttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTg                                                                                                          <  1:4384197/150‑1 (MQ=255)
                                         tGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgc                                                                                                         >  2:2685094/1‑150 (MQ=255)
                                          gACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcg                                                                                                        <  1:4934481/150‑1 (MQ=255)
                                           aCAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcgg                                                                                                       <  1:1152039/150‑1 (MQ=255)
                                            cAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggc                                                                                                      >  2:345100/1‑150 (MQ=255)
                                             aaGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggcg                                                                                                     >  1:5029638/1‑150 (MQ=255)
                                             aaGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggcg                                                                                                     >  1:444524/1‑150 (MQ=255)
                                              aGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGgcgc                                                                                                    >  2:3831488/1‑150 (MQ=255)
                                                 gCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCggg                                                                                                 <  1:586677/150‑1 (MQ=255)
                                                  ccATTTACCTACTGCTGCAGTGTAATAACAATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAATCGGGCGGTAGCTCGACCGGTCATTATTGCGGCGCGGGa                                                                                                >  2:573047/1‑150 (MQ=255)
                                                  ccATTGACCTACTGCTGCAGTGTAAGAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGa                                                                                                <  2:4242428/150‑1 (MQ=255)
                                                   cATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                               <  1:2216524/150‑1 (MQ=255)
                                                       tACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                           <  1:4758523/150‑1 (MQ=255)
                                                       tACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                           >  2:4036608/1‑150 (MQ=255)
                                                         ccTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                         <  1:4000410/150‑1 (MQ=255)
                                                         ccTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                         <  2:1266904/150‑1 (MQ=255)
                                                          cTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCgg                                                                                        >  2:4963064/1‑150 (MQ=255)
                                                           tACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGa                                                                                       <  1:3735810/150‑1 (MQ=255)
                                                             ctgctgCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaaa                                                                                     <  2:4186422/150‑1 (MQ=255)
                                                               gctgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAaca                                                                                   >  2:384039/1‑150 (MQ=255)
                                                                 tgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACAcc                                                                                 >  1:1080291/1‑150 (MQ=255)
                                                                 tgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACAcc                                                                                 >  1:2393818/1‑150 (MQ=255)
                                                                    aGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtg                                                                              >  1:4438253/1‑150 (MQ=255)
                                                                       gtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                           >  2:643751/1‑150 (MQ=255)
                                                                         aataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCca                                                                         >  1:1631570/1‑150 (MQ=255)
                                                                          ataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcac                                                                        <  2:1764048/150‑1 (MQ=255)
                                                                           taaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcaca                                                                       <  1:1993864/150‑1 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGc                                                                   <  2:2981101/150‑1 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGc                                                                   <  2:3320495/150‑1 (MQ=255)
                                                                                aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGcc                                                                  <  2:1486115/150‑1 (MQ=255)
                                                                                 tGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCg                                                                 >  1:2528258/1‑150 (MQ=255)
                                                                                  gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGt                                                                >  2:2282480/1‑150 (MQ=255)
                                                                                  gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGt                                                                <  2:2237207/150‑1 (MQ=255)
                                                                                        ggAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaa                                                          >  1:1346099/1‑150 (MQ=255)
                                                                                        ggAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaa                                                          <  2:1335642/150‑1 (MQ=255)
                                                                                            taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACg                                                      >  2:339138/1‑150 (MQ=255)
                                                                                             aatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGt                                                     >  2:1933964/1‑150 (MQ=255)
                                                                                             aatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGt                                                     >  2:1235508/1‑150 (MQ=255)
                                                                                                    gCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAATCGTACGCCtt                                              >  2:4540803/1‑150 (MQ=255)
                                                                                                     cAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTc                                             >  1:2407698/1‑150 (MQ=255)
                                                                                                      aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCa                                            >  2:1592601/1‑150 (MQ=255)
                                                                                                      aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCa                                            >  2:1459288/1‑150 (MQ=255)
                                                                                                          aaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa                                        >  2:4578912/1‑150 (MQ=255)
                                                                                                           aGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                       >  1:2063679/1‑150 (MQ=255)
                                                                                                                   ttAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaaa                               <  1:3894480/150‑1 (MQ=255)
                                                                                                                   ttAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaaa                               >  2:2394094/1‑150 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                              >  1:2025546/1‑150 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                              >  2:3685204/1‑150 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                              <  1:4162043/150‑1 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                              >  1:4332543/1‑150 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                              <  1:1584157/150‑1 (MQ=255)
                                                                                                                     aGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAgg                             <  2:1161980/150‑1 (MQ=255)
                                                                                                                        ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                          >  2:2498282/1‑150 (MQ=255)
                                                                                                                        ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                          <  1:1051958/150‑1 (MQ=255)
                                                                                                                           aTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                       <  1:2814279/150‑1 (MQ=255)
                                                                                                                            tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCg                      <  2:2594530/150‑1 (MQ=255)
                                                                                                                               cATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCg                   >  2:749632/1‑150 (MQ=255)
                                                                                                                                 tCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGtt                 <  2:1663753/150‑1 (MQ=255)
                                                                                                                                 tCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGtt                 <  1:1293952/150‑1 (MQ=255)
                                                                                                                                  ccTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg                <  2:3471906/150‑1 (MQ=255)
                                                                                                                                  ccTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg                >  1:3075018/1‑150 (MQ=255)
                                                                                                                                  ccTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg                >  1:4611009/1‑150 (MQ=255)
                                                                                                                                   cTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa               >  2:330977/1‑150 (MQ=255)
                                                                                                                                   cTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa               >  1:3825604/1‑150 (MQ=255)
                                                                                                                                      aGCATCGGTATAGGAATAGACACCAACGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt            <  2:5125075/150‑1 (MQ=255)
                                                                                                                                      aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt            >  1:1891735/1‑150 (MQ=255)
                                                                                                                                      aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt            <  1:693568/150‑1 (MQ=255)
                                                                                                                                      aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt            >  1:114803/1‑150 (MQ=255)
                                                                                                                                        cATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttttt          <  2:3183091/150‑1 (MQ=255)
                                                                                                                                         aTCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg         >  2:901563/1‑150 (MQ=255)
                                                                                                                                          tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtgt        <  2:4582247/150‑1 (MQ=255)
                                                                                                                                          tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtgt        <  1:3049563/150‑1 (MQ=255)
                                                                                                                                            ggTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc      >  1:853562/1‑150 (MQ=255)
                                                                                                                                            ggTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc      >  2:1053610/1‑150 (MQ=255)
                                                                                                                                             gTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGct     >  2:3005113/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
CGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAGCT  >  CP009273/1576180‑1576473

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: