Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,759,323 C→T 9.9% P1261S (CCC→TCC)  rhsA → Rhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,759,3230CT 9.9% 861.3 / 57.7 264P1261S (CCC→TCC) rhsARhs family protein, putative polymorphic toxin; putative polysaccharide synthesis/export protein; putative neighboring cell growth inhibitor
Reads supporting (aligned to +/- strand):  ref base C (114/123);  new base T (14/12);  total (129/135)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 269 total aligned reads displayed.

ACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGT  >  CP009273/3759174‑3759472
                                                                                                                                                     |                                                                                                                                                     
aCCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGc                                                                                                                                                       <  2:958712/150‑1 (MQ=255)
aCCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGc                                                                                                                                                       >  2:802465/1‑150 (MQ=255)
   cGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCa                                                                                                                                                    <  1:1461975/150‑1 (MQ=255)
        cAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGc                                                                                                                                               <  1:4419607/150‑1 (MQ=255)
        cAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGc                                                                                                                                               <  2:2465551/150‑1 (MQ=255)
           gggCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCt                                                                                                                                            <  1:4477093/150‑1 (MQ=255)
           gggCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCt                                                                                                                                            <  2:1313062/150‑1 (MQ=255)
              cGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGca                                                                                                                                         <  1:4216390/150‑1 (MQ=255)
              cGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGca                                                                                                                                         <  1:1985514/150‑1 (MQ=255)
                    atCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCAcagcag                                                                                                                                   >  2:2496638/1‑150 (MQ=255)
                     tCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACagcagc                                                                                                                                  >  1:773537/1‑150 (MQ=255)
                       aCTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGTCCAAAAGCCCTGCACAGCagcaa                                                                                                                                >  1:3076326/1‑150 (MQ=255)
                       aCTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCagcaa                                                                                                                                >  2:1431146/1‑150 (MQ=255)
                         tccggttCCGATTGGGCTGAAGGGGGGATGGATTTTTTATCAGTATCCGTTGAATCGAGTTACTAATACAGATCCTCTGTGGTTAGAAGTTTTTCCTAGACCATTCCCCTGGCGAATTCCATGGCCCAAAAGCCCTGCACAGCagcaagc                                                                                                                              <  2:754061/144‑1 (MQ=38)
                          cAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCagcaagca                                                                                                                             <  2:4225900/150‑1 (MQ=255)
                           aGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAg                                                                                                                            <  1:2650640/150‑1 (MQ=255)
                             gATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAgat                                                                                                                          <  1:2576149/150‑1 (MQ=255)
                                 cGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATa                                                                                                                      >  1:3711996/1‑150 (MQ=255)
                                 cGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATa                                                                                                                      >  2:3595161/1‑150 (MQ=255)
                                   aTTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAt                                                                                                                    >  1:50303/1‑150 (MQ=255)
                                   aTTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAt                                                                                                                    >  2:4322270/1‑150 (MQ=255)
                                    ttGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAtg                                                                                                                   <  1:3014340/150‑1 (MQ=255)
                                     tGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAtgc                                                                                                                  >  2:3328807/1‑150 (MQ=255)
                                      gggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAtgct                                                                                                                 <  1:1235724/150‑1 (MQ=255)
                                      gggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAtgct                                                                                                                 <  2:162316/150‑1 (MQ=255)
                                      gggCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAAtgct                                                                                                                 <  1:921899/150‑1 (MQ=255)
                                           gAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGGTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCaaa                                                                                                            <  1:486650/150‑1 (MQ=255)
                                            aaGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCaaaa                                                                                                           <  2:4024938/150‑1 (MQ=255)
                                            aaGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCaaaa                                                                                                           >  1:2126437/1‑150 (MQ=255)
                                             aGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAg                                                                                                          <  2:1465735/150‑1 (MQ=255)
                                              ggggggATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGc                                                                                                         <  2:1715599/150‑1 (MQ=255)
                                               gggggATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCa                                                                                                        >  2:3219076/1‑150 (MQ=255)
                                                  ggATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTg                                                                                                     >  1:2628414/1‑150 (MQ=255)
                                                   gATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGa                                                                                                    <  2:608483/150‑1 (MQ=255)
                                                   gATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGa                                                                                                    >  2:499697/1‑150 (MQ=255)
                                                      ggAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACaa                                                                                                 <  2:2390165/150‑1 (MQ=255)
                                                        aaTTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACaaaa                                                                                               >  1:1093244/1‑150 (MQ=255)
                                                          tttttATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtg                                                                                             >  1:2425106/1‑150 (MQ=255)
                                                           ttttATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtgg                                                                                            >  1:2392213/1‑150 (MQ=255)
                                                           ttttATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtgg                                                                                            <  1:3397062/150‑1 (MQ=255)
                                                            tttATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtggt                                                                                           >  1:3636619/1‑150 (MQ=255)
                                                              tATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAAtggtgg                                                                                         <  2:2120093/150‑1 (MQ=255)
                                                                 cAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAAc                                                                                      >  2:2314833/1‑150 (MQ=255)
                                                                  aGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACg                                                                                     >  1:981014/1‑150 (MQ=255)
                                                                   gTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGa                                                                                    >  2:317336/1‑150 (MQ=255)
                                                                   gTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCATAGACCATTCACCATGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCATATGAAAATACTTCAAAAGCATTGACAAAAATGTGGAACGa                                                                                    >  2:1656996/1‑150 (MQ=255)
                                                                    tATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGAt                                                                                   <  1:4083142/150‑1 (MQ=255)
                                                                      tCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATAc                                                                                 <  2:237960/150‑1 (MQ=255)
                                                                       ccGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACa                                                                                <  2:3895028/150‑1 (MQ=255)
                                                                       ccGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACa                                                                                <  2:418107/150‑1 (MQ=255)
                                                                       ccGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTATTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACa                                                                                <  2:3122583/150‑1 (MQ=255)
                                                                        cGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAg                                                                               >  1:1654022/1‑150 (MQ=255)
                                                                         gTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGc                                                                              >  1:172559/1‑150 (MQ=255)
                                                                          ttGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCa                                                                             <  2:346155/150‑1 (MQ=255)
                                                                          ttGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCa                                                                             <  2:1693649/150‑1 (MQ=255)
                                                                             aaTCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATca                                                                          <  1:3117594/150‑1 (MQ=255)
                                                                              aTCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATcac                                                                         >  1:119977/1‑150 (MQ=255)
                                                                                 cAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACaaa                                                                      >  2:4202469/1‑150 (MQ=255)
                                                                                   gTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGa                                                                    >  2:1465449/1‑150 (MQ=255)
                                                                                   gTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGa                                                                    >  2:2426691/1‑150 (MQ=255)
                                                                                     tACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAat                                                                  <  1:2910010/150‑1 (MQ=255)
                                                                                       cGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAatat                                                                <  2:3783465/150‑1 (MQ=255)
                                                                                       cGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAAGTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAAAGATACAGCATCACAAAGAatat                                                                <  2:1947701/150‑1 (MQ=255)
                                                                                        gAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATAtt                                                               >  2:2866047/1‑150 (MQ=255)
                                                                                        gAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATAtt                                                               >  1:1035685/1‑150 (MQ=255)
                                                                                         aaTACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATAttt                                                              <  1:4496428/150‑1 (MQ=255)
                                                                                          aTACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTg                                                             <  1:2794445/150‑1 (MQ=255)
                                                                                              aGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGActc                                                         >  1:439366/1‑150 (MQ=255)
                                                                                                 tCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGActctct                                                      <  1:4485484/150‑1 (MQ=255)
                                                                                                     ctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAata                                                  <  1:401368/150‑1 (MQ=255)
                                                                                                     ctGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAata                                                  >  2:3017145/1‑150 (MQ=255)
                                                                                                      tGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAatat                                                 >  2:2044354/1‑150 (MQ=255)
                                                                                                      tGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAatat                                                 <  2:610161/150‑1 (MQ=255)
                                                                                                        gggTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGTATCACAAAGAATATTTGACTCTCTAATATTg                                               >  1:1435723/1‑150 (MQ=255)
                                                                                                         ggTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGa                                              >  1:1158261/1‑150 (MQ=255)
                                                                                                          gTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGaa                                             >  2:4218216/1‑150 (MQ=255)
                                                                                                            tAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGaata                                           >  1:2010009/1‑150 (MQ=255)
                                                                                                                 gTTTTTCCTAGACCATTCCCCTTGCAATTTCCATGGTCCAAAAGCCCTGCACAGCAGCATGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCg                                      <  2:179231/150‑1 (MQ=255)
                                                                                                                     ttCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGAc                                  >  2:277717/1‑150 (MQ=255)
                                                                                                                     ttCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGAc                                  >  2:4551967/1‑150 (MQ=255)
                                                                                                                      tCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACt                                 >  1:4087582/1‑150 (MQ=255)
                                                                                                                       ccTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTa                                <  2:212583/150‑1 (MQ=255)
                                                                                                                        cTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCATCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGTAAGATACAGCATCACAAAGATTATTTGACTCTCTAATATTGAATAATCCGGGACTAg                               <  1:348878/150‑1 (MQ=255)
                                                                                                                         tAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGc                              >  2:2893739/1‑150 (MQ=255)
                                                                                                                                ttCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGat                       <  1:3336135/150‑1 (MQ=255)
                                                                                                                                 tCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGata                      <  1:1554925/150‑1 (MQ=255)
                                                                                                                                      ttGCCAATTCCATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACa                 >  1:1397581/1‑150 (MQ=255)
                                                                                                                                      ttGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGTCTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAAcc                 <  2:111131/150‑2 (MQ=255)
                                                                                                                                       tGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACaa                <  1:3270290/150‑1 (MQ=255)
                                                                                                                                        gCCAATTCCATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAAt               >  1:678287/1‑150 (MQ=255)
                                                                                                                                        gCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAAt               <  1:562899/150‑1 (MQ=255)
                                                                                                                                         ccAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATg              >  2:269220/1‑150 (MQ=255)
                                                                                                                                           aaTTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGAt            <  1:3164028/150‑1 (MQ=255)
                                                                                                                                           aaTTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGAt            <  2:878872/150‑1 (MQ=255)
                                                                                                                                            aTTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATa           <  1:2292054/150‑1 (MQ=255)
                                                                                                                                                cATGGTCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCtt       <  1:3239597/150‑1 (MQ=255)
                                                                                                                                                 aTGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTtc      >  1:1866335/1‑150 (MQ=255)
                                                                                                                                                   ggCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTtctc    >  2:1287352/1‑150 (MQ=255)
                                                                                                                                                   ggCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTtctc    <  2:3711996/150‑1 (MQ=255)
                                                                                                                                                     cccAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAATAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGt  <  2:1985947/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACCCGCTGCAGGGGCGGTATATCACTCAGGATCCGATTGGGCTGAAGGGGGGATGGAATTTTTATCAGTATCCGTTGAATCCAGTTACGAATACAGATCCTCTGGGGTTAGAAGTTTTTCCTAGACCATTCCCCTTGCCAATTCCATGGCCCAAAAGCCCTGCACAGCAGCAAGCAGATGATAATGCTGCAAAAGCATTGACAAAATGGTGGAACGATACAGCATCACAAAGAATATTTGACTCTCTAATATTGAATAATCCGGGACTAGCATTAGATATAACAATGATAGCTTCTCGT  >  CP009273/3759174‑3759472

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