Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,122,653 C→T 100% W160* (TGG→TAG)  wcaE ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,122,6530CT97.4% 760.7 / 6.7 230W160* (TGG→TAG) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base T (113/111);  total (118/112)
Fisher's exact test for biased strand distribution p-value = 2.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 230 total aligned reads displayed.

ACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGTT  >  CP009273/2122507‑2122798
                                                                                                                                                  |                                                                                                                                                 
acaATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAt                                                                                                                                                 <  1:2972934/149‑1 (MQ=255)
 caATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtt                                                                                                                                                >  2:2650247/1‑149 (MQ=255)
  aaTTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAttt                                                                                                                                               <  1:2434066/149‑1 (MQ=255)
  aaTTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAttt                                                                                                                                               <  1:43698/149‑1 (MQ=255)
   aTTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtttt                                                                                                                                              <  1:516321/149‑1 (MQ=255)
   aTTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTAtttt                                                                                                                                              <  1:3276240/149‑1 (MQ=255)
       cATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGACGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCa                                                                                                                                          >  2:1535864/1‑149 (MQ=255)
          attattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGc                                                                                                                                       <  2:1677655/149‑1 (MQ=255)
           ttattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGcc                                                                                                                                      >  2:3527916/1‑149 (MQ=255)
            tattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATACCCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCg                                                                                                                                     <  2:2445838/149‑1 (MQ=255)
            tattGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCg                                                                                                                                     >  2:3120865/1‑149 (MQ=255)
                     tAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATCACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAggg                                                                                                                            <  2:2117071/149‑1 (MQ=255)
                        aTACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGaaa                                                                                                                         <  1:3124079/149‑1 (MQ=255)
                        aTACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGaaa                                                                                                                         >  1:2193854/1‑149 (MQ=255)
                             ccacccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAatat                                                                                                                    >  1:1569017/1‑149 (MQ=255)
                                cccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGc                                                                                                                 >  1:1654796/1‑149 (MQ=255)
                                cccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGc                                                                                                                 >  2:214288/1‑149 (MQ=255)
                                cccaTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGc                                                                                                                 <  1:3955673/149‑1 (MQ=255)
                                  caTGGAAAATTCAGACACCATGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCt                                                                                                               <  2:371222/149‑1 (MQ=255)
                                  caTGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCt                                                                                                               >  2:4085191/1‑149 (MQ=255)
                                    tGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCtga                                                                                                             <  1:1667663/149‑1 (MQ=255)
                                       aaaaTTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCTGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCtgatga                                                                                                          <  2:3774358/149‑1 (MQ=255)
                                        aaaTTCAGACACCAGTCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCAGCCAGCGCGTAGTCGTAATAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGAc                                                                                                         <  2:3272702/149‑1 (MQ=255)
                                        aaaTTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGAc                                                                                                         <  2:1605415/149‑1 (MQ=255)
                                         aaTTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACt                                                                                                        <  2:3208814/149‑1 (MQ=255)
                                             cAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCg                                                                                                    >  1:362782/1‑149 (MQ=255)
                                              aGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCgg                                                                                                   >  2:1891276/1‑149 (MQ=255)
                                                   ccAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAgg                                                                                              >  1:3845470/1‑149 (MQ=255)
                                                    cAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGc                                                                                             <  1:1830023/149‑1 (MQ=255)
                                                    cAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGc                                                                                             <  1:1512733/149‑1 (MQ=255)
                                                    cAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGc                                                                                             <  1:709240/149‑1 (MQ=255)
                                                      ggCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCtg                                                                                           >  2:573443/1‑149 (MQ=255)
                                                          aTTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGAt                                                                                       >  1:963390/1‑149 (MQ=255)
                                                             gAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATaaa                                                                                    <  1:1038563/149‑1 (MQ=255)
                                                             gAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATaaa                                                                                    >  1:2349334/1‑149 (MQ=255)
                                                               gTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAata                                                                                  >  2:1829851/1‑149 (MQ=255)
                                                                ttttttAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCTGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAatat                                                                                 <  2:2305526/149‑1 (MQ=255)
                                                                 tttttAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAAtata                                                                                >  2:1911535/1‑149 (MQ=255)
                                                                   tttAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATAcc                                                                              >  2:1309904/1‑149 (MQ=255)
                                                                      aaaTGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGc                                                                           >  2:736928/1‑149 (MQ=255)
                                                                      aaaTGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGc                                                                           >  2:2767342/1‑149 (MQ=255)
                                                                       aaTGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGcc                                                                          <  1:2782307/149‑1 (MQ=255)
                                                                       aaTGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGcc                                                                          >  2:2309972/1‑149 (MQ=255)
                                                                             tAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGttt                                                                    <  1:1060287/149‑1 (MQ=255)
                                                                             tAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGttt                                                                    <  1:3589424/149‑1 (MQ=255)
                                                                              aaCCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTg                                                                   <  1:1696051/149‑1 (MQ=255)
                                                                               aCCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAAAATACCAGCCCGGTTTgg                                                                  >  2:3696856/1‑149 (MQ=255)
                                                                                  ggCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGgcgc                                                               >  1:908858/1‑149 (MQ=255)
                                                                                  ggCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGgcgc                                                               <  1:2016245/149‑1 (MQ=255)
                                                                                   gCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCt                                                              <  1:803336/149‑1 (MQ=255)
                                                                                     tttATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTAc                                                            >  1:1315875/1‑149 (MQ=255)
                                                                                     tttATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTAc                                                            <  1:2369329/149‑1 (MQ=255)
                                                                                       tataCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGt                                                          <  1:4061896/149‑1 (MQ=255)
                                                                                         taCATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGttt                                                        >  2:149938/1‑149 (MQ=255)
                                                                                            aTTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAt                                                     >  2:124861/1‑149 (MQ=255)
                                                                                             ttttGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAAtt                                                    <  2:204933/149‑1 (MQ=255)
                                                                                                tGGCTGCCAGCGCGTAGTCTTAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTa                                                 <  2:4000055/149‑1 (MQ=255)
                                                                                                 ggCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTAt                                                <  1:2941641/149‑1 (MQ=255)
                                                                                                  gCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATg                                               <  1:1911535/149‑1 (MQ=255)
                                                                                                   cTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGc                                              <  1:1533586/149‑1 (MQ=255)
                                                                                                      ccAGCGCGTAGTCTGAAGAAAATTTATTTTCCATGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTTATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCg                                           <  2:3136268/149‑1 (MQ=255)
                                                                                                      ccAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCg                                           >  2:3104066/1‑149 (MQ=255)
                                                                                                      ccAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCg                                           >  1:1187772/1‑149 (MQ=255)
                                                                                                       cAGCGCGTAGTCTGAATAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGt                                          <  2:3340901/149‑1 (MQ=255)
                                                                                                               aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCaaaa                                  <  1:3527916/149‑1 (MQ=255)
                                                                                                               aGTCGGAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCaaaa                                  <  1:1322353/149‑1 (MQ=255)
                                                                                                                  cGGAAGAAACTTTATATTTCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATcc                               <  1:89330/149‑1 (MQ=255)
                                                                                                                   ggAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCca                              >  2:756706/1‑149 (MQ=255)
                                                                                                                   ggAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCca                              >  2:2219436/1‑149 (MQ=255)
                                                                                                                   ggAAGAAACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCca                              >  2:3979958/1‑149 (MQ=255)
                                                                                                                     aagaaACTTTATATTCCTGGTCATAACGCCATTTTTTCAAGCCGGATACAGTGAAAAATATCGCCTGATGACTGGCGTGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagc                            <  2:482525/149‑1 (MQ=255)
                                                                                                                     aagaaACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagc                            >  1:3173991/1‑149 (MQ=255)
                                                                                                                     aagaaACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagc                            <  1:144816/149‑1 (MQ=255)
                                                                                                                      agaaACTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCcagca                           <  1:3978545/149‑1 (MQ=255)
                                                                                                                         aaCTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAgcg                        >  2:3376599/1‑149 (MQ=255)
                                                                                                                          aCTTTATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAgcgc                       >  2:2672713/1‑149 (MQ=255)
                                                                                                                             ttATATTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATc                    >  2:2869366/1‑149 (MQ=255)
                                                                                                                               atatTCCAGTTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTTATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGc                  <  2:3389406/149‑1 (MQ=255)
                                                                                                                                tatTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGcc                 <  2:776835/149‑1 (MQ=255)
                                                                                                                                 atTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCATCAGCGCATCGCCg                >  2:255463/1‑149 (MQ=255)
                                                                                                                                 atTCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCg                >  2:10305/1‑149 (MQ=255)
                                                                                                                                 atTCCAGGTCATAACGATATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCg                <  2:2976042/149‑1 (MQ=255)
                                                                                                                                   tCCAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGt              >  1:721248/1‑149 (MQ=255)
                                                                                                                                    ccAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTg             >  2:532443/1‑149 (MQ=255)
                                                                                                                                     cAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGa            >  1:231071/1‑149 (MQ=255)
                                                                                                                                     cAGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGa            <  1:805268/149‑1 (MQ=255)
                                                                                                                                      aGGTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGat           <  1:3939094/149‑1 (MQ=255)
                                                                                                                                       gttCATAACGCTATTTTTTCAAGCCGGATACAGGTAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatc          <  2:2243456/147‑1 (MQ=255)
                                                                                                                                       ggTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatc          >  2:1399802/1‑149 (MQ=255)
                                                                                                                                        gTCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatca         <  2:1216340/149‑1 (MQ=255)
                                                                                                                                         tCATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatcat        <  1:3221107/149‑1 (MQ=255)
                                                                                                                                          cATAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGatcatc       >  2:2835198/1‑149 (MQ=255)
                                                                                                                                           aTAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGAtcatca      <  2:1652597/149‑1 (MQ=255)
                                                                                                                                           aTAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGAtcatca      >  2:611051/1‑149 (MQ=255)
                                                                                                                                            tAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCAc     >  1:2684211/1‑149 (MQ=255)
                                                                                                                                            tAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCAc     >  1:2482697/1‑149 (MQ=255)
                                                                                                                                            tAACGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCAc     <  1:145303/149‑1 (MQ=255)
                                                                                                                                              aCGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGt   <  1:3489858/149‑1 (MQ=255)
                                                                                                                                              aCGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGt   <  1:532443/149‑1 (MQ=255)
                                                                                                                                              aCGCTATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGt   >  2:206538/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ACAATTCCATATTATTGGTGGTAGATACCCCACCCATGGAAAATTCAGACACCAGGCCATTGAGTTTTTTAAATGCATAACCGGCTTTATACATTTTGGCTGCCAGCGCGTAGTCGGAAGAAACTTTATATTCCAGGTCATAACGCCATTTTTTCAAGCCGGATACAGGGAAAAATATCGCCTGATGACTGGCGGGCAGGCTGTGATAAATATACCAGCCCGGTTTGGCGCTACGTTTAATTTTATGCCCGTCGCCAAAATCCAGCAGCGCATCGCCGGTGATCATCACGTT  >  CP009273/2122507‑2122798

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: