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,3250TA97.2% 577.6 / 4.4 180N153Y (AAT→TAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (3/2);  new base A (101/74);  total (104/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 180 total aligned reads displayed.

CGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576180‑1576469
                                                                                                                                                 |                                                                                                                                                
 gcgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGa                                                                                                                                               <  1:1681567/148‑1 (MQ=255)
  cgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTCTAGGaa                                                                                                                                              >  2:681035/1‑148 (MQ=255)
     ttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                           >  1:2555721/1‑148 (MQ=255)
     ttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                           >  1:2095299/1‑148 (MQ=255)
      tgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCACGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTCTAGGAATAGa                                                                                                                                          >  1:3073028/1‑148 (MQ=255)
           aTCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCa                                                                                                                                     >  2:1277816/1‑148 (MQ=255)
            tCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaa                                                                                                                                    <  1:3419509/148‑1 (MQ=255)
             cTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaaa                                                                                                                                   <  2:1373451/148‑1 (MQ=255)
             cTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaaa                                                                                                                                   >  2:2564008/1‑148 (MQ=255)
              ttCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAg                                                                                                                                  <  1:2813683/148‑1 (MQ=255)
                ccGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcg                                                                                                                                <  2:425599/148‑1 (MQ=255)
                  gggaCCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAgcggg                                                                                                                              <  2:3073134/145‑1 (MQ=255)
                         gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       >  1:1293942/1‑148 (MQ=255)
                         gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       <  2:889753/148‑1 (MQ=255)
                          cacaTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                      >  1:2934349/1‑148 (MQ=255)
                             aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGa                                                                                                                   <  2:3089375/148‑1 (MQ=255)
                              ttACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAc                                                                                                                  <  1:3079172/148‑1 (MQ=255)
                                 cTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCgg                                                                                                               <  1:60489/148‑1 (MQ=255)
                                  tGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGACTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGt                                                                                                              >  2:3504716/1‑148 (MQ=255)
                                    aTTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCa                                                                                                            <  2:2290911/148‑1 (MQ=255)
                                        ttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCattat                                                                                                        >  1:2765681/1‑148 (MQ=255)
                                           aCAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgc                                                                                                     <  2:388402/148‑1 (MQ=255)
                                             aaGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcgg                                                                                                   <  1:316418/148‑1 (MQ=255)
                                              aGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggc                                                                                                  >  1:690992/1‑148 (MQ=255)
                                                 gCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGcgcg                                                                                               >  2:2739111/1‑148 (MQ=255)
                                                 gCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGcgcg                                                                                               >  1:2989164/1‑148 (MQ=255)
                                                  ccATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCgg                                                                                              <  2:3372164/148‑1 (MQ=255)
                                                   cATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTCCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCggg                                                                                             >  2:1725351/1‑148 (MQ=255)
                                                      ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCGTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  2:3432969/1‑148 (MQ=255)
                                                      ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  2:860867/1‑148 (MQ=255)
                                                         ccTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                       <  1:3181794/148‑1 (MQ=255)
                                                               gctgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaaa                                                                                 >  1:479262/1‑148 (MQ=255)
                                                                   cAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACAcc                                                                             >  1:1011987/1‑148 (MQ=255)
                                                                      tgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtg                                                                          <  1:2433471/148‑1 (MQ=255)
                                                                         aataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                       >  2:3076491/1‑148 (MQ=255)
                                                                         aataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                       >  2:2650974/1‑148 (MQ=255)
                                                                         aataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                       >  2:1555702/1‑148 (MQ=255)
                                                                              ccATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACAt                                                                  <  1:1739754/148‑1 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATa                                                                 >  1:3363752/1‑148 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATa                                                                 >  1:2303472/1‑148 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATa                                                                 >  2:3191173/1‑148 (MQ=255)
                                                                               cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATa                                                                 >  2:3191174/1‑148 (MQ=255)
                                                                                aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTATTAGGTTAATTCCATCCTGAGCATCGGTATATTAATAGACACCAAATCGGGCGGGAGCTCGTCCGGTCATTATTGAGGGGCGGTAGGGGCCGGAAACTCCGTGTGCCAAATAg                                                                <  2:2999869/148‑1 (MQ=255)
                                                                                 tGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGc                                                               <  2:3157446/148‑1 (MQ=255)
                                                                                  gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGcc                                                              >  1:1768311/1‑148 (MQ=255)
                                                                                      cTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTa                                                          <  1:2646536/148‑1 (MQ=255)
                                                                                       tGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAg                                                         >  1:890543/1‑148 (MQ=255)
                                                                                        ggAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGt                                                        >  2:2127614/1‑148 (MQ=255)
                                                                                        ggAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGt                                                        >  2:953267/1‑148 (MQ=255)
                                                                                         gAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTa                                                       <  1:825117/148‑1 (MQ=255)
                                                                                         gAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTa                                                       <  1:3342498/148‑1 (MQ=255)
                                                                                         gAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTa                                                       >  2:2373661/1‑148 (MQ=255)
                                                                                         gAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTa                                                       >  1:454748/1‑148 (MQ=255)
                                                                                           ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaa                                                     <  1:625860/148‑1 (MQ=255)
                                                                                            taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    >  1:874784/1‑148 (MQ=255)
                                                                                              atTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACg                                                  >  1:413620/1‑148 (MQ=255)
                                                                                              atTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACg                                                  <  2:3167090/148‑1 (MQ=255)
                                                                                              atTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATAGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACg                                                  >  2:1673986/1‑148 (MQ=255)
                                                                                                tCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                >  2:2859938/1‑148 (MQ=255)
                                                                                                  aGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACg                                              >  1:592480/1‑148 (MQ=255)
                                                                                                      aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGNCCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCtt                                          <  1:1808924/148‑1 (MQ=255)
                                                                                                       aGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTc                                         <  2:2832045/148‑1 (MQ=255)
                                                                                                       aGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTc                                         <  2:1828481/148‑1 (MQ=255)
                                                                                                         aaaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCAt                                       >  1:1412454/1‑148 (MQ=255)
                                                                                                           aGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATATCCGTTAGTAAAACGTACGCCTTCATcc                                     >  2:1269673/1‑148 (MQ=255)
                                                                                                             tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   >  2:3459589/1‑148 (MQ=255)
                                                                                                             tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAt                                   >  1:1431863/1‑148 (MQ=255)
                                                                                                              ttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAtt                                  >  2:3299111/1‑148 (MQ=255)
                                                                                                              ttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAtt                                  >  1:1150688/1‑148 (MQ=255)
                                                                                                                cTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGCGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTaa                                >  2:1889655/1‑148 (MQ=255)
                                                                                                                 tGTTAGCGGAATACCATCCTGAGCATCGGTATAGTAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAAt                               <  1:1521459/148‑1 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                            <  2:2989164/148‑1 (MQ=255)
                                                                                                                    tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                            >  1:3465720/1‑148 (MQ=255)
                                                                                                                     aGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaaa                           <  1:1007894/148‑1 (MQ=255)
                                                                                                                      gCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                          >  2:2120400/1‑148 (MQ=255)
                                                                                                                      gCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                          >  2:2049484/1‑148 (MQ=255)
                                                                                                                       cGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAgg                         >  1:3111477/1‑148 (MQ=255)
                                                                                                                        ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGa                        <  1:2650978/148‑1 (MQ=255)
                                                                                                                         gAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAg                       <  2:820596/148‑1 (MQ=255)
                                                                                                                          aaTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      <  2:1197065/148‑1 (MQ=255)
                                                                                                                          aaTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  2:2938665/1‑148 (MQ=255)
                                                                                                                          aaTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  2:2519624/1‑148 (MQ=255)
                                                                                                                          aaTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  2:1847514/1‑148 (MQ=255)
                                                                                                                            tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGt                    <  1:463510/148‑1 (MQ=255)
                                                                                                                             aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   <  1:860867/148‑1 (MQ=255)
                                                                                                                             aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   >  1:1664141/1‑148 (MQ=255)
                                                                                                                             aCCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                   >  2:2700615/1‑148 (MQ=255)
                                                                                                                                  ccTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGt              <  2:661819/148‑1 (MQ=255)
                                                                                                                                   cTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGtt             >  1:1636780/1‑148 (MQ=255)
                                                                                                                                    tGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg            <  1:1211286/148‑1 (MQ=255)
                                                                                                                                     gAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  1:2562099/1‑148 (MQ=255)
                                                                                                                                     gAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  1:1363373/1‑148 (MQ=255)
                                                                                                                                     gAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  1:1363366/1‑148 (MQ=255)
                                                                                                                                      aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt          <  1:3162991/148‑1 (MQ=255)
                                                                                                                                        cATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        >  1:511269/1‑148 (MQ=255)
                                                                                                                                        cATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGAACGCCTTCATCCAGAAATCAAAGGAGAGTCGGCGtggtatg        >  1:908730/1‑142 (MQ=255)
                                                                                                                                          tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttttt      >  2:104080/1‑148 (MQ=255)
                                                                                                                                           cGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg     >  1:2221012/1‑148 (MQ=255)
                                                                                                                                           cGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg     >  1:3315549/1‑148 (MQ=255)
                                                                                                                                             gTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtgtg   <  2:2681734/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576180‑1576469

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: ‑