Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,038,158 C→A 29.3% A503A (GCG→GCT)  etk ← tyrosine‑protein kinase, role in O‑antigen capsule formation

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,038,1580CA29.3% 448.8 / 156.6 224A503A (GCG→GCT) etktyrosine‑protein kinase, role in O‑antigen capsule formation
Reads supporting (aligned to +/- strand):  ref base C (68/90);  new base A (25/41);  total (93/131)
Fisher's exact test for biased strand distribution p-value = 5.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

GAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCT  >  CP009273/1038013‑1038305
                                                                                                                                                 |                                                                                                                                                   
gagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                                                                  <  2:753338/149‑1 (MQ=255)
gagTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGcc                                                                                                                                                  >  2:3738252/1‑149 (MQ=255)
 agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCa                                                                                                                                                 <  2:1942420/149‑1 (MQ=255)
 agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCa                                                                                                                                                 >  1:3075836/1‑149 (MQ=255)
 agTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCa                                                                                                                                                 >  2:656602/1‑149 (MQ=255)
  gTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAg                                                                                                                                                >  2:518425/1‑149 (MQ=255)
    tGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGa                                                                                                                                              <  1:2384641/149‑1 (MQ=255)
       aCTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAgg                                                                                                                                           <  2:3140381/149‑1 (MQ=255)
           aCAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGat                                                                                                                                       <  1:2997679/149‑1 (MQ=255)
           aCAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGat                                                                                                                                       <  1:940438/149‑1 (MQ=255)
            cAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGata                                                                                                                                      <  2:1025766/149‑1 (MQ=255)
              aaCGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAtt                                                                                                                                    <  1:2752256/149‑1 (MQ=255)
                cGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAtttt                                                                                                                                  <  1:2924329/149‑1 (MQ=255)
                 gTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATAttttt                                                                                                                                 <  1:917433/149‑1 (MQ=255)
                  ttttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTa                                                                                                                                <  2:1347178/149‑1 (MQ=255)
                    ttACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGt                                                                                                                              <  2:1333013/149‑1 (MQ=255)
                      aCCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTAc                                                                                                                            <  1:1255947/149‑1 (MQ=255)
                        cACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGa                                                                                                                          <  1:1421301/149‑1 (MQ=255)
                        cACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGa                                                                                                                          <  1:710861/149‑1 (MQ=255)
                              cTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGgcgc                                                                                                                    <  1:4068615/149‑1 (MQ=255)
                               tGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgct                                                                                                                   <  1:2888123/149‑1 (MQ=255)
                               tGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCgct                                                                                                                   <  2:3043544/149‑1 (MQ=255)
                                 gCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                                 <  1:1099533/149‑1 (MQ=255)
                                 gCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                                 >  2:3352614/1‑149 (MQ=255)
                                 gCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgc                                                                                                                 <  2:471963/149‑1 (MQ=255)
                                  cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgct                                                                                                                <  1:1835653/149‑1 (MQ=255)
                                  cGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCgctgct                                                                                                                <  1:2747647/149‑1 (MQ=255)
                                   gTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                                                               <  1:3460236/149‑1 (MQ=255)
                                   gTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                                                               <  1:992070/149‑1 (MQ=255)
                                   gTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGctgctg                                                                                                               >  2:834040/1‑149 (MQ=255)
                                    tCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGa                                                                                                              <  2:848153/149‑1 (MQ=255)
                                     cGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAt                                                                                                             <  2:2824624/149‑1 (MQ=255)
                                      gCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGAtt                                                                                                            >  2:3224608/1‑149 (MQ=255)
                                          cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaa                                                                                                        <  1:2165678/149‑1 (MQ=255)
                                          cGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaa                                                                                                        <  1:2545540/149‑1 (MQ=255)
                                           ggTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaa                                                                                                       <  2:706524/149‑1 (MQ=255)
                                            gTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCAAGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGaaaa                                                                                                      <  2:32780/149‑1 (MQ=255)
                                              gATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaata                                                                                                    <  2:2221642/149‑1 (MQ=255)
                                              gATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaata                                                                                                    >  2:3869810/1‑149 (MQ=255)
                                               aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataa                                                                                                   >  2:3376132/1‑149 (MQ=255)
                                               aTCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataa                                                                                                   <  2:3516339/149‑1 (MQ=255)
                                                tCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAaaataaa                                                                                                  >  2:3598091/1‑149 (MQ=255)
                                                  aTTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtt                                                                                                <  2:2008797/149‑1 (MQ=255)
                                                   ttAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               >  2:3770674/1‑149 (MQ=255)
                                                   ttAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               >  2:1727222/1‑149 (MQ=255)
                                                   ttAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAttt                                                                                               <  2:2325414/149‑1 (MQ=255)
                                                    tAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAAtttt                                                                                              <  1:1428951/149‑1 (MQ=255)
                                                      gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTc                                                                                            >  1:3156690/1‑149 (MQ=255)
                                                      gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTc                                                                                            >  1:823233/1‑149 (MQ=255)
                                                      gAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTc                                                                                            <  1:2766183/149‑1 (MQ=255)
                                                        aTGTTAGTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCtt                                                                                          <  1:3097509/149‑1 (MQ=255)
                                                         tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                         >  1:1949264/1‑149 (MQ=255)
                                                         tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                         <  1:2119440/149‑1 (MQ=255)
                                                         tGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTa                                                                                         <  1:449820/149‑1 (MQ=255)
                                                          gTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTAc                                                                                        <  1:3782629/149‑1 (MQ=255)
                                                           ttATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACg                                                                                       >  1:769000/1‑149 (MQ=255)
                                                            tATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGc                                                                                      <  2:156616/149‑1 (MQ=255)
                                                             aTTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCa                                                                                     <  1:2145785/149‑1 (MQ=255)
                                                              ttCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAg                                                                                    >  1:2652271/1‑149 (MQ=255)
                                                                ctcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGAc                                                                                  >  1:3268100/1‑149 (MQ=255)
                                                                 tcCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTAc                                                                                        >  2:3189730/1‑142 (MQ=255)
                                                                   cGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACATCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACggg                                                                               >  1:1154305/1‑149 (MQ=255)
                                                                   cGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACggg                                                                               <  2:3986183/149‑1 (MQ=255)
                                                                     tctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTg                                                                             <  1:353066/149‑1 (MQ=255)
                                                                      ctcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGc                                                                            <  1:164088/149‑1 (MQ=255)
                                                                       tcCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCg                                                                           >  1:4023589/1‑149 (MQ=255)
                                                                        ccATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGt                                                                          >  2:459885/1‑149 (MQ=255)
                                                                            catAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTAt                                                                      <  1:3679083/149‑1 (MQ=255)
                                                                              tAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATcc                                                                    <  1:3192299/149‑1 (MQ=255)
                                                                              tAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATcc                                                                    >  2:807931/1‑149 (MQ=255)
                                                                                 cGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGc                                                                 <  2:1900528/149‑1 (MQ=255)
                                                                                      tGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTc                                                            <  2:1884509/149‑1 (MQ=255)
                                                                                       gCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCg                                                           <  1:880560/149‑1 (MQ=255)
                                                                                         aGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGa                                                         >  1:3683726/1‑149 (MQ=255)
                                                                                         aGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATCTTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGa                                                         >  2:4079106/1‑149 (MQ=255)
                                                                                          gACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                                                        >  1:2880141/1‑149 (MQ=255)
                                                                                          gACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGAc                                                        <  1:1871396/149‑1 (MQ=255)
                                                                                              ggTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTg                                                    >  2:2523814/1‑149 (MQ=255)
                                                                                               gTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTgg                                                   >  1:2807166/1‑149 (MQ=255)
                                                                                                   gTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATa                                               <  2:3692575/149‑1 (MQ=255)
                                                                                                    tAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAg                                              >  1:4106942/1‑149 (MQ=255)
                                                                                                    tAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTATAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAg                                              <  2:3075854/149‑1 (MQ=255)
                                                                                                         cACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCa                                         >  2:1821432/1‑149 (MQ=255)
                                                                                                               cGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACg                                   >  2:1471719/1‑149 (MQ=255)
                                                                                                                ggCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                  >  1:2564476/1‑149 (MQ=255)
                                                                                                                ggCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                  >  2:3614289/1‑149 (MQ=255)
                                                                                                                ggCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGc                                  <  1:1441676/149‑1 (MQ=255)
                                                                                                                 gCTTCCACAGCAGAATCCGCCGGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCt                                 >  1:848609/1‑149 (MQ=255)
                                                                                                                   ttCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGa                               <  1:2741709/149‑1 (MQ=255)
                                                                                                                      cacaGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGc                            <  1:2305979/149‑1 (MQ=255)
                                                                                                                            agAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGct                      <  1:1821030/149‑1 (MQ=255)
                                                                                                                             gAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGctc                     >  2:539654/1‑149 (MQ=255)
                                                                                                                                  cgccgGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCa                <  2:3683706/149‑1 (MQ=255)
                                                                                                                                    ccgGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgt              >  2:329279/1‑149 (MQ=255)
                                                                                                                                     cgGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgtt             <  1:3598091/149‑1 (MQ=255)
                                                                                                                                     cgGGTTATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgtt             >  1:100668/1‑149 (MQ=255)
                                                                                                                                       ggTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgttgt           >  1:783182/1‑149 (MQ=255)
                                                                                                                                       ggTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAgttgt           >  2:2720721/1‑149 (MQ=255)
                                                                                                                                         ttATCCACAGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTc         <  1:3352614/149‑1 (MQ=255)
                                                                                                                                                ccGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCt  <  2:264440/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GAGTTGAACTGACAAACGTTTTACCACTGTCTGGCGTCGCACCGGTGATCATTAGAATGTTATTCTCCGTCTCCATCATAGCGAAATGCAGACTGGTTCGTAGCGCACGTACGGCTTCCACAGCAGAATCCGCCGGGTTATCCACCGCCAGGAAGGGGATATTTTTAGTACGATGGCGCTGCTGATTAGAAAATAAATTTTTCTTACGCAGACGGGTGCGTTTATCCAGCCACTCGGACATTGGGATAGTGGCATAAACGCTGATGCCGTGCTCTTCCAGTTGTTCCGGGGCT  >  CP009273/1038013‑1038305

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