Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,150,005 G→T 5.9% A146A (GCG→GCT yceG → septation protein, ampicillin sensitivity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,150,0050GT 5.9% 940.7 / 22.1 250A146A (GCG→GCTyceGseptation protein, ampicillin sensitivity
Reads supporting (aligned to +/- strand):  ref base G (121/114);  new base T (7/8);  total (128/122)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 256 total aligned reads displayed.

GAGATGCTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGC  >  CP009273/1149856‑1150153
                                                                                                                                                     |                                                                                                                                                    
gagaTGCTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAgcg                                                                                                                                                      >  1:4290844/1‑150 (MQ=255)
  gaTGCTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCa                                                                                                                                                    <  2:1192612/150‑1 (MQ=255)
  gaTGCTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCa                                                                                                                                                    <  2:3042700/150‑1 (MQ=255)
      cTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCa                                                                                                                                                <  1:2404089/150‑1 (MQ=255)
       tGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCAc                                                                                                                                               >  1:1859819/1‑150 (MQ=255)
         aaaTTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACtt                                                                                                                                             <  1:783234/150‑1 (MQ=255)
          aaTTGCTGGAAAGCTGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTg                                                                                                                                            <  2:3436057/150‑1 (MQ=255)
           aTTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGa                                                                                                                                           >  1:409712/1‑150 (MQ=255)
             tGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAAc                                                                                                                                         <  2:285567/150‑1 (MQ=255)
              gCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACt                                                                                                                                        >  1:3471146/1‑150 (MQ=255)
                  gAAAGCGGTAAAGAAGCACAGTTCCCTCTGGGACTGGTAGAAGTGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGaa                                                                                                                                    >  1:4461709/1‑150 (MQ=255)
                  gAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGaa                                                                                                                                    >  2:2219096/1‑150 (MQ=255)
                  gAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATAACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGaa                                                                                                                                    >  1:4330362/1‑150 (MQ=255)
                        ggTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAGAACCCg                                                                                                                              <  2:4322554/150‑1 (MQ=255)
                        ggTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCg                                                                                                                              <  2:2548828/150‑1 (MQ=255)
                         gTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGACAACCCgg                                                                                                                             >  2:3337201/1‑150 (MQ=255)
                          tAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGa                                                                                                                            >  2:2667236/1‑150 (MQ=255)
                            aagaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGt                                                                                                                          <  2:656290/150‑1 (MQ=255)
                             agaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTg                                                                                                                         <  1:4203920/150‑1 (MQ=255)
                             agaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTg                                                                                                                         >  1:3614007/1‑150 (MQ=255)
                             agaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTg                                                                                                                         <  1:293457/150‑1 (MQ=255)
                              gaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGCGTgg                                                                                                                        >  2:4077251/1‑150 (MQ=255)
                              gaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTgg                                                                                                                        <  2:2241369/150‑1 (MQ=255)
                              gaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTgg                                                                                                                        <  2:1739010/150‑1 (MQ=255)
                              gaagCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTgg                                                                                                                        <  1:155110/150‑1 (MQ=255)
                                agCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGAt                                                                                                                      <  1:2715647/150‑1 (MQ=255)
                                  cacaGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTg                                                                                                                    <  1:3089513/150‑1 (MQ=255)
                                   acaGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCACGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGa                                                                                                                   >  2:262840/1‑150 (MQ=255)
                                   acaGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGa                                                                                                                   >  2:626068/1‑150 (MQ=255)
                                     aGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAg                                                                                                                 >  1:1456725/1‑150 (MQ=255)
                                     aGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAg                                                                                                                 >  2:237354/1‑150 (MQ=255)
                                     aGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAg                                                                                                                 <  1:869611/150‑1 (MQ=255)
                                     aGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAg                                                                                                                 <  2:2512119/150‑1 (MQ=255)
                                      gTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAgg                                                                                                                <  2:2327566/150‑1 (MQ=255)
                                        tCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAggtt                                                                                                              <  2:4427127/150‑1 (MQ=255)
                                           ctctGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAggttggt                                                                                                           >  1:2767252/1‑150 (MQ=255)
                                           ctctGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAggttggt                                                                                                           >  1:3582034/1‑150 (MQ=255)
                                                  aCTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGcc                                                                                                    >  2:788434/1‑150 (MQ=255)
                                                  aCTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGcc                                                                                                    >  2:788453/1‑150 (MQ=255)
                                                    tGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAg                                                                                                  <  1:4487905/150‑1 (MQ=255)
                                                          aaGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCt                                                                                            >  1:3718457/1‑150 (MQ=255)
                                                          aaGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCt                                                                                            >  1:3286915/1‑150 (MQ=255)
                                                           aGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTg                                                                                           <  2:3805773/150‑1 (MQ=255)
                                                           aGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTg                                                                                           <  2:3806316/150‑1 (MQ=255)
                                                           aGGGATGCGTCTGAGCGATTACCTCAAGCAAGTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTg                                                                                           <  2:3806344/150‑1 (MQ=255)
                                                              gATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGAt                                                                                        <  1:3859952/150‑1 (MQ=255)
                                                                 gCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGta                                                                                     <  1:1984683/150‑1 (MQ=255)
                                                                  cGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGtat                                                                                    <  1:4203975/150‑1 (MQ=255)
                                                                   gTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGaata                                                                                   >  2:897938/1‑150 (MQ=255)
                                                                     cTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATAcc                                                                                 >  2:2021892/1‑150 (MQ=255)
                                                                      tGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCg                                                                                <  2:4484554/150‑1 (MQ=255)
                                                                        aGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGcc                                                                              >  1:2841779/1‑150 (MQ=255)
                                                                          cGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCaa                                                                            <  1:250538/150‑1 (MQ=255)
                                                                             ttACCTCACGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGAGTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATac                                                                         >  2:151680/1‑150 (MQ=255)
                                                                             ttACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATac                                                                         >  2:1719559/1‑150 (MQ=255)
                                                                               aCCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATacca                                                                       >  1:2844003/1‑150 (MQ=255)
                                                                               aCCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATacca                                                                       <  2:683614/150‑1 (MQ=255)
                                                                               aCCTCAAGCAATTGCGTGAGGCCCCGTAGATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATacca                                                                       <  2:683695/150‑1 (MQ=255)
                                                                                ccTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATaccac                                                                      >  2:1687388/1‑150 (MQ=255)
                                                                                ccTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATaccac                                                                      >  2:1686724/1‑150 (MQ=255)
                                                                                ccTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATaccac                                                                      <  1:3342018/150‑1 (MQ=255)
                                                                                 cTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCGAGACACATGCATGTATACCGCCAATaccacc                                                                     <  1:3687236/150‑1 (MQ=255)
                                                                                  tCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCg                                                                    <  2:4025875/150‑1 (MQ=255)
                                                                                  tCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCg                                                                    <  2:3276006/150‑1 (MQ=255)
                                                                                   cAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGTCACCTGGATGTATACCGCCAATACCACCGa                                                                   >  2:1863333/1‑150 (MQ=255)
                                                                                       cAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTc                                                               >  1:3205346/1‑150 (MQ=255)
                                                                                         aTTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTcgc                                                             >  2:560582/1‑150 (MQ=255)
                                                                                          ttGCGTGAGGCACCGTATATCAAGCATACGCTGAGCGATGATAAGTACTCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATAGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTcgcg                                                            >  1:1747552/1‑150 (MQ=255)
                                                                                               tGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACt                                                       <  1:1036043/150‑1 (MQ=255)
                                                                                                    ccccGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAagc                                                  <  1:1508418/150‑1 (MQ=255)
                                                                                                      ccGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCTGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAagcga                                                >  2:4378033/1‑150 (MQ=255)
                                                                                                      ccGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAagcga                                                >  1:215819/1‑150 (MQ=255)
                                                                                                        gTATATCAAGCATACGCTGAGCGATGATTAGTACGCCACCGTTGCGCAGGCACTTGAACTGTAACACCCGGAGTGGTTGGACGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAagcgagc                                              <  1:883380/150‑1 (MQ=255)
                                                                                                           tatCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGca                                           >  1:810647/1‑150 (MQ=255)
                                                                                                             tCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGcaca                                         >  1:4062968/1‑150 (MQ=255)
                                                                                                              cAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGAGGTATACCGCCAATACCACCGATGTCGCGTGACTCAAGCGAGCGCCCaa                                        >  2:3520012/1‑150 (MQ=255)
                                                                                                               aaGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAg                                       <  2:4296848/150‑1 (MQ=255)
                                                                                                                aGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGa                                      <  2:2607661/150‑1 (MQ=255)
                                                                                                                aGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGa                                      <  2:312474/150‑1 (MQ=255)
                                                                                                                  cATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGaaa                                    <  1:2954069/150‑1 (MQ=255)
                                                                                                                    tACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAAt                                  <  2:1064352/150‑1 (MQ=255)
                                                                                                                        cTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTg                              >  2:3719807/1‑150 (MQ=255)
                                                                                                                        cTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTg                              >  2:1550554/1‑150 (MQ=255)
                                                                                                                         tGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGa                             >  1:1380360/1‑150 (MQ=255)
                                                                                                                         tGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGAGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAACTGGTGa                             >  2:50468/1‑150 (MQ=255)
                                                                                                                            gCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAg                          <  2:3302833/150‑1 (MQ=255)
                                                                                                                               atgatAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCgg                       <  2:1409835/150‑1 (MQ=255)
                                                                                                                               atgatAAGTACGCCACCGGAGCGCAGGCACTTGAACTGGAATACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAATATGGTGAAAGCgg                       <  2:2837535/150‑1 (MQ=255)
                                                                                                                                tgatAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGt                      <  1:1040635/150‑1 (MQ=255)
                                                                                                                                 gatAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTc                     >  1:2168240/1‑150 (MQ=255)
                                                                                                                                 gatAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTc                     <  2:3934300/150‑1 (MQ=255)
                                                                                                                                 gatAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTc                     <  2:3934318/150‑1 (MQ=255)
                                                                                                                                  atAAGTACGCCACCGTAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCg                    <  1:3331218/150‑1 (MQ=255)
                                                                                                                                  atAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCg                    >  1:158676/1‑150 (MQ=255)
                                                                                                                                  atAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCg                    >  2:1898251/1‑150 (MQ=255)
                                                                                                                                   tAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGa                   >  2:1233577/1‑150 (MQ=255)
                                                                                                                                     aGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATa                 <  1:1414816/150‑1 (MQ=255)
                                                                                                                                           ccACCGCAGCGCACGCACCTGAACTGGACAACCCGGAGTGGCTTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCt           <  1:2343257/150‑1 (MQ=255)
                                                                                                                                                 tAGCTCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGTTAGCGCCTGGGAAg     >  2:4000449/1‑150 (MQ=255)
                                                                                                                                                    cgcAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGc  <  1:4034925/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GAGATGCTGAAATTGCTGGAAAGCGGTAAAGAAGCACAGTTCCCTCTGCGACTGGTAGAAGGGATGCGTCTGAGCGATTACCTCAAGCAATTGCGTGAGGCCCCGTATATCAAGCATACGCTGAGCGATGATAAGTACGCCACCGTAGCGCAGGCACTTGAACTGGAAAACCCGGAGTGGATTGAAGGTTGGTTCTGGCCAGACACCTGGATGTATACCGCCAATACCACCGATGTCGCGTTACTCAAGCGAGCGCACAAGAAAATGGTGAAAGCGGTCGATAGCGCCTGGGAAGGGC  >  CP009273/1149856‑1150153

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: