Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,389,233 T→G 17.6% S129R (AGT→CGT)  yhdP ← DUF3971‑AsmA2 domains protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,389,2330TG17.6% 775.5 / 187.6 327S129R (AGT→CGT) yhdPDUF3971‑AsmA2 domains protein
Reads supporting (aligned to +/- strand):  ref base T (135/134);  new base G (29/28);  total (165/162)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

Only 100 of 330 total aligned reads displayed.

CGACGTGGATCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCAT  >  CP009273/3389085‑3389382
                                                                                                                                                    |                                                                                                                                                     
 gACGTGGATCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGt                                                                                                                                                     <  1:1028059/150‑1 (MQ=255)
 gACGTGGATCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGt                                                                                                                                                     <  1:999157/150‑1 (MQ=255)
        aTCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACt                                                                                                                                              <  2:1739477/150‑1 (MQ=255)
        aTCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACt                                                                                                                                              <  1:1718015/150‑1 (MQ=255)
        aTCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACt                                                                                                                                              <  2:3100944/150‑1 (MQ=255)
             tCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTaccac                                                                                                                                         >  1:3194632/1‑150 (MQ=255)
              cAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTaccacc                                                                                                                                        >  2:2288425/1‑150 (MQ=255)
               aGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCg                                                                                                                                       >  2:2864090/1‑150 (MQ=255)
                 ccAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCt                                                                                                                                     <  1:4843432/150‑1 (MQ=255)
                   aGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTgg                                                                                                                                   >  1:4292549/1‑150 (MQ=255)
                   aGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTgg                                                                                                                                   <  2:2046892/150‑1 (MQ=255)
                   aGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTgg                                                                                                                                   >  2:3201880/1‑150 (MQ=255)
                   aGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTgg                                                                                                                                   >  2:779524/1‑150 (MQ=255)
                     gTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTg                                                                                                                                 <  1:2255531/150‑1 (MQ=255)
                       gAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGAt                                                                                                                               <  1:115619/150‑1 (MQ=255)
                        aGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATa                                                                                                                              >  2:894761/1‑150 (MQ=255)
                            gtgGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAg                                                                                                                          <  1:3381648/150‑1 (MQ=255)
                             tgGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAgt                                                                                                                         <  1:540467/150‑1 (MQ=255)
                             tgGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAgt                                                                                                                         <  2:3064251/150‑1 (MQ=255)
                              gggATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAgtg                                                                                                                        <  2:4035627/150‑1 (MQ=255)
                                gATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGtt                                                                                                                      <  1:2169803/150‑1 (MQ=255)
                                gATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGtt                                                                                                                      <  2:912812/150‑1 (MQ=255)
                                 aTCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTg                                                                                                                     >  1:1610847/1‑150 (MQ=255)
                                   cGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGt                                                                                                                   <  1:757890/150‑1 (MQ=255)
                                    gCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTg                                                                                                                  <  2:3802612/150‑1 (MQ=255)
                                      cAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCg                                                                                                                >  1:3873900/1‑150 (MQ=255)
                                          tCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAg                                                                                                            >  2:469757/1‑150 (MQ=255)
                                           cGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAgc                                                                                                           <  1:1170583/150‑1 (MQ=255)
                                           cGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAgc                                                                                                           <  1:2534838/150‑1 (MQ=255)
                                            ggCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAgcg                                                                                                          >  1:3021399/1‑150 (MQ=255)
                                              cgcgcTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCa                                                                                                        >  1:807320/1‑150 (MQ=255)
                                              cgcgcTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCa                                                                                                        >  1:807242/1‑150 (MQ=255)
                                               gcgcTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAg                                                                                                       <  2:4839563/150‑1 (MQ=255)
                                                    gACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGcc                                                                                                  >  1:1079943/1‑150 (MQ=255)
                                                     aCCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCa                                                                                                 <  1:486969/150‑1 (MQ=255)
                                                        ggATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGaa                                                                                              >  1:4567848/1‑150 (MQ=255)
                                                         gATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGaaa                                                                                             >  1:2305702/1‑150 (MQ=255)
                                                           tGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGt                                                                                           >  2:1103473/1‑150 (MQ=255)
                                                               gTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCCGAAAGTGAgg                                                                                       >  1:2251714/1‑150 (MQ=255)
                                                                 cAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTc                                                                                     <  1:3630707/150‑1 (MQ=255)
                                                                    gAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTcgcg                                                                                  >  2:3877393/1‑150 (MQ=255)
                                                                     aaaCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGa                                                                                 <  2:1377227/150‑1 (MQ=255)
                                                                     aaaCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGa                                                                                 <  1:779524/150‑1 (MQ=255)
                                                                     aaaCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGa                                                                                 <  1:4742602/150‑1 (MQ=255)
                                                                     aaaCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGa                                                                                 >  2:1814853/1‑150 (MQ=255)
                                                                      aaCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGaa                                                                                <  1:976067/150‑1 (MQ=255)
                                                                       aCTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGaaa                                                                               <  1:4076228/150‑1 (MQ=255)
                                                                        cTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGTGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAAc                                                                              >  2:3992591/1‑150 (MQ=255)
                                                                            cTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGcc                                                                          >  1:1962498/1‑150 (MQ=255)
                                                                             ttCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCa                                                                         >  2:70213/1‑150 (MQ=255)
                                                                             ttCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCa                                                                         <  2:4704545/150‑1 (MQ=255)
                                                                              tCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACAGCCAg                                                                        >  1:3102350/1‑150 (MQ=255)
                                                                                aCTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAgcg                                                                      >  1:4702788/1‑150 (MQ=255)
                                                                                 cTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAgcgc                                                                     >  1:560638/1‑150 (MQ=255)
                                                                                 cTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAgcgc                                                                     >  2:2603657/1‑150 (MQ=255)
                                                                                  tGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCa                                                                    >  2:3261931/1‑150 (MQ=255)
                                                                                        gAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTa                                                              <  1:1110919/150‑1 (MQ=255)
                                                                                        gAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTa                                                              >  2:1548018/1‑150 (MQ=255)
                                                                                        gAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTa                                                              >  2:1162002/1‑150 (MQ=255)
                                                                                         aaGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTaa                                                             >  2:1862056/1‑150 (MQ=255)
                                                                                         aaGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTaa                                                             <  2:603754/150‑1 (MQ=255)
                                                                                          aGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAAc                                                            >  2:2611301/1‑150 (MQ=255)
                                                                                          aGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAAc                                                            <  2:693620/150‑1 (MQ=255)
                                                                                           gATCGAAATTGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACa                                                           <  1:658115/150‑1 (MQ=255)
                                                                                            aTCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAg                                                          <  2:2174295/150‑1 (MQ=255)
                                                                                               gAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGct                                                       >  2:2525816/1‑150 (MQ=255)
                                                                                                aaaTGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGctc                                                      >  2:3952172/1‑150 (MQ=255)
                                                                                                aaaTGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGctc                                                      >  1:352614/1‑150 (MQ=255)
                                                                                                 aaTGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGctct                                                     <  1:1200917/150‑1 (MQ=255)
                                                                                                 aaTGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGctct                                                     <  2:3323992/150‑1 (MQ=255)
                                                                                                  aTGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTg                                                    >  1:478657/1‑150 (MQ=255)
                                                                                                       cAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGa                                               <  2:2886264/150‑1 (MQ=255)
                                                                                                       cAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGa                                               >  1:1853337/1‑150 (MQ=255)
                                                                                                          aTTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACAt                                            <  1:2890650/150‑1 (MQ=255)
                                                                                                               cGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAgc                                       <  2:4375967/150‑1 (MQ=255)
                                                                                                               cGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAgc                                       <  2:1925948/150‑1 (MQ=255)
                                                                                                                gaagaaACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAgcg                                      >  2:3887369/1‑150 (MQ=255)
                                                                                                                  agaaACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGcc                                    >  1:2712800/1‑150 (MQ=255)
                                                                                                                  agaaACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGcc                                    >  2:3638717/1‑150 (MQ=255)
                                                                                                                    aaaCAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCag                                  <  1:963188/150‑1 (MQ=255)
                                                                                                                      aCAGATCGCTGATGTGACTGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCagag                                >  2:1690555/1‑150 (MQ=255)
                                                                                                                         gATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTaa                             <  1:2070077/150‑1 (MQ=255)
                                                                                                                          aTCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAAc                            <  1:1812934/150‑1 (MQ=255)
                                                                                                                          aTCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAAc                            >  2:2444208/1‑150 (MQ=255)
                                                                                                                            cGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAAcgc                          <  1:3933972/150‑1 (MQ=255)
                                                                                                                              cTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGt                        >  2:3460630/1‑150 (MQ=255)
                                                                                                                               tGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGtt                       >  2:2803259/1‑150 (MQ=255)
                                                                                                                                    tgACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAAcc                  <  1:2440577/150‑1 (MQ=255)
                                                                                                                                     gACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCg                 <  2:1997234/150‑1 (MQ=255)
                                                                                                                                     gACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCg                 >  2:1252694/1‑150 (MQ=255)
                                                                                                                                      aCTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGa                <  2:4643585/150‑1 (MQ=255)
                                                                                                                                        tGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGaaa              <  1:3281267/150‑1 (MQ=255)
                                                                                                                                        tGGCTTCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGaaa              <  1:4130356/150‑1 (MQ=255)
                                                                                                                                         ggCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGaaaa             <  1:140191/150‑1 (MQ=255)
                                                                                                                                          gCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAAt            <  1:827212/150‑1 (MQ=255)
                                                                                                                                          gCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACGGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAAt            <  2:4528101/150‑1 (MQ=255)
                                                                                                                                            ttCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTc          <  1:3747431/150‑1 (MQ=255)
                                                                                                                                            ttCCAGACGGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTc          >  1:3892768/1‑150 (MQ=255)
                                                                                                                                                   cTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCa   <  1:4226920/150‑1 (MQ=255)
                                                                                                                                                    tGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCAt  <  2:4292384/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGACGTGGATCGTTCAGCCAGGTGAGTTGTGGGATCGCCAGCTCGGCGCGCTGACCGGATGGCGTCAGGAAACTGACTTCACTGTCGCGAAGATCGAAATGGTCAAATTGACGAAGAAACAGATCGCTGATGTGACTGGCTTCCAGACTGTCATCACTACCACCGCTGGTGATAGGAGTGTTGGTGCGAAAGCGCAGCTGCCAGAAAGTGAGGTCGCGAAACTGCCAGCGCATATGTAACAGGCTCTGCCAGACATCCAGCGCCAGAGTAACGCGTTTAACCGAAAATTCGCCGCCAT  >  CP009273/3389085‑3389382

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: