Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 554,868 G→T 10.5% E147* (GAA→TAA)  sfmC → putative periplasmic pilus chaperone

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273554,8680GT10.5% 754.2 / 43.0 228E147* (GAA→TAA) sfmCputative periplasmic pilus chaperone
Reads supporting (aligned to +/- strand):  ref base G (106/98);  new base T (17/7);  total (123/105)
Fisher's exact test for biased strand distribution p-value = 8.75e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

Only 100 of 229 total aligned reads displayed.

GCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCAGAG  >  CP009273/554723‑555015
                                                                                                                                                 |                                                                                                                                                   
gCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGaa                                                                                                                                                   >  2:891041/1‑148 (MQ=255)
   ggCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTa                                                                                                                                                <  2:3439385/148‑1 (MQ=255)
      agcagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTAccc                                                                                                                                             >  2:237499/1‑148 (MQ=255)
        cagATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGc                                                                                                                                           >  1:1791095/1‑148 (MQ=255)
             cGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAg                                                                                                                                      <  2:1791095/148‑1 (MQ=255)
              gTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAgc                                                                                                                                     >  1:3086966/1‑148 (MQ=255)
                  gTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCg                                                                                                                                 >  1:2296727/1‑148 (MQ=255)
                    ctctGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGa                                                                                                                               <  2:779827/148‑1 (MQ=255)
                       tGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGacac                                                                                                                            <  2:2492906/148‑1 (MQ=255)
                            tGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCa                                                                                                                       <  2:2127615/148‑1 (MQ=255)
                                 gAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGttt                                                                                                                  >  2:1868819/1‑148 (MQ=255)
                                 gAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGttt                                                                                                                  <  1:411869/148‑1 (MQ=255)
                                 gAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGttt                                                                                                                  >  2:986639/1‑148 (MQ=255)
                                    tGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTcg                                                                                                               >  1:1401542/1‑148 (MQ=255)
                                    tGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTcg                                                                                                               <  2:3341008/148‑1 (MQ=255)
                                        aaGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGAt                                                                                                           <  1:777144/148‑1 (MQ=255)
                                         aGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATc                                                                                                          <  1:46428/148‑1 (MQ=255)
                                          gACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATcc                                                                                                         <  1:2600448/148‑1 (MQ=255)
                                           aCGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCg                                                                                                        <  1:560454/148‑1 (MQ=255)
                                            cGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCgg                                                                                                       >  2:199397/1‑148 (MQ=255)
                                             gATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGt                                                                                                      >  2:675802/1‑148 (MQ=255)
                                               tCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTaa                                                                                                    >  2:3276530/1‑148 (MQ=255)
                                                cccTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAAc                                                                                                   <  1:1971269/148‑1 (MQ=255)
                                                 ccTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACt                                                                                                  <  1:2619578/148‑1 (MQ=255)
                                                 ccTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACt                                                                                                  >  2:1138899/1‑148 (MQ=255)
                                                  cTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACta                                                                                                 >  2:2849370/1‑148 (MQ=255)
                                                   ttCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACtat                                                                                                <  1:3548501/148‑1 (MQ=255)
                                                    tCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACtata                                                                                               <  2:54882/148‑1 (MQ=255)
                                                        tAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCaa                                                                                           >  2:2290572/1‑148 (MQ=255)
                                                         aGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAAt                                                                                          <  1:2341453/148‑1 (MQ=255)
                                                         aGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAAt                                                                                          >  2:3101566/1‑148 (MQ=255)
                                                          gATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATg                                                                                         >  2:2175499/1‑148 (MQ=255)
                                                          gATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATg                                                                                         <  2:300964/148‑1 (MQ=255)
                                                          gATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATg                                                                                         >  2:1830404/1‑148 (MQ=255)
                                                            tAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGtt                                                                                       >  1:101855/1‑148 (MQ=255)
                                                                 aTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATaa                                                                                  >  1:509407/1‑148 (MQ=255)
                                                                 aTGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCAGTCCAATTAAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATaa                                                                                  <  1:602681/148‑1 (MQ=255)
                                                                   gCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATc                                                                                <  1:3306711/148‑1 (MQ=255)
                                                                    cATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATcc                                                                               <  1:2511261/148‑1 (MQ=255)
                                                                     aTTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCa                                                                              >  1:2946066/1‑148 (MQ=255)
                                                                     aTTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCa                                                                              >  1:2096890/1‑148 (MQ=255)
                                                                      ttGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCAt                                                                             <  1:1716916/148‑1 (MQ=255)
                                                                        gAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCa                                                                           >  2:152871/1‑148 (MQ=255)
                                                                        gAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCa                                                                           >  2:2693251/1‑148 (MQ=255)
                                                                            ggCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCtt                                                                       <  1:987136/148‑1 (MQ=255)
                                                                                ggAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTa                                                                   <  2:1915017/148‑1 (MQ=255)
                                                                                  aaTGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATg                                                                 >  1:12725/1‑148 (MQ=255)
                                                                                  aaTGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGGCCAATACAATTACAAGAATAACCCCCCGAAGCGCCGGACACACGCAAGTTTTCCCGAGCCGGTAACTATAACAATGTGCAAAAACCATCACCTTTTTATg                                                                 >  1:1956956/1‑148 (MQ=255)
                                                                                     gTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCa                                                              >  1:1094229/1‑148 (MQ=255)
                                                                                     gTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCa                                                              <  1:1120252/148‑1 (MQ=255)
                                                                                     gTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCa                                                              <  1:2300260/148‑1 (MQ=255)
                                                                                      ttttGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCAc                                                             >  2:1841610/1‑148 (MQ=255)
                                                                                       tttGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCAcc                                                            >  2:2848677/1‑148 (MQ=255)
                                                                                       tttGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCAcc                                                            >  2:1987908/1‑148 (MQ=255)
                                                                                        ttGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCAccc                                                           <  2:1218504/148‑1 (MQ=255)
                                                                                          gCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTg                                                         >  2:2897946/1‑148 (MQ=255)
                                                                                          gCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTg                                                         >  2:1509639/1‑148 (MQ=255)
                                                                                             aCTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGtt                                                      >  1:1244283/1‑148 (MQ=255)
                                                                                               tGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTaa                                                    >  2:2053678/1‑148 (MQ=255)
                                                                                                ggCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAAc                                                   >  1:162080/1‑148 (MQ=255)
                                                                                                ggCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAAc                                                   <  1:2761732/148‑1 (MQ=255)
                                                                                                ggCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAAc                                                   >  1:3224720/1‑148 (MQ=255)
                                                                                                 gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  <  1:882146/148‑1 (MQ=255)
                                                                                                 gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  >  2:3366969/1‑148 (MQ=255)
                                                                                                 gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  >  1:1521711/1‑148 (MQ=255)
                                                                                                 gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  <  2:1834760/148‑1 (MQ=255)
                                                                                                 gCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACt                                                  <  1:986639/148‑1 (MQ=255)
                                                                                                  cGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACtt                                                 >  2:1388482/1‑148 (MQ=255)
                                                                                                  cGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACtt                                                 <  2:1430570/148‑1 (MQ=255)
                                                                                                  cGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACtt                                                 <  1:299415/148‑1 (MQ=255)
                                                                                                     ttttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACa                                              >  2:272088/1‑148 (MQ=255)
                                                                                                      tttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACaa                                             <  2:936757/148‑1 (MQ=255)
                                                                                                       ttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAg                                            <  1:1185186/148‑1 (MQ=255)
                                                                                                       ttATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAg                                            <  1:3385475/148‑1 (MQ=255)
                                                                                                         aTCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTg                                          <  1:1496781/148‑1 (MQ=255)
                                                                                                          tCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTgg                                         <  1:2351958/148‑1 (MQ=255)
                                                                                                           cgcgCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTggg                                        <  2:12725/148‑1 (MQ=255)
                                                                                                            gcgcATGAAATTATTTCTCCGTGCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGc                                       <  1:66930/148‑1 (MQ=255)
                                                                                                            gcgcATGAAATTATTTCTCCGTCCAATTCAATTACAATAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGc                                       <  1:2276359/148‑1 (MQ=255)
                                                                                                                aTGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGcc                                   >  1:1775402/1‑148 (MQ=255)
                                                                                                                  gAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaa                                 >  2:1464126/1‑148 (MQ=255)
                                                                                                                  gAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaa                                 >  1:3424682/1‑148 (MQ=255)
                                                                                                                    aaTTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCaaaa                               >  2:1644154/1‑148 (MQ=255)
                                                                                                                     attattTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAg                              >  1:3176319/1‑148 (MQ=255)
                                                                                                                        attTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTg                           >  2:2336033/1‑148 (MQ=255)
                                                                                                                         tttCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTgg                          >  1:2854022/1‑148 (MQ=255)
                                                                                                                          ttCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATAAATGTTCATAATCCATCACCTATTTATGTCACCATCGTTAACTTACAAGTCTGCAGCCAAAAGTTggg                         >  1:1871405/1‑148 (MQ=255)
                                                                                                                           tctcCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTgggg                        <  1:838833/148‑1 (MQ=255)
                                                                                                                            ctcCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGa                       >  2:2113007/1‑148 (MQ=255)
                                                                                                                             tcCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGaa                      <  1:1138911/148‑1 (MQ=255)
                                                                                                                              ccGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAAt                     >  2:2139145/1‑148 (MQ=255)
                                                                                                                                gTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGc                   >  2:2550811/1‑148 (MQ=255)
                                                                                                                                  ccAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTa                 >  2:1073560/1‑148 (MQ=255)
                                                                                                                                  ccAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTa                 <  2:105394/148‑1 (MQ=255)
                                                                                                                                   caattcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTAt                <  1:2158877/148‑1 (MQ=255)
                                                                                                                                      ttcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGc             >  1:2937136/1‑148 (MQ=255)
                                                                                                                                      ttcaattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGc             >  2:1682548/1‑148 (MQ=255)
                                                                                                                                        caattACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTg           >  1:177000/1‑148 (MQ=255)
                                                                                                                                            tACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCAcc       >  1:2717590/1‑148 (MQ=255)
                                                                                                                                                 gAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCagag  >  1:325566/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCTGGCAGCAGATCGTGAGTCTCTGTTCTGGATGAATGTTAAGACGATCCCTTCGGTAGATAAAAATGCATTGAACGGCAGGAATGTTTTGCAACTGGCGATTTTATCGCGCATGAAATTATTTCTCCGTCCAATTCAATTACAAGAATTACCCGCAGAAGCGCCGGACACACTCAAGTTTTCGCGATCCGGTAACTATATCAATGTTCATAATCCATCACCTTTTTATGTCACCCTGGTTAACTTACAAGTGGGCAGCCAAAAGTTGGGGAATGCTATGGCTGCACCCAGAG  >  CP009273/554723‑555015

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: