Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 725,270 T→A 5.5% Y78N (TAC→AAC)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273725,2700TA 5.5% 832.9 / 7.2 352Y78N (TAC→AAC) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base T (153/169);  new base A (4/17);  total (159/193)
Fisher's exact test for biased strand distribution p-value = 1.23e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 356 total aligned reads displayed.

TCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTT  >  CP009273/725121‑725406
                                                                                                                                                     |                                                                                                                                        
tCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCt                                                                                                                                          <  1:2642165/150‑1 (MQ=38)
 cGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTa                                                                                                                                         <  2:596695/150‑1 (MQ=38)
  ggTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAc                                                                                                                                        >  2:3350454/1‑150 (MQ=37)
  ggTGCCCCCACCGGCGTGGCCTGGTCGGTGTGCCCCGGCGGGGTGACGGCCGGCCCTCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAc                                                                                                                                        >  1:573887/1‑150 (MQ=255)
    tGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcag                                                                                                                                      <  2:894243/150‑1 (MQ=37)
    tGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcag                                                                                                                                      >  1:736729/1‑150 (MQ=37)
     gCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCACTCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGACACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagc                                                                                                                                     <  2:2424933/150‑1 (MQ=255)
      cccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagca                                                                                                                                    <  1:163313/150‑1 (MQ=37)
       ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagcag                                                                                                                                   >  1:1026216/1‑150 (MQ=37)
       ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagcag                                                                                                                                   >  2:239533/1‑150 (MQ=37)
       ccccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTAcagcag                                                                                                                                   >  1:2043047/1‑150 (MQ=37)
         ccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGtt                                                                                                                                 >  2:1872131/1‑150 (MQ=37)
         ccACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGtt                                                                                                                                 >  2:4301098/1‑150 (MQ=37)
          cACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTa                                                                                                                                <  1:3190987/150‑1 (MQ=37)
             cGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCg                                                                                                                             <  1:2301386/150‑1 (MQ=37)
             cGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCg                                                                                                                             >  2:3906232/1‑150 (MQ=37)
             cGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCg                                                                                                                             <  1:2929762/150‑1 (MQ=37)
             cGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCg                                                                                                                             <  2:2033637/150‑1 (MQ=37)
              ggCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCgg                                                                                                                            <  2:1353686/150‑1 (MQ=37)
                cGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGAc                                                                                                                          <  1:2389440/150‑1 (MQ=37)
                 gTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACa                                                                                                                         <  2:61305/150‑1 (MQ=37)
                     ccTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACaaaaa                                                                                                                     >  1:254556/1‑150 (MQ=37)
                       tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                                                                   <  1:3068701/150‑1 (MQ=37)
                       tGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                                                                   <  2:3920658/150‑1 (MQ=37)
                        gTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGc                                                                                                                  >  2:1929071/1‑150 (MQ=37)
                        gTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGc                                                                                                                  <  2:3473622/150‑1 (MQ=37)
                         ttCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGcc                                                                                                                 <  1:4666133/150‑1 (MQ=37)
                                    ccGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTgg                                                                                                      >  2:4591258/1‑150 (MQ=37)
                                      ggcggGGTGACGTCGGGCAATCCGGTAAATCCGCTGCTCGGGGCAACGGTCCTTCCCGGTGAAACCGACAGCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCGCCCGAACATACCGCATTTACGGGACAAAAACGCCCGCGCCGGTgggg                                                                                                    <  2:1451336/150‑1 (MQ=11)
                                       ggggggggACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCGTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCTGTTCCCGGACAAACACGCCCGCGCCGGTGGGGa                                                                                                   <  2:655956/148‑1 (MQ=11)
                                       gcggGGTGACGTCCGGCCAGCCGGTCAATCCGCTGCTCGGTGCAACGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTCCAGCAGTTCCCGGACAAACACGCCCGCGCCGGTGGGGa                                                                                                   <  2:4094535/150‑1 (MQ=11)
                                         ggggTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTGCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGc                                                                                                 >  1:959219/1‑150 (MQ=37)
                                            gTGCCGTCCGGCCTTCCGGTCACTCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTc                                                                                              <  2:3544363/150‑1 (MQ=11)
                                            gTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTc                                                                                              >  1:3211196/1‑150 (MQ=34)
                                             tgccGTCCGCCCATCCGGTCAAGCCGCTGCCCGTGGAAAAGGTCGTGCGCGGTTAAACCGAAACCCCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCCACCGCCTTTACCGTACAAACACGCCCGCGCCGGTGGGGAGCCTCg                                                                                             <  2:443132/147‑1 (MQ=2)
                                             tgaggcccGGCAATCCGTTCAATCCGCGGCGCGGTGCAACGGTCCGCCCCGGGCAACCCGACTTCGCCCGGCCGGGCCCGCTGCCTTACATGCACCCCCCCCCCACCAGCAGTTCCCGCACAAACACGCCCGCGCCGGTGGGGAGCCTCg                                                                                             <  2:2826651/144‑1 (MQ=2)
                                             tGACGTCCGGCCATCCGGTCAATCCGCTGCTCGCTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAGTCTCTCCCGCACCAACAGCAGTTACCGGACCAAAACGCCCGCCCCGGGGGGGAGCCACg                                                                                             >  1:3359673/1‑150 (MQ=11)
                                                   ccGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCg                                                                                       >  2:150288/1‑150 (MQ=34)
                                                   ccGGCCAGCGGTTCACGCCGCTGCGCGGTGCAACGGTCCTCCCGGGTCACCCCGACACCGCCCGGCCCGCCCCGCTGCCGTACTTTCACACCCCCACCAACCGCTTTCACCGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCg                                                                                       <  2:1580851/150‑1 (MQ=2)
                                                    cGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                                                                                      >  2:3812489/1‑150 (MQ=34)
                                                    cGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCgg                                                                                      >  2:1405809/1‑150 (MQ=34)
                                                     ggCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGc                                                                                     >  2:4018394/1‑150 (MQ=34)
                                                            cGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGaaaa                                                                              >  2:3408779/1‑150 (MQ=21)
                                                             ggTCAATCCGCTGCTGGGTGCAAAGGTCCTTCCCGGTGAACCCGACATCGCCCGGCCCGGCCGGCTGCGGTTCATTCTCTCCCGCCCCTACCGCTTTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAAt                                                                             <  2:4485171/150‑1 (MQ=11)
                                                              gTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATg                                                                            >  1:328628/1‑150 (MQ=34)
                                                                cAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGcc                                                                          >  2:2334921/1‑150 (MQ=34)
                                                                  aTCCGCTGCTGGGTGCACTGGTCCTTCCGGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCCCACCTACCGCAGTTCCCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTg                                                                        <  2:1687049/150‑1 (MQ=11)
                                                                    ccGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCg                                                                      >  2:1842905/1‑150 (MQ=34)
                                                                      gctgctCGGTGCAATGGTCCTGCCCGGTGAACCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCCCCTCCAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGa                                                                    <  2:4317156/150‑1 (MQ=11)
                                                                       ctgctgGGTGAAACGGTCCTCCCGGGTGAAACCCAAATCGCCCGGCCCGGCCCGCTGCGTTTATTTCACCCCCCCACATCCAGCTGTTCCCGTAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGat                                                                   <  1:626323/150‑1 (MQ=2)
                                                                         gctGGGTGAAACGGTCCTTCCGGGTTAAACCTACCTCGCCCGCCCGGGCCCGCGGCGTTTCTTTCTCCCCCCCACAAACCCCTTTTCCCGTAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGatat                                                                 <  1:4157872/150‑1 (MQ=2)
                                                                          ctCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATc                                                                >  1:441699/1‑150 (MQ=32)
                                                                           tgggTGCAATGGTCCTTCCGGGTGAAACCGACATCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                                               <  1:364272/148‑1 (MQ=11)
                                                                           tCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                                               >  2:2946946/1‑150 (MQ=32)
                                                                           ggggggCAATGGTCCTTCCGGGTGAACCCGACATCGCCCTGCCCGGCCCGCTGCCGTACATTCACTCCCGCACCTCCACCTGTTACCGCAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATcc                                                               <  1:1775287/145‑1 (MQ=11)
                                                                            cGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCg                                                              >  2:153493/1‑150 (MQ=32)
                                                                             ggTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc                                                             >  2:2309024/1‑150 (MQ=32)
                                                                               tGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCATTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCtt                                                           <  2:119995/150‑1 (MQ=12)
                                                                                  aaaGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa                                                        >  1:3927247/1‑150 (MQ=32)
                                                                                  aaaGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACa                                                        >  1:345974/1‑150 (MQ=32)
                                                                                   acgggccTGCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACCGCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg                                                       <  2:3770726/145‑1 (MQ=255)
                                                                                   aaGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg                                                       >  1:3991272/1‑150 (MQ=17)
                                                                                   aaGGTCAGTCCCGGTGAAACCGACCTCGCCCGGCCCGCCCCGCTGCCTTTCAGTCTCTCCCGCACCTACCGCTGTTACCGGACAACAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAg                                                       <  2:168893/150‑1 (MQ=11)
                                                                                    aGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGc                                                      >  1:1732137/1‑150 (MQ=32)
                                                                                        ccTTCCCGGTGAACCCGACATCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTCCCGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                                  <  2:4642818/150‑1 (MQ=255)
                                                                                        ccTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                                  <  2:2263743/150‑1 (MQ=17)
                                                                                        ccTTCCCGGTGAAACCGACAGCGCCCGGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcg                                                  <  2:1429144/150‑1 (MQ=12)
                                                                                         cTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTgcgc                                                 >  1:1454066/1‑150 (MQ=32)
                                                                                           tCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGCGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGa                                               >  1:3703826/1‑150 (MQ=17)
                                                                                            ccgggTGAAACCGACATCGCCCTGCCCGGCCCGCGGCCTTTCATTCTCTCCCGCACCTCCCGCTGTTCCCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGAt                                              <  2:539515/147‑1 (MQ=11)
                                                                                             ccGGTGAAACCGACATCGCCCGGCCCGGCCCGCTGCCGTACATTCTCTCCCGCACCTACCGCAGTTACCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATa                                             <  2:222310/150‑1 (MQ=255)
                                                                                               ggTGAACCCGACAGCGCCGGGCCCGGCCCGCGGCCGTTCATTCTCGCCCGCACCTCCCGCTGTTACCGGAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                           <  2:3650596/150‑1 (MQ=11)
                                                                                               ggTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                           >  1:2474256/1‑150 (MQ=32)
                                                                                               ggTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                           >  2:1087196/1‑150 (MQ=32)
                                                                                               ggTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac                                           >  1:98802/1‑150 (MQ=32)
                                                                                                 tGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGAGATCCGCTTACAGCTGCGCGATAacac                                         >  2:3456210/1‑150 (MQ=17)
                                                                                                   acaCCGACAGCCCCCGGCCGGCCCCGCGGCCTTTAATTCTCCCCCCCACCCACCGCTGTCCCCGCAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacacac                                       <  1:1235727/148‑1 (MQ=11)
                                                                                                   aaaCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAacacac                                       >  2:3447410/1‑150 (MQ=32)
                                                                                                    aaCCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACt                                      >  1:1892089/1‑150 (MQ=32)
                                                                                                       tgACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAt                                   >  2:1645740/2‑150 (MQ=17)
                                                                                                       tgACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGAt                                   >  2:1610179/2‑150 (MQ=17)
                                                                                                         aCATCGCCGGGCCCGCCCCGCTGCTTTAATTTATCTCCCGCACCTACACCAGTTACGGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATAc                                 <  1:304446/150‑1 (MQ=11)
                                                                                                         aCATCGCCCTGCCGGGCCCGCTGCCGTTAATTATCTCCCGCACCTACAGATGTTACCGCAAAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATAc                                 <  1:2676805/150‑1 (MQ=11)
                                                                                                          cagcGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACt                                <  2:3004528/147‑1 (MQ=17)
                                                                                                          cATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCAACAGCAGTTACCGGACCAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCCGCGGGAAAAAGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACAGATACt                                >  1:3599879/1‑150 (MQ=11)
                                                                                                              gCCCTGTCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGt                            >  2:681866/1‑150 (MQ=17)
                                                                                                              gCCCTGTCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGt                            >  2:3051190/1‑150 (MQ=17)
                                                                                                               cccTGCCGGGCCCGATCCCGTTCTTTCTCTCCCGCACCCACCGAAGTTACCGTACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCACGGCTTGAAAATGCCTGCGTATATCCGCTTACAGCTGCTCGATAACACACTGATACTCAGTg                           <  2:4719346/150‑1 (MQ=11)
                                                                                                                ccTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGa                          >  1:1400759/1‑150 (MQ=17)
                                                                                                                     ccggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGAAAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcg                     <  2:1792503/150‑1 (MQ=17)
                                                                                                                     ccggcccgCTGCCGTTCATTCTCGCCCGCACCTACAGCTGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcg                     <  2:4178519/150‑1 (MQ=17)
                                                                                                                      cggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcgg                    >  1:1082107/1‑150 (MQ=32)
                                                                                                                      cggcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcgg                    >  1:3545231/1‑150 (MQ=32)
                                                                                                                        gcccgCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAAcggcg                  <  2:3137310/150‑1 (MQ=17)
                                                                                                                           cgCTGCGTTTAATTCTCCCCCCCACATACCGCTTTCCCCGTACAAACCCGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCa               <  1:4525607/150‑1 (MQ=11)
                                                                                                                           cgCTGCCTTTCATTCTCTCCCGCCCCTCCAGCAGTTACCTGACAAAAACGCCCGCGCCGGTGGTGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCa               <  2:2728044/150‑1 (MQ=11)
                                                                                                                           cgCTGCCGTTCATTCTCACCCGCACCTACAGCATTTACCGGACAAACACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCa               <  2:3729715/150‑1 (MQ=12)
                                                                                                                             cTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGa             >  1:4642817/1‑150 (MQ=32)
                                                                                                                             cTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGa             >  2:946800/1‑150 (MQ=32)
                                                                                                                                     ctttCTCTCCCGCACCTACAGCTGTTACCGCACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTa     <  1:1134832/148‑1 (MQ=9)
                                                                                                                                                     |                                                                                                                                        
TCGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGTGACGTCCGGCCATCCGGTCAATCCGCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATTCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAACACACTGATACTCAGTGATAACGGCGGCAGAAGCCTGTATTT  >  CP009273/725121‑725406

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