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.6% 972.1 / 1.0 297N153Y (AAT→TAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base A (167/125);  total (169/128)
Fisher's exact test for biased strand distribution p-value = 6.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 298 total aligned reads displayed.

ACGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  CP009273/1576179‑1576471
                                                                                                                                                  |                                                                                                                                                  
aCGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGa                                                                                                                                                 <  2:880304/150‑1 (MQ=255)
  gcgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAAt                                                                                                                                               <  2:822884/150‑1 (MQ=255)
   cgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATa                                                                                                                                              >  1:1497336/1‑150 (MQ=255)
    gtttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                             <  2:941112/150‑1 (MQ=255)
    gtttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAg                                                                                                                                             >  1:907875/1‑150 (MQ=255)
     tttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGa                                                                                                                                            >  1:446472/1‑150 (MQ=255)
      ttgtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGac                                                                                                                                           <  1:2041494/150‑1 (MQ=255)
        gtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGacac                                                                                                                                         >  1:4616724/1‑150 (MQ=255)
        gtttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGacac                                                                                                                                         >  2:3854823/1‑150 (MQ=255)
         tttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACAcc                                                                                                                                        <  2:1058058/150‑1 (MQ=255)
              cTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgc                                                                                                                                   >  2:2734936/1‑150 (MQ=255)
              cTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgc                                                                                                                                   >  2:4715575/1‑150 (MQ=255)
               ttCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAgcg                                                                                                                                  <  2:702226/150‑1 (MQ=255)
                   ggTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcg                                                                                                                              >  2:1430882/1‑150 (MQ=255)
                      aCCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGa                                                                                                                           >  2:1954952/1‑150 (MQ=255)
                         ggCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                        >  2:5203908/1‑150 (MQ=255)
                         ggCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                        <  1:5267675/150‑1 (MQ=255)
                         ggCACATTACAGATTTATGACAAGTGCCAATTACCTACTGCAGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                        >  1:1700048/1‑150 (MQ=255)
                               ttACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCg                                                                                                                  >  2:4267597/1‑150 (MQ=255)
                               ttACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCg                                                                                                                  >  1:3037633/1‑150 (MQ=255)
                                 aCTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGt                                                                                                                >  2:4403062/1‑150 (MQ=255)
                                  cTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTc                                                                                                               >  1:5237858/1‑150 (MQ=255)
                                  cTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTAGAGGAATAGACACCAAATCGGGCGGGAGCTCGACCGGTc                                                                                                               >  2:3582275/1‑150 (MQ=255)
                                   tGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCa                                                                                                              >  2:4991943/1‑150 (MQ=255)
                                    gATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCat                                                                                                             >  2:3792702/1‑150 (MQ=255)
                                     aTTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGACTAGACACCAAAGCGGGCGGGAGCTCGACCGGTCatt                                                                                                            >  2:4957772/1‑150 (MQ=255)
                                          tGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgc                                                                                                       <  2:1741174/150‑1 (MQ=255)
                                           gACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcg                                                                                                      <  2:550765/150‑1 (MQ=255)
                                                 tGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCgg                                                                                                >  1:4572790/1‑150 (MQ=255)
                                                 tGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCgg                                                                                                >  1:1756463/1‑150 (MQ=255)
                                                    cATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                             >  1:2615232/1‑150 (MQ=255)
                                                     aTTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                            >  2:5187603/1‑150 (MQ=255)
                                                     aTTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                            >  2:219146/1‑150 (MQ=255)
                                                     aTTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                            >  2:4080404/1‑150 (MQ=255)
                                                       ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTGAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgggg                                                                                          <  2:3934533/150‑1 (MQ=255)
                                                          ccTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                       >  2:238673/1‑150 (MQ=255)
                                                           cTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCgg                                                                                      >  2:2345800/1‑150 (MQ=255)
                                                            tACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGa                                                                                     >  2:3083908/1‑150 (MQ=255)
                                                            tACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGa                                                                                     <  1:3555876/150‑1 (MQ=255)
                                                                  tgcAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACAcc                                                                               <  2:1775138/150‑1 (MQ=255)
                                                                    cAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTGGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgt                                                                             <  2:1512347/150‑1 (MQ=255)
                                                                    cAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgt                                                                             >  2:4787954/1‑150 (MQ=255)
                                                                     aGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtg                                                                            >  1:3976291/1‑150 (MQ=255)
                                                                      gtgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtgt                                                                           >  2:3001510/1‑150 (MQ=255)
                                                                       tgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGtgtg                                                                          <  1:2852286/150‑1 (MQ=255)
                                                                        gtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                         <  1:1170555/150‑1 (MQ=255)
                                                                         taataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGcc                                                                        >  1:705573/1‑150 (MQ=255)
                                                                           ataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcac                                                                      >  2:2893421/1‑150 (MQ=255)
                                                                            taaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcaca                                                                     >  1:4876740/1‑150 (MQ=255)
                                                                               ccATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAg                                                                  >  1:4811777/1‑150 (MQ=255)
                                                                               ccATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAg                                                                  <  2:4275035/150‑1 (MQ=255)
                                                                               ccATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAg                                                                  >  1:939393/1‑150 (MQ=255)
                                                                                cATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGc                                                                 >  1:4372945/1‑150 (MQ=255)
                                                                                 aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGcc                                                                >  2:58899/1‑150 (MQ=255)
                                                                                   gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGt                                                              <  2:1658187/150‑1 (MQ=255)
                                                                                   gATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGt                                                              >  2:3905016/1‑150 (MQ=255)
                                                                                     ttCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTa                                                            <  2:4547855/150‑1 (MQ=255)
                                                                                      tCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGAGTGCCACATAGCCGTTAg                                                           >  2:1191865/1‑150 (MQ=255)
                                                                                      tCTGGAATAATTCAGGCAAGAAAGTGTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAg                                                           <  2:4015687/150‑1 (MQ=255)
                                                                                         ggAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaa                                                        >  2:2297701/1‑150 (MQ=255)
                                                                                            ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                     >  2:1492197/1‑150 (MQ=255)
                                                                                             taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACg                                                    >  2:3546345/1‑150 (MQ=255)
                                                                                              aatTCAGGCAAGAAAGTTTCTGTTAGCGGATTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGt                                                   <  2:3362548/150‑1 (MQ=255)
                                                                                               atTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                  <  1:2370050/150‑1 (MQ=255)
                                                                                               atTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                  >  2:4325773/1‑150 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                 <  1:2812114/150‑1 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                 <  1:1430882/150‑1 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                                 <  1:558181/150‑1 (MQ=255)
                                                                                                    ggCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCt                                             >  2:2124072/1‑150 (MQ=255)
                                                                                                    ggCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCt                                             >  1:4165765/1‑150 (MQ=255)
                                                                                                      catgAAAGTTTCTGTTAGCTGTTGACGATCCTGAGCATCGGTATTGGTATAGACACGATAGCGGGCGGGAGCTCGTCCGGTCATTTTTGGGGCGCGGGAGGGGGCTGAATCTCCGTGTGGCACATAGCCGTTAGTTAAACGTACGCCTTc                                           <  2:4690563/147‑1 (MQ=255)
                                                                                                       aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCa                                          >  2:5034428/1‑150 (MQ=255)
                                                                                                          aaaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATcc                                       <  1:3835357/150‑1 (MQ=255)
                                                                                                          aaaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATcc                                       >  2:5322159/1‑150 (MQ=255)
                                                                                                          aaaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATcc                                       <  1:3292687/150‑1 (MQ=255)
                                                                                                           aaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa                                      >  1:1830121/1‑150 (MQ=255)
                                                                                                           aaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa                                      >  1:1233769/1‑150 (MQ=255)
                                                                                                                   gTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAATGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                              >  2:1792832/1‑150 (MQ=255)
                                                                                                                   gTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                              >  2:2464873/1‑150 (MQ=255)
                                                                                                                     tAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAg                            <  1:1258379/150‑1 (MQ=255)
                                                                                                                      aGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAgg                           >  2:1563152/1‑150 (MQ=255)
                                                                                                                      aGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAgg                           >  2:2824536/1‑150 (MQ=255)
                                                                                                                         ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                        <  1:1011268/150‑1 (MQ=255)
                                                                                                                         ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                        <  1:5329409/150‑1 (MQ=255)
                                                                                                                         ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                        <  2:3714021/150‑1 (MQ=255)
                                                                                                                          gAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCg                       >  2:1339573/1‑150 (MQ=255)
                                                                                                                            aTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTc                     <  1:3084820/150‑1 (MQ=255)
                                                                                                                             tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGTCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCg                    <  2:265445/150‑1 (MQ=255)
                                                                                                                             tACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCg                    <  2:2227362/150‑1 (MQ=255)
                                                                                                                                cATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCg                 >  2:2556568/1‑150 (MQ=255)
                                                                                                                                   ccTGAGCATCGGTATAGGAATAGACACCACAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg              <  2:3460946/150‑1 (MQ=255)
                                                                                                                                   ccTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTg              <  2:4116097/150‑1 (MQ=255)
                                                                                                                                    cTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa             >  2:4234987/1‑150 (MQ=255)
                                                                                                                                     tGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt            >  1:367245/1‑150 (MQ=255)
                                                                                                                                      gAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAtt           <  2:5102312/150‑1 (MQ=255)
                                                                                                                                       aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt          >  2:31902/1‑150 (MQ=255)
                                                                                                                                          aTCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg       <  2:3670730/150‑1 (MQ=255)
                                                                                                                                          aTCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg       <  2:1880759/150‑1 (MQ=255)
                                                                                                                                             ggTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc    <  2:1453998/150‑1 (MQ=255)
                                                                                                                                              gTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCa   >  1:3455410/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
ACGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGCAG  >  CP009273/1576179‑1576471

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