Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,454,297 T→G 64.2% E39D (GAA→GAC yfcZ ← UPF0381 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,454,2970TG64.2% 470.3 / inf 310E39D (GAA→GACyfcZUPF0381 family protein
Reads supporting (aligned to +/- strand):  ref base T (54/57);  new base G (88/111);  total (142/168)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 313 total aligned reads displayed.

AAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCAAGGTG  >  CP009273/2454149‑2454445
                                                                                                                                                    |                                                                                                                                                    
aaaTCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCt                                                                                                                                                      <  2:590053/149‑1 (MQ=255)
  aTCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTc                                                                                                                                                    <  1:4278754/149‑1 (MQ=255)
  aTCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTc                                                                                                                                                    >  2:1589816/1‑149 (MQ=255)
    cAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTg                                                                                                                                                  >  1:778096/1‑149 (MQ=255)
     aGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTgc                                                                                                                                                 <  1:2202762/149‑1 (MQ=255)
           tCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGtt                                                                                                                                           <  2:253250/149‑1 (MQ=255)
            cAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGttt                                                                                                                                          >  2:4474874/1‑149 (MQ=255)
            cAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGttt                                                                                                                                          <  2:1443231/149‑1 (MQ=255)
                ttCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCaa                                                                                                                                      >  2:3768536/1‑149 (MQ=255)
                ttCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGATCTGCGTCTGCGCGGATTGCaa                                                                                                                                      >  2:1009913/1‑149 (MQ=255)
                 tCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCaaa                                                                                                                                     <  2:3458366/149‑1 (MQ=255)
                  cGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAac                                                                                                                                    >  1:38212/1‑149 (MQ=255)
                    cAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAacac                                                                                                                                  <  2:3797131/149‑1 (MQ=255)
                    cAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAacac                                                                                                                                  <  1:1041127/149‑1 (MQ=255)
                      ggCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcgc                                                                                                                                >  2:260182/1‑149 (MQ=255)
                       gCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACAcgcg                                                                                                                               >  2:2335115/1‑149 (MQ=255)
                         aaaaGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGc                                                                                                                             <  1:3099680/149‑1 (MQ=255)
                            aGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGt                                                                                                                          <  1:2726953/149‑1 (MQ=255)
                             gTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTa                                                                                                                         >  1:753949/1‑149 (MQ=255)
                                aaaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGa                                                                                                                      <  2:3921567/149‑1 (MQ=255)
                                 aaTCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGAc                                                                                                                     >  1:1637135/1‑149 (MQ=255)
                                     gATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCgg                                                                                                                 >  1:2183744/1‑149 (MQ=255)
                                      atatCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGt                                                                                                                >  1:4389776/1‑149 (MQ=255)
                                      atatCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGt                                                                                                                <  1:3615379/149‑1 (MQ=255)
                                        atCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGc                                                                                                              >  1:3455962/1‑149 (MQ=255)
                                         tCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCa                                                                                                             <  1:551387/149‑1 (MQ=255)
                                         tCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCGGTACGACGCGGTGCa                                                                                                             >  2:4369479/1‑149 (MQ=255)
                                           cAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGt                                                                                                           >  1:1421221/1‑149 (MQ=255)
                                           cAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGt                                                                                                           <  1:1568056/149‑1 (MQ=255)
                                           cAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGt                                                                                                           <  1:3130837/149‑1 (MQ=255)
                                           cAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGt                                                                                                           <  2:2850116/149‑1 (MQ=255)
                                            aGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTc                                                                                                          <  2:3539358/149‑1 (MQ=255)
                                             ggCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCg                                                                                                         >  1:2997475/1‑149 (MQ=255)
                                               cgcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGa                                                                                                       <  2:4350351/149‑1 (MQ=255)
                                                gcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAg                                                                                                      >  1:2596656/1‑149 (MQ=255)
                                                gcACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAg                                                                                                      <  1:1533388/149‑1 (MQ=255)
                                                   cacCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTg                                                                                                   <  2:3110970/149‑1 (MQ=255)
                                                     ccGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTc                                                                                                 >  2:494562/1‑149 (MQ=255)
                                                      cGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTcc                                                                                                <  1:1849261/149‑1 (MQ=255)
                                                        tCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCAt                                                                                              <  1:1252620/149‑1 (MQ=255)
                                                          aGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATaa                                                                                            <  1:691848/149‑1 (MQ=255)
                                                          aGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGACAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATaa                                                                                            >  2:261437/1‑149 (MQ=255)
                                                           gACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAAt                                                                                           <  2:1750663/149‑1 (MQ=255)
                                                             ctctTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATgg                                                                                         >  1:3710204/1‑149 (MQ=255)
                                                              tctTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGt                                                                                        >  1:3663401/1‑149 (MQ=255)
                                                                 tCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGcc                                                                                     <  2:4140238/149‑1 (MQ=255)
                                                                   aGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCaa                                                                                   <  1:2100116/149‑1 (MQ=255)
                                                                     tGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACa                                                                                 <  1:3543525/149‑1 (MQ=255)
                                                                      gAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACAt                                                                                >  2:190776/1‑149 (MQ=255)
                                                                         aGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCa                                                                             >  1:4127457/1‑149 (MQ=255)
                                                                         aGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCa                                                                             >  1:3385524/1‑149 (MQ=255)
                                                                           tCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATa                                                                           >  2:2799333/1‑149 (MQ=255)
                                                                              gggTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATAcag                                                                        >  1:1358342/1‑149 (MQ=255)
                                                                                gTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATAcagca                                                                      <  2:3251982/149‑1 (MQ=255)
                                                                                gTGATTTTGCAAGGTTCAGAGTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATAcagca                                                                      <  1:1606824/149‑1 (MQ=255)
                                                                                 tGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATAcagcag                                                                     >  1:4030111/1‑149 (MQ=255)
                                                                                   aTTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAgcagca                                                                   >  2:1962762/1‑149 (MQ=255)
                                                                                     tttGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCaaa                                                                 <  2:809448/149‑1 (MQ=255)
                                                                                      ttGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAAc                                                                <  1:2888315/149‑1 (MQ=255)
                                                                                       tGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAAcc                                                               <  2:2430732/149‑1 (MQ=255)
                                                                                          aaGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCggg                                                            >  2:361778/1‑149 (MQ=255)
                                                                                          aaGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCggg                                                            >  1:2656269/1‑149 (MQ=255)
                                                                                           aGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCgggg                                                           >  1:4174543/1‑149 (MQ=255)
                                                                                           aGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCgggg                                                           <  2:4201459/149‑1 (MQ=255)
                                                                                            ggTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGt                                                          >  2:1690247/1‑149 (MQ=255)
                                                                                              ttCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGttt                                                        <  2:3039743/149‑1 (MQ=255)
                                                                                                cAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtca                                                      >  1:1319572/1‑149 (MQ=255)
                                                                                                 aGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcat                                                     >  2:457434/1‑149 (MQ=255)
                                                                                                 aGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcat                                                     <  1:4521435/149‑1 (MQ=255)
                                                                                                  gATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatc                                                    <  2:2919046/149‑1 (MQ=255)
                                                                                                   aTTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTtcatca                                                   >  2:3313750/1‑149 (MQ=255)
                                                                                                    ttCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAg                                                  >  2:315438/1‑149 (MQ=255)
                                                                                                      ccACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCa                                                >  2:2955581/1‑149 (MQ=255)
                                                                                                       cACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCAc                                               <  2:334261/149‑1 (MQ=255)
                                                                                                        aCGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACt                                              >  1:4108711/1‑149 (MQ=255)
                                                                                                             aCGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACAtt                                         >  1:718204/1‑149 (MQ=255)
                                                                                                             aCGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACAtt                                         <  2:2283179/149‑1 (MQ=255)
                                                                                                                ggCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTAc                                      >  1:3732540/1‑149 (MQ=255)
                                                                                                                 gCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACt                                     <  2:919004/149‑1 (MQ=255)
                                                                                                                  cTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTc                                    <  1:2399605/149‑1 (MQ=255)
                                                                                                                  cTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTc                                    <  2:4127457/149‑1 (MQ=255)
                                                                                                                    ttttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCAt                                  <  1:4117627/149‑1 (MQ=255)
                                                                                                                     tttCAGTCAGCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATc                                 <  1:3590229/149‑1 (MQ=255)
                                                                                                                          gtcagCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGa                            <  1:3225604/149‑1 (MQ=255)
                                                                                                                           tcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAt                           <  2:1211032/149‑1 (MQ=255)
                                                                                                                           tcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAt                           <  1:4350036/149‑1 (MQ=255)
                                                                                                                           tcagCGCCGCCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAt                           >  2:4048042/1‑149 (MQ=255)
                                                                                                                            cagCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAtt                          <  1:1135464/149‑1 (MQ=255)
                                                                                                                             agCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGAttt                         <  1:3059129/149‑1 (MQ=255)
                                                                                                                               cgccgcCAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTcc                       <  2:210102/149‑1 (MQ=255)
                                                                                                                                 ccgccAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCctc                     >  1:2278112/1‑149 (MQ=255)
                                                                                                                                   gccAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTg                   <  2:214663/149‑1 (MQ=255)
                                                                                                                                    ccAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGt                  >  2:3996453/1‑149 (MQ=255)
                                                                                                                                     cAGCGTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTa                 <  1:406829/149‑1 (MQ=255)
                                                                                                                                         gTTTGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTAtgtg             <  2:2731075/149‑1 (MQ=255)
                                                                                                                                          tttGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGc            >  1:832924/1‑149 (MQ=255)
                                                                                                                                          tttGCTCTGCGTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGc            >  1:3281309/1‑149 (MQ=255)
                                                                                                                                           ttGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCa           <  1:3916741/149‑1 (MQ=255)
                                                                                                                                               tctGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCTCTTCCTTTACTCATCGTTGATTTCCTCTGTATGTGCACCCa       <  2:223079/149‑1 (MQ=255)
                                                                                                                                                   cTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCAAGGt   <  2:778096/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AAATCAGCATTTCAGCTTCGCAGGCAAAAGTGAAATCGATATCCAGGCGCACACCGTCAGACTCTTCAGTGAAAGTCGGGGTGATTTTGCAAGGTTCAGATTCCACGCTACGGGCTTTTTCAGTCAGCGCCGCCAGCGTTTGCTCTGCTTCTGCGCGGTTTGCAAACACGCGGCTGTAAGACGCGGTGCAGTCGGAGTTGTCCATAATGGTGCCAACATCCATACAGCAGCAAACCGGGGTTTCATCAGCACTACATTTACTCATCGTTGATTTCCTCTGTATGTGCACCCAAGGTG  >  CP009273/2454149‑2454445

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: