Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 522,859 T→G 19.7% S1381R (AGT→AGG rhsD → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273522,8590TG19.7% 474.6 / 155.5 238S1381R (AGT→AGGrhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base T (103/88);  new base G (21/26);  total (124/114)
Fisher's exact test for biased strand distribution p-value = 3.28e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 239 total aligned reads displayed.

GATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCT  >  CP009273/522710‑523007
                                                                                                                                                     |                                                                                                                                                    
gATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGt                                                                                                                                                      >  2:596490/1‑150 (MQ=255)
  tAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTaa                                                                                                                                                    >  1:2103772/1‑150 (MQ=255)
   aaTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTaaa                                                                                                                                                   <  1:447008/150‑1 (MQ=255)
    aTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAg                                                                                                                                                  >  1:320500/1‑150 (MQ=255)
     tttAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGa                                                                                                                                                 >  1:1031364/1‑150 (MQ=255)
     tttAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGa                                                                                                                                                 >  1:1031300/1‑150 (MQ=255)
      ttAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGat                                                                                                                                                >  1:550015/1‑150 (MQ=255)
              ttattCATATCCTATTAAAATAAGTTGTCCCACATGTTCAGATTGATTTGCTGATTATAATACTCTCGGTTAAGATATTAATGGCTTTTATTTTGATGATTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAATTAAAGATAATAATCt                                                                                                                                        <  2:4251449/146‑1 (MQ=255)
                aGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTg                                                                                                                                      >  1:939520/1‑150 (MQ=255)
                   cATATCCTCTTAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACTGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAg                                                                                                                                   <  2:2091674/150‑1 (MQ=255)
                   cATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAg                                                                                                                                   >  1:3559049/1‑150 (MQ=255)
                   cATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTAAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAg                                                                                                                                   >  1:965628/1‑150 (MQ=255)
                    atatCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGc                                                                                                                                  <  1:1092909/150‑1 (MQ=255)
                      atCCTCTGAAAAGAAAATTTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCAt                                                                                                                                >  2:1179863/1‑150 (MQ=255)
                       tCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCAtt                                                                                                                               >  2:2889750/1‑150 (MQ=255)
                         ctctGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCAtttt                                                                                                                             >  2:2730288/1‑150 (MQ=255)
                          tctGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTa                                                                                                                            <  1:3778315/150‑1 (MQ=255)
                           ctGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAt                                                                                                                           <  1:3248598/150‑1 (MQ=255)
                           ctGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAt                                                                                                                           >  1:3373284/1‑150 (MQ=255)
                            tGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtc                                                                                                                          >  1:1327028/1‑150 (MQ=255)
                             taaaataaaatGTACCAAAGGTAAAGATATTGTTGCTGGTTATCATAATCACGGTGCAGATAGTCATGGCGATTATTTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtct                                                                                                                         <  2:2930051/149‑1 (MQ=255)
                             gaaaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAATTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtct                                                                                                                         >  2:2443686/1‑150 (MQ=255)
                              aaaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtctc                                                                                                                        >  1:3486348/1‑150 (MQ=255)
                                 agaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCa                                                                                                                     >  2:2432257/1‑150 (MQ=255)
                                 agaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGAGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCa                                                                                                                     <  1:3977316/150‑1 (MQ=255)
                                    aaaTGTCCCACAGGTACAGATATACAAGCTGCTTATCATACTCAAGGAGCAGATAGACATGTGGATTATGTTGATCAAATTTTAACAAGTAGCGAAAAAAATCTACTAAGAAGTAAATAAAATAAACTTGAAGCAAATTATCTCGCAaca                                                                                                                  >  1:1826436/1‑150 (MQ=255)
                                     aaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAacac                                                                                                                 <  1:4499103/150‑1 (MQ=255)
                                       tGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCt                                                                                                               <  1:1653740/150‑1 (MQ=255)
                                        gTCCCACAGGTACAGATAGAGATGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACAACTg                                                                                                              >  1:2961524/1‑150 (MQ=255)
                                         tcacacaGGTACATATAGAGTTGCTGATTATAATACTCACGTTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGa                                                                                                             <  2:1704834/147‑1 (MQ=255)
                                         tCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGa                                                                                                             >  2:4321381/1‑150 (MQ=255)
                                         tCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGa                                                                                                             >  2:1025377/1‑150 (MQ=255)
                                           ccACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTGGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATg                                                                                                           >  2:2572557/1‑150 (MQ=255)
                                              caGTTACAGATTGAGTTGATTCTTATCATACTCACGTTGCAGATAGTCATGGCTATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGAc                                                                                                        <  2:888211/150‑1 (MQ=255)
                                               aGGTACATATAGAGTTGCTGCTTATCATAATCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACg                                                                                                       <  2:2638092/150‑1 (MQ=255)
                                                  tACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGAtt                                                                                                    <  1:174760/150‑1 (MQ=255)
                                                    cAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTg                                                                                                  >  2:1383451/1‑150 (MQ=255)
                                                    cAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCGTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTg                                                                                                  >  2:3803973/1‑150 (MQ=255)
                                                     aGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGa                                                                                                 >  1:2672998/1‑150 (MQ=255)
                                                     aGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGa                                                                                                 >  1:4612013/1‑150 (MQ=255)
                                                      gATAGAGTTTCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAg                                                                                                <  2:72998/150‑1 (MQ=255)
                                                      gATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAg                                                                                                <  2:1116696/150‑1 (MQ=255)
                                                         agatTTGCTGATTATAATACTCACGGTGCATATAGTCATGGCTATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGgcg                                                                                             <  2:4448144/150‑1 (MQ=255)
                                                         agagTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGgcg                                                                                             >  1:4123809/1‑150 (MQ=255)
                                                          gagTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGAGTTGAGgcgc                                                                                            >  2:2757918/1‑150 (MQ=255)
                                                            gTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTTGCAACACCTGATGGACGATTTGAGGCGCtt                                                                                          <  2:1947490/150‑1 (MQ=255)
                                                             ttGCTTATTATTATACTCACTGTGCAGATTGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATATTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTta                                                                                         <  2:3379990/150‑1 (MQ=255)
                                                                ctgctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaata                                                                                      <  1:4070463/150‑1 (MQ=255)
                                                                ctgctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaata                                                                                      >  2:2218616/1‑150 (MQ=255)
                                                                  tctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTttataat                                                                                    <  2:1996493/149‑1 (MQ=255)
                                                                  gctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataat                                                                                    <  2:2266048/150‑1 (MQ=255)
                                                                   ctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTTATGAATTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataataa                                                                                  <  2:2322159/150‑1 (MQ=255)
                                                                   ctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataata                                                                                   <  2:2742312/150‑1 (MQ=255)
                                                                   ctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataataa                                                                                  <  2:2322139/150‑1 (MQ=255)
                                                                    ttATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataataa                                                                                  <  1:4471784/150‑1 (MQ=255)
                                                                     tATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATaaa                                                                                 <  1:3874657/150‑1 (MQ=255)
                                                                     tATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATaaa                                                                                 <  1:419142/150‑1 (MQ=255)
                                                                      aTCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAg                                                                                >  2:1935926/1‑150 (MQ=255)
                                                                      aTCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAg                                                                                >  2:4102547/1‑150 (MQ=255)
                                                                          tACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGgaga                                                                            >  1:2057004/1‑150 (MQ=255)
                                                                               aCGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAatata                                                                       <  1:824415/150‑1 (MQ=255)
                                                                                    gCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATAttttt                                                                  >  2:3550070/1‑150 (MQ=255)
                                                                                     cAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTa                                                                 <  1:2277620/150‑1 (MQ=255)
                                                                                                    aTTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTcc                                                  <  2:1264607/150‑1 (MQ=255)
                                                                                                     ttATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTccc                                                 <  1:1903153/150‑1 (MQ=255)
                                                                                                      tATTTTTATTAATTTTTTTCAATTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATATCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCg                                                <  2:4572764/150‑1 (MQ=255)
                                                                                                       aTGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCgg                                               >  2:833793/1‑150 (MQ=255)
                                                                                                          ttGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGAt                                            <  1:3416957/150‑1 (MQ=255)
                                                                                                              ttaaTTTTTTTCAAGTAGATATAAAAATCTTTTAAGAATGAAATATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAg                                        <  2:3611145/148‑1 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                      >  2:1196566/1‑150 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                      <  2:1833740/150‑1 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                      >  2:141020/1‑150 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                      <  1:3652441/150‑1 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                      <  1:1308459/150‑1 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGGGTCCCGGGATTGAGtt                                      >  1:2329560/1‑150 (MQ=255)
                                                                                                                 aTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTc                                     >  1:168252/1‑150 (MQ=255)
                                                                                                                 aTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTc                                     >  2:1618550/1‑150 (MQ=255)
                                                                                                                    tttttCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGt                                  >  1:3635864/1‑150 (MQ=255)
                                                                                                                     ttttCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTa                                 >  1:2400520/1‑150 (MQ=255)
                                                                                                                          aaGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCAt                            <  1:1246619/150‑1 (MQ=255)
                                                                                                                                cGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTa                      >  1:153560/1‑150 (MQ=255)
                                                                                                                                    aaaaaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCTCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCAt                  >  1:2592407/1‑150 (MQ=255)
                                                                                                                                     aaaaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATg                 >  2:3310875/1‑150 (MQ=255)
                                                                                                                                      aaaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGa                >  2:4018846/1‑150 (MQ=255)
                                                                                                                                       aaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGAt               <  2:2128615/150‑1 (MQ=255)
                                                                                                                                        aTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGAtt              <  1:1559666/150‑1 (MQ=255)
                                                                                                                                        aTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGAtt              >  2:1396164/1‑150 (MQ=255)
                                                                                                                                         tattGTAATAATTAAAGATAATTATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGTAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTa             <  2:235549/148‑1 (MQ=255)
                                                                                                                                          cTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTaa            <  1:2191750/150‑1 (MQ=255)
                                                                                                                                           ttGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAAt           <  2:1439793/150‑1 (MQ=255)
                                                                                                                                             gTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAAttt         >  2:3119147/1‑150 (MQ=255)
                                                                                                                                              tAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAAtttt        >  2:355290/1‑150 (MQ=255)
                                                                                                                                               aagaagGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTa       <  1:1152624/150‑1 (MQ=255)
                                                                                                                                               aagaagGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTa       <  1:1152673/150‑1 (MQ=255)
                                                                                                                                                agaagTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAg      >  2:3569717/1‑150 (MQ=255)
                                                                                                                                                agaagTAAAGATAATAATCTTGAAGCATTTTATCTAGCAACACATTATTGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAg      <  1:1179863/150‑1 (MQ=255)
                                                                                                                                                 gaagTAAAGATAATAATCTTGAAGCATTTTATCTCTCAACACCTTATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATAAGTAATAGTGTCCCGGGATTGAGTTAAGTATGCATACCGTATCATGATTAATTTTAGt     <  2:1381357/150‑1 (MQ=255)
                                                                                                                                                   agTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGc   >  2:4555350/1‑150 (MQ=255)
                                                                                                                                                   agTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGc   >  2:2500778/1‑150 (MQ=255)
                                                                                                                                                    gTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCt  <  2:926021/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCT  >  CP009273/522710‑523007

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: