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,3250TA96.1% 808.2 / 16.6 253N153Y (AAT→TAT) ydeNputative Ser‑type periplasmic non‑aryl sulfatase
Reads supporting (aligned to +/- strand):  ref base T (2/8);  new base A (137/106);  total (139/114)
Fisher's exact test for biased strand distribution p-value = 4.63e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 253 total aligned reads displayed.

ACGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576179‑1576469
                                                                                                                                                  |                                                                                                                                                
aCGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAgg                                                                                                                                                >  2:3743300/1‑149 (MQ=255)
   cgTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAAt                                                                                                                                             <  1:2621491/149‑1 (MQ=255)
       tgtttATCTTCCGGGACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGac                                                                                                                                         <  1:2704102/149‑1 (MQ=255)
         tttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGacac                                                                                                                                       >  2:3673306/1‑149 (MQ=255)
          ttATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACAcc                                                                                                                                      <  1:2252683/149‑1 (MQ=255)
           tATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCa                                                                                                                                     >  2:1663522/1‑149 (MQ=255)
            aTCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCaa                                                                                                                                    <  2:2518356/149‑1 (MQ=255)
               ttCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgc                                                                                                                                 >  1:1207620/1‑149 (MQ=255)
                   ggTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggc                                                                                                                             <  1:1678871/149‑1 (MQ=255)
                    gTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcg                                                                                                                            >  2:3163495/1‑149 (MQ=255)
                     tACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAgcgggcgg                                                                                                                           <  1:3664308/149‑1 (MQ=255)
                       ccGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGa                                                                                                                         <  2:2762047/149‑1 (MQ=255)
                       ccGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGa                                                                                                                         <  2:566133/149‑1 (MQ=255)
                         ggCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       <  2:2800709/149‑1 (MQ=255)
                         ggCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGc                                                                                                                       <  1:72358/149‑1 (MQ=255)
                          gCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATTATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCt                                                                                                                      <  2:2678274/149‑1 (MQ=255)
                            acaTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCg                                                                                                                    <  2:436913/149‑1 (MQ=255)
                             caTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGa                                                                                                                   >  2:1243558/1‑149 (MQ=255)
                             caTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGa                                                                                                                   >  1:3559477/1‑149 (MQ=255)
                              aTTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGAc                                                                                                                  >  2:1750411/1‑149 (MQ=255)
                                   tGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTc                                                                                                             >  2:3965951/1‑149 (MQ=255)
                                    gATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCa                                                                                                            >  2:2228013/1‑149 (MQ=255)
                                         ttGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCattatt                                                                                                       <  2:896328/149‑1 (MQ=255)
                                           gACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgc                                                                                                     >  1:1902381/1‑149 (MQ=255)
                                             cAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCACGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgtgg                                                                                                   >  2:1588705/1‑149 (MQ=255)
                                             cAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcgg                                                                                                   >  2:586176/1‑149 (MQ=255)
                                              aaGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTgcggc                                                                                                  >  1:4059657/1‑149 (MQ=255)
                                                 tGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGcgcg                                                                                               <  2:2977760/149‑1 (MQ=255)
                                                 tGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGcgcg                                                                                               <  2:279744/149‑1 (MQ=255)
                                                     aTTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAg                                                                                           >  1:92928/1‑149 (MQ=255)
                                                      tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  1:3206790/1‑149 (MQ=255)
                                                      tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  1:789961/1‑149 (MQ=255)
                                                      tttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgg                                                                                          >  1:787032/1‑149 (MQ=255)
                                                       ttACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAggg                                                                                         <  1:2958328/149‑1 (MQ=255)
                                                        tACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgggg                                                                                        >  2:3950229/1‑149 (MQ=255)
                                                        tACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgggg                                                                                        >  2:3949547/1‑149 (MQ=255)
                                                        tACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAgggg                                                                                        >  2:3274202/1‑149 (MQ=255)
                                                         aCCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                       <  2:3908565/149‑1 (MQ=255)
                                                         aCCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                       >  2:1931821/1‑149 (MQ=255)
                                                         aCCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGc                                                                                       >  2:1932223/1‑149 (MQ=255)
                                                           cTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCg                                                                                     >  1:1755948/1‑149 (MQ=255)
                                                             aCTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGa                                                                                   >  2:720214/1‑149 (MQ=255)
                                                              ctgctgCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaa                                                                                  >  1:141590/1‑149 (MQ=255)
                                                              ctgctgCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGaa                                                                                  >  2:1088615/1‑149 (MQ=255)
                                                                      gtgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtg                                                                          >  2:3494792/1‑149 (MQ=255)
                                                                       tgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtgt                                                                         >  2:3065337/1‑149 (MQ=255)
                                                                       tgtAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCgtgt                                                                         >  2:2256676/1‑149 (MQ=255)
                                                                         taataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                       >  1:1371112/1‑149 (MQ=255)
                                                                         taataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGc                                                                       >  1:1370887/1‑149 (MQ=255)
                                                                          aataaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGcc                                                                      >  1:97298/1‑149 (MQ=255)
                                                                            taaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcac                                                                    <  1:1732201/149‑1 (MQ=255)
                                                                             aaCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCcaca                                                                   <  1:3965951/149‑1 (MQ=255)
                                                                              aCCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACAt                                                                  <  1:1534905/149‑1 (MQ=255)
                                                                                 aTGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGc                                                               >  1:2130527/1‑149 (MQ=255)
                                                                                  tGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGcc                                                              >  2:2140681/1‑149 (MQ=255)
                                                                                    aTTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGt                                                            <  1:3626925/149‑1 (MQ=255)
                                                                                     ttCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGtt                                                           >  2:420889/1‑149 (MQ=255)
                                                                                     ttCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGtt                                                           >  2:2448775/1‑149 (MQ=255)
                                                                                       cTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAg                                                         >  1:2051415/1‑149 (MQ=255)
                                                                                       cTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAg                                                         <  2:2782104/149‑1 (MQ=255)
                                                                                            ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    <  2:3033870/149‑1 (MQ=255)
                                                                                            ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    <  2:773960/149‑1 (MQ=255)
                                                                                            ataatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTaaaa                                                    <  1:1425662/149‑1 (MQ=255)
                                                                                             taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                   <  2:257087/149‑1 (MQ=255)
                                                                                             taatTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAAc                                                   <  2:2875147/149‑1 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                >  1:3751213/1‑149 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                <  2:749219/149‑1 (MQ=255)
                                                                                                ttCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTa                                                <  1:1243558/149‑1 (MQ=255)
                                                                                                 tCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTAc                                               <  2:1346273/149‑1 (MQ=255)
                                                                                                   aGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGc                                             >  1:3277999/1‑149 (MQ=255)
                                                                                                    ggCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGcc                                            <  2:819115/149‑1 (MQ=255)
                                                                                                    ggCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGcc                                            >  2:4009441/1‑149 (MQ=255)
                                                                                                       aaGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTc                                         <  1:2524448/149‑1 (MQ=255)
                                                                                                           aaGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATcc                                     >  2:619937/1‑149 (MQ=255)
                                                                                                            aGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCa                                    <  1:4005991/149‑1 (MQ=255)
                                                                                                              tttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCAtt                                  >  2:1809366/1‑149 (MQ=255)
                                                                                                               ttCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTa                                 >  2:3535353/1‑149 (MQ=255)
                                                                                                                tCTGTTAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTaa                                <  1:3132422/149‑1 (MQ=255)
                                                                                                                    ttAGCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGaa                            >  2:7905/1‑149 (MQ=255)
                                                                                                                       gCGGAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAgg                         <  2:1118148/149‑1 (MQ=255)
                                                                                                                        ctgAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGa                        <  1:539395/147‑1 (MQ=255)
                                                                                                                         ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAg                       <  1:2802755/149‑1 (MQ=255)
                                                                                                                         ggAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAg                       <  2:1619568/149‑1 (MQ=255)
                                                                                                                          gAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  1:1541668/1‑149 (MQ=255)
                                                                                                                          gAATACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGc                      >  1:2576521/1‑149 (MQ=255)
                                                                                                                            aTACCATCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGt                    >  1:2223705/1‑149 (MQ=255)
                                                                                                                                  tCCTGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGt              <  1:44069/149‑1 (MQ=255)
                                                                                                                                     tGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  2:1490314/1‑149 (MQ=255)
                                                                                                                                     tGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           <  2:1987167/149‑1 (MQ=255)
                                                                                                                                     tGAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGa           >  2:2645504/1‑149 (MQ=255)
                                                                                                                                      gAGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAt          <  1:4026341/149‑1 (MQ=255)
                                                                                                                                       aGCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAtt         <  2:2588440/149‑1 (MQ=255)
                                                                                                                                        gCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        <  1:1088615/149‑1 (MQ=255)
                                                                                                                                        gCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        <  1:3352027/149‑1 (MQ=255)
                                                                                                                                        gCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        >  1:3054311/1‑149 (MQ=255)
                                                                                                                                        gCATCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGAttt        <  1:1168047/149‑1 (MQ=255)
                                                                                                                                           tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGTTTTTtg     >  2:3051347/1‑149 (MQ=255)
                                                                                                                                           tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg     >  2:2081187/1‑149 (MQ=255)
                                                                                                                                           tCGGTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTtg     <  1:1820666/149‑1 (MQ=255)
                                                                                                                                              gTATAGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGc  >  2:3817234/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACGCGTTTGTTTATCTTCCGGTACCGGCACATTACTGATTTTTGACAAGTGCCATTTACCTACTGCTGCAGTGTAATAACCATGATTCTGGAATAATTCAGGCAAGAAAGTTTCTGTTAGCGGAATACCATCCTGAGCATCGGTATTGGAATAGACACCAAAGCGGGCGGGAGCTCGACCGGTCATTATTGCGGCGCGGGAGGGGCCGGAAACACCGTGTGCCACATAGCCGTTAGTAAAACGTACGCCTTCATCCATTAATGAAAGGAGCGTCGGCGTTGATTTTTGTGC  >  CP009273/1576179‑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: ‑