Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,850,839 G→T 20.9% P309Q (CCG→CAG)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,850,8390GT20.9% 825.9 / 180.4 300P309Q (CCG→CAG) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base G (106/131);  new base T (28/35);  total (134/166)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 306 total aligned reads displayed.

TGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGT  >  CP009273/3850691‑3850986
                                                                                                                                                    |                                                                                                                                                   
 ggCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCgg                                                                                                                                                    >  1:3364202/1‑149 (MQ=255)
   ctctTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGaa                                                                                                                                                  <  1:2781756/149‑1 (MQ=255)
   ctctTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGaa                                                                                                                                                  >  1:3588098/1‑149 (MQ=255)
         cACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCggg                                                                                                                                            <  1:3793028/149‑1 (MQ=255)
          aCGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGt                                                                                                                                           >  2:1977771/1‑149 (MQ=255)
                aGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATgtc                                                                                                                                     >  2:2254295/1‑149 (MQ=255)
                aGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATgtc                                                                                                                                     <  2:2674308/149‑1 (MQ=255)
                    tatTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCa                                                                                                                                 >  2:3613056/1‑149 (MQ=255)
                       tttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATAc                                                                                                                              <  2:2495948/149‑1 (MQ=255)
                       tttCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATAc                                                                                                                              >  2:3943056/1‑149 (MQ=255)
                       tgtCCCCCGGCAGAATCTCTGGATTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATAc                                                                                                                              <  2:2512017/147‑1 (MQ=255)
                           ccccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcg                                                                                                                          >  1:2831923/1‑149 (MQ=255)
                            cccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcgc                                                                                                                         >  2:3082655/1‑149 (MQ=255)
                             ccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATcgcc                                                                                                                        >  2:498727/1‑149 (MQ=255)
                             ccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcgcc                                                                                                                        >  2:4631641/1‑149 (MQ=255)
                             ccGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATcgcc                                                                                                                        >  2:298487/1‑149 (MQ=255)
                                gCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCgccgcc                                                                                                                     >  1:5031260/1‑149 (MQ=255)
                                gCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCgccgcc                                                                                                                     >  2:1187551/1‑149 (MQ=255)
                                gCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCgccgcc                                                                                                                     >  1:988953/1‑149 (MQ=255)
                                 cAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGccc                                                                                                                    >  1:5091634/1‑149 (MQ=255)
                                   gAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCt                                                                                                                  <  2:1508782/149‑1 (MQ=255)
                                    aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCtt                                                                                                                 <  2:801501/149‑1 (MQ=255)
                                    aaTCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCtt                                                                                                                 >  1:1991515/1‑149 (MQ=255)
                                       ctctGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTAc                                                                                                              <  2:2637204/149‑1 (MQ=255)
                                          tGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGa                                                                                                           <  1:1805115/149‑1 (MQ=255)
                                           ggTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGAt                                                                                                          <  2:901461/149‑1 (MQ=255)
                                             tttttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCa                                                                                                        >  1:2424701/1‑149 (MQ=255)
                                               tttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGc                                                                                                      <  2:1707017/149‑1 (MQ=255)
                                                ttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCt                                                                                                     >  2:2782104/1‑149 (MQ=255)
                                                ttCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCt                                                                                                     <  1:4263226/149‑1 (MQ=255)
                                                 tCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCtt                                                                                                    <  1:1254596/149‑1 (MQ=255)
                                                 tCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCtt                                                                                                    >  2:1063114/1‑149 (MQ=255)
                                                     tatCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgt                                                                                                >  1:2168504/1‑149 (MQ=255)
                                                      atCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAtgtg                                                                                               <  2:3854185/149‑1 (MQ=255)
                                                        cTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTgcg                                                                                             >  2:434303/1‑149 (MQ=255)
                                                               cgcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatca                                                                                      <  1:4788018/149‑1 (MQ=255)
                                                                gcCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCcatcat                                                                                     >  2:3934879/1‑149 (MQ=255)
                                                                  catcatTTTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGACGCTCCCCGTGCGGCGAACGGATGATAAGCGGAATGCGGGTGATGTCGTCATACATCGTCGCCCCTTTACTTTTCAGCTTATGTGCGCCCTTCAttt                                                                                   <  2:2383213/149‑1 (MQ=255)
                                                                  catcatTGTCGGCAGTACATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAttt                                                                                   >  1:1714348/1‑149 (MQ=255)
                                                                  catcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCAttt                                                                                   <  2:2819192/149‑1 (MQ=255)
                                                                   atcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTc                                                                                  >  1:694014/1‑149 (MQ=255)
                                                                    tcatTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTACCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGGGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCg                                                                                 >  2:3556364/1‑149 (MQ=255)
                                                                      atTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGcc                                                                               >  1:24363/1‑149 (MQ=255)
                                                                      atTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTAGGTGCGCCCATCATTTCGcc                                                                               >  2:2138278/1‑149 (MQ=255)
                                                                       ttGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCg                                                                              >  1:2229789/1‑149 (MQ=255)
                                                                       ttGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCg                                                                              >  1:247051/1‑149 (MQ=255)
                                                                        tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                             <  1:1429865/149‑1 (MQ=255)
                                                                        tGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGt                                                                             >  2:5188864/1‑149 (MQ=255)
                                                                           cGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGAt                                                                          <  2:4719803/149‑1 (MQ=255)
                                                                             gCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCg                                                                        <  2:2883760/149‑1 (MQ=255)
                                                                               aGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGa                                                                      >  1:1644313/1‑149 (MQ=255)
                                                                                 tAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAgg                                                                    >  1:1565877/1‑149 (MQ=255)
                                                                                  aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGACTGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGt                                                                   >  2:2612206/1‑149 (MQ=255)
                                                                                  aaaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGt                                                                   >  2:3874817/1‑149 (MQ=255)
                                                                                   aaTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGta                                                                  >  1:4628319/1‑149 (MQ=255)
                                                                                    aTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                 >  1:2020629/1‑149 (MQ=255)
                                                                                    aTCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGtat                                                                 >  2:1629451/1‑149 (MQ=255)
                                                                                      cGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataa                                                               <  2:1742781/149‑1 (MQ=255)
                                                                                      cGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataa                                                               <  2:4336921/149‑1 (MQ=255)
                                                                                       gATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaa                                                              <  2:3535987/149‑1 (MQ=255)
                                                                                         tatGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaata                                                            >  1:1550044/1‑149 (MQ=255)
                                                                                         tatGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTataaata                                                            >  1:891334/1‑149 (MQ=255)
                                                                                           tgactgacTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAAc                                                          <  2:4152502/149‑1 (MQ=255)
                                                                                                gactgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACg                                                     <  2:331369/149‑1 (MQ=255)
                                                                                                 actgGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATTCATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGt                                                    <  2:1962008/149‑1 (MQ=255)
                                                                                                    ggCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAtt                                                 >  1:1602690/1‑149 (MQ=255)
                                                                                                     gCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTAttt                                                <  1:2591580/149‑1 (MQ=255)
                                                                                                          tCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGAGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACg                                           >  2:436022/1‑149 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                          <  1:1015376/149‑1 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                          <  1:578191/149‑1 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGt                                          <  1:4836218/149‑1 (MQ=255)
                                                                                                           cGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGAATTTTCACGt                                          >  1:3960217/1‑149 (MQ=255)
                                                                                                            gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                         <  1:3215839/149‑1 (MQ=255)
                                                                                                            gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                         <  1:4745428/149‑1 (MQ=255)
                                                                                                            gACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGtt                                         >  1:4108021/1‑149 (MQ=255)
                                                                                                             aCCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTg                                        >  1:2339426/1‑149 (MQ=255)
                                                                                                               cTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGct                                      <  1:2372256/149‑1 (MQ=255)
                                                                                                               cTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGct                                      <  2:160973/149‑1 (MQ=255)
                                                                                                               cTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGct                                      <  1:4041904/149‑1 (MQ=255)
                                                                                                                tGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctc                                     <  2:4862444/149‑1 (MQ=255)
                                                                                                                 gTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGctct                                    <  2:400984/149‑1 (MQ=255)
                                                                                                                  tCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTg                                   >  2:3303882/1‑149 (MQ=255)
                                                                                                                   cgccgcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTgg                                  <  2:1374053/149‑1 (MQ=255)
                                                                                                                    gccgcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGc                                 >  2:3485550/1‑149 (MQ=255)
                                                                                                                       gcTCCCCTTGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGtt                              <  2:1276237/149‑1 (MQ=255)
                                                                                                                       gcTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGtt                              >  1:2147264/1‑149 (MQ=255)
                                                                                                                          ccccTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAg                           <  1:3999971/149‑1 (MQ=255)
                                                                                                                          ccccTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAg                           <  1:2670996/149‑1 (MQ=255)
                                                                                                                             cTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCa                        <  2:335774/149‑1 (MQ=255)
                                                                                                                             cTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCa                        <  1:3613118/149‑1 (MQ=255)
                                                                                                                              ttGCGGCGAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCAt                       >  2:2074864/1‑149 (MQ=255)
                                                                                                                              ttGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCAt                       <  1:1848842/149‑1 (MQ=255)
                                                                                                                                 cggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtga                    <  1:109834/149‑1 (MQ=255)
                                                                                                                                 cggcgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtga                    <  2:4512338/149‑1 (MQ=255)
                                                                                                                                  ggcgAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtgat                   <  2:4237288/149‑1 (MQ=255)
                                                                                                                                    cgAACGGATGATCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATtgatga                 <  2:2006396/149‑1 (MQ=255)
                                                                                                                                          gatgatCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTc           <  2:62845/149‑1 (MQ=255)
                                                                                                                                            tgatCAGCTGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCg         >  2:943803/1‑149 (MQ=255)
                                                                                                                                             gatCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGa        <  1:2553735/149‑1 (MQ=255)
                                                                                                                                                   cGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGt  <  2:811038/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGCTCTTTCACGGCAAGGATATTTTCCCCCGGCAGAATCTCTGGTTTTTCAATATCTGCCAGCGCCATCATTGTCGGCAGTAAATCGATATGACTGACTGGCGTATCGACCTGTCGCCGCTCCCCTTGCGGCGAACGGATGATCAGCGGAATGCGGGTGATGTCGTCATACATCGCCGCCCCTTTACTGATCAGCTTATGTGCGCCCATCATTTCGCCGTGATCGGAGGTATAAATAACCCACGTATTTTCACGTTGCTCTGGCGTTAAGGCATTGATGACCCGTCCGATTTGGT  >  CP009273/3850691‑3850986

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: