Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 522,859 T→G 12.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,8590TG12.7% 592.8 / 85.8 233S1381R (AGT→AGGrhsDRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base T (106/97);  new base G (16/14);  total (122/111)
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 236 total aligned reads displayed.

GATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTT  >  CP009273/522710‑523008
                                                                                                                                                     |                                                                                                                                                     
gATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGt                                                                                                                                                       <  1:3070850/150‑1 (MQ=255)
 aTAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTa                                                                                                                                                      <  1:400064/150‑1 (MQ=255)
 aTAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTa                                                                                                                                                      >  2:2040436/1‑150 (MQ=255)
 aTAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTa                                                                                                                                                      <  2:1771588/150‑1 (MQ=255)
   aaTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTaaa                                                                                                                                                    >  2:4043694/1‑150 (MQ=255)
    aTTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAg                                                                                                                                                   >  2:1914597/1‑150 (MQ=255)
       tAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGata                                                                                                                                                >  2:3129390/1‑150 (MQ=255)
       tAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGata                                                                                                                                                >  2:4487/1‑150 (MQ=255)
        aaGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGataa                                                                                                                                               <  1:2968859/150‑1 (MQ=255)
                 gTCATATAATCTTAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTTCAGATAGTAATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGa                                                                                                                                      <  2:202916/150‑1 (MQ=255)
                  tCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGATGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTAAGTTGATGAATTTTTTTAAAGTAGCGATAAAAATATTGTAAGAAGTAAAGATAATAATCTAGaa                                                                                                                                     >  2:2749733/1‑150 (MQ=255)
                    atatCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGc                                                                                                                                   >  1:1888135/1‑150 (MQ=255)
                     tatCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCa                                                                                                                                  >  1:3941025/1‑150 (MQ=255)
                     tatCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCa                                                                                                                                  >  2:3613640/1‑150 (MQ=255)
                         ctctTAAAAGAAAATGTCCCACAGTTACAGATAGAGTTGCTGATTATCATAATCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCAtttt                                                                                                                              <  2:3558942/150‑1 (MQ=255)
                         ctctGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCAtttt                                                                                                                              <  1:2751972/150‑1 (MQ=255)
                             taaaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtct                                                                                                                          <  2:234609/149‑1 (MQ=255)
                             gaaaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtct                                                                                                                          >  1:1398007/1‑150 (MQ=255)
                             gaaaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTAtct                                                                                                                          >  1:736457/1‑150 (MQ=255)
                               aaagaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCg                                                                                                                        >  1:2957062/1‑150 (MQ=255)
                                  gaaaaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCaa                                                                                                                     >  1:283038/1‑150 (MQ=255)
                                    aaaTGTCACACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAaca                                                                                                                   <  2:1846899/150‑1 (MQ=255)
                                     aaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATTAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAAACTTGAAGCATTTTATCTCGCAaaac                                                                                                                  >  1:401482/1‑150 (MQ=255)
                                     aaTGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAacac                                                                                                                  >  1:2718640/1‑150 (MQ=255)
                                        gTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTg                                                                                                               <  1:1789379/150‑1 (MQ=255)
                                         tCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGa                                                                                                              >  1:2429307/1‑150 (MQ=255)
                                          cccACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGAt                                                                                                             >  1:4189004/1‑150 (MQ=255)
                                          cccACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGAt                                                                                                             <  2:2957134/150‑1 (MQ=255)
                                            cacaTGTAAAGATATAGTTTATTATTATCATACTCACTGTTCATATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATgg                                                                                                           <  2:3756821/150‑1 (MQ=255)
                                              caGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGAc                                                                                                         <  1:1407801/150‑1 (MQ=255)
                                              caGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGAc                                                                                                         >  1:2096801/1‑150 (MQ=255)
                                              caGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATATCGCAACACCTGATGGAc                                                                                                         >  1:4427776/1‑150 (MQ=255)
                                                ggTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGa                                                                                                       >  2:716391/1‑150 (MQ=255)
                                                  tACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGAtt                                                                                                     <  1:729710/150‑1 (MQ=255)
                                                   aCATATAGATTTTCTGATTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGAttt                                                                                                    <  2:1385071/150‑1 (MQ=255)
                                                   aCAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATGATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGAttt                                                                                                    <  1:2554127/150‑1 (MQ=255)
                                                    cAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTg                                                                                                   <  1:1980114/150‑1 (MQ=255)
                                                     aGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGa                                                                                                  <  2:1568094/150‑1 (MQ=255)
                                                       aTAGAGTTGCTGATTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAgg                                                                                                <  1:1698800/150‑1 (MQ=255)
                                                        tAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGgc                                                                                               >  2:3608157/1‑150 (MQ=255)
                                                        tAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGAATATGTAGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCTCAACACCTGATGGACGATTTGAGgc                                                                                               >  2:345774/1‑150 (MQ=255)
                                                         agagTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGgcg                                                                                              <  1:2238947/150‑1 (MQ=255)
                                                           agTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCt                                                                                            <  1:369354/150‑1 (MQ=255)
                                                           agTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCt                                                                                            <  2:1129769/150‑1 (MQ=255)
                                                           agTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCt                                                                                            >  1:2918387/1‑150 (MQ=255)
                                                            tttGCTGCTTATCATACTCTCGGTGCATATTTTCATGGAGATTATGTTGATTAATTTTTTTCAAGTAGAGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCtt                                                                                           <  2:224178/149‑1 (MQ=255)
                                                             ttGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCNCAACACCTGATGGACGATTTGAGGCGCTta                                                                                          <  1:4284230/150‑1 (MQ=255)
                                                             ttGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTta                                                                                          <  1:3096699/150‑1 (MQ=255)
                                                             ttGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTta                                                                                          <  1:3249802/150‑1 (MQ=255)
                                                               gctgctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaat                                                                                        <  1:1587612/150‑1 (MQ=255)
                                                               gctgctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaat                                                                                        >  1:1921290/1‑150 (MQ=255)
                                                               gctgctTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaat                                                                                        <  1:372822/150‑1 (MQ=255)
                                                                cttctTATCATACTCACGGTTCAGATAGTAATGGAGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaata                                                                                       <  2:3327687/147‑1 (MQ=255)
                                                                    ttATCATACTCACGGTTCATATAGTCATGGCGATTATTTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTtaataataa                                                                                   <  2:1305954/150‑1 (MQ=255)
                                                                     tATCATACTCACGGTGCATATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATaaa                                                                                  <  2:1081510/150‑1 (MQ=255)
                                                                     tATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATaaa                                                                                  <  1:4517398/150‑1 (MQ=255)
                                                                       tCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAgg                                                                                >  1:331036/1‑150 (MQ=255)
                                                                         aTACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGgag                                                                              <  1:3639996/150‑1 (MQ=255)
                                                                           aCTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGaa                                                                            >  2:30892/1‑150 (MQ=255)
                                                                                     cAGATAGTCATGGAGATTATGTTGTTGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTa                                                                  <  2:4262479/150‑1 (MQ=255)
                                                                                      aGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTAt                                                                 >  1:3968757/1‑150 (MQ=255)
                                                                                        aTAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCa                                                               >  2:1391147/1‑150 (MQ=255)
                                                                                        aTAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCa                                                               <  2:1142636/150‑1 (MQ=255)
                                                                                        aTAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCa                                                               >  1:1762016/1‑150 (MQ=255)
                                                                                        aTAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCa                                                               >  1:4091205/1‑150 (MQ=255)
                                                                                          aGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGa                                                             >  2:837680/1‑150 (MQ=255)
                                                                                            tCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTTAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGaaa                                                           >  2:2083475/1‑150 (MQ=255)
                                                                                              aTGTCGTTTATGTTTATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATa                                                         <  2:2797931/150‑1 (MQ=255)
                                                                                                ggCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAgt                                                       >  2:4168623/1‑150 (MQ=255)
                                                                                                  cGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAgtgt                                                     <  2:1875810/150‑1 (MQ=255)
                                                                                                    aTTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTcc                                                   >  1:426234/1‑150 (MQ=255)
                                                                                                     ttttGTTGATGAATTTTTTTCAAGTAGCTATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCTCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTccc                                                  <  2:3885515/147‑1 (MQ=255)
                                                                                                     ttATTTTGATTAATTTTTTTAAAGTAGCGATAAAAATATTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATATCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTccc                                                  <  2:2872900/150‑1 (MQ=255)
                                                                                                     ttATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAAGAGTGTccc                                                  >  2:3814178/1‑150 (MQ=255)
                                                                                                       aTGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCgg                                                <  1:2474298/150‑1 (MQ=255)
                                                                                                           tgatgaATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGAtt                                            <  1:3129194/150‑1 (MQ=255)
                                                                                                             atgaATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGa                                          >  2:2940998/1‑150 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                       <  1:1323972/150‑1 (MQ=255)
                                                                                                                aaTTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGtt                                       >  1:2970612/1‑150 (MQ=255)
                                                                                                                    tttttCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGt                                   <  1:801537/150‑1 (MQ=255)
                                                                                                                     ttttCAAGTAGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTa                                  >  1:478930/1‑150 (MQ=255)
                                                                                                                       ttCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATg                                >  1:1056376/1‑150 (MQ=255)
                                                                                                                         cAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCa                              <  2:1724936/150‑1 (MQ=255)
                                                                                                                          aaGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCAt                             >  2:3509135/1‑150 (MQ=255)
                                                                                                                           aGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGTGGCGCTTAATAATAAAGGAGAATATATTTTTGTCAGAGATAGTGGCCCGGGATTGAGTTCAGTATGCAtg                            >  1:2537926/1‑149 (MQ=255)
                                                                                                                             tAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATAcc                          >  2:597100/1‑150 (MQ=255)
                                                                                                                              aGCGATAAAAATCTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCg                         >  2:1422810/1‑150 (MQ=255)
                                                                                                                                cGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTa                       >  2:1509920/1‑150 (MQ=255)
                                                                                                                                 gATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTAt                      >  1:3469217/1‑150 (MQ=255)
                                                                                                                                 gATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTAt                      >  2:2338104/1‑150 (MQ=255)
                                                                                                                                    aaaaaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCAt                   <  2:2118861/150‑1 (MQ=255)
                                                                                                                                       aaTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGAt                >  1:3933796/1‑150 (MQ=255)
                                                                                                                                        aTCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGAtt               <  1:575827/150‑1 (MQ=255)
                                                                                                                                          cTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTaa             >  1:2700098/1‑150 (MQ=255)
                                                                                                                                          cTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTaa             >  1:1154811/1‑150 (MQ=255)
                                                                                                                                          cTTGTAAGAAGGAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTaa             >  2:3020172/1‑150 (MQ=255)
                                                                                                                                               aagaagTAAAGATAATAATCTTGAAGCATTGTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTa        >  1:2381375/1‑150 (MQ=255)
                                                                                                                                                agaagTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAg       >  2:4364222/1‑150 (MQ=255)
                                                                                                                                                   agTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGc    >  1:1940426/1‑150 (MQ=255)
                                                                                                                                                     tAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGGTCAGTATGCATACCGTATCATGATTAATTTTAGTGCtt  >  2:2278676/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GATAATTTAAGAAAGGAGTCATATCCTCTGAAAAGAAAATGTCCCACAGGTACAGATAGAGTTGCTGCTTATCATACTCACGGTGCAGATAGTCATGGCGATTATGTTGATGAATTTTTTTCAAGTAGCGATAAAAATCTTGTAAGAAGTAAAGATAATAATCTTGAAGCATTTTATCTCGCAACACCTGATGGACGATTTGAGGCGCTTAATAATAAAGGAGAATATATTTTTATCAGAAATAGTGTCCCGGGATTGAGTTCAGTATGCATACCGTATCATGATTAATTTTAGTGCTT  >  CP009273/522710‑523008

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: