Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,012,054 C→T 5.7% E568K (GAA→AAA)  ycbZ ← putative peptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,012,0540CT 5.7% 962.2 / 29.1 265E568K (GAA→AAA) ycbZputative peptidase
Reads supporting (aligned to +/- strand):  ref base C (116/134);  new base T (8/7);  total (124/141)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 281 total aligned reads displayed.

GGATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCCTTCTTCTACC  >  CP009273/1011906‑1012199
                                                                                                                                                    |                                                                                                                                                 
ggATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttc                                                                                                                                                   <  1:2898279/149‑1 (MQ=255)
ggATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttc                                                                                                                                                   >  1:1793414/1‑149 (MQ=255)
ggATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttc                                                                                                                                                   <  2:2430392/149‑1 (MQ=255)
ggATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttc                                                                                                                                                   <  1:1952842/149‑1 (MQ=255)
  aTTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttcct                                                                                                                                                 >  1:1323856/1‑149 (MQ=255)
  aTTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttcct                                                                                                                                                 <  2:3153574/149‑1 (MQ=255)
  aTTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGtccttcct                                                                                                                                                 >  2:1136365/1‑149 (MQ=255)
   ttCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTg                                                                                                                                                >  1:969641/1‑149 (MQ=255)
     cgcgTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGtt                                                                                                                                              <  2:3805680/149‑1 (MQ=255)
     cgcgTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGtt                                                                                                                                              <  2:3047866/149‑1 (MQ=255)
         tttATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGa                                                                                                                                          <  1:1178395/149‑1 (MQ=255)
           tATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATg                                                                                                                                        >  2:2707609/1‑149 (MQ=255)
            aTCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGc                                                                                                                                       >  2:175924/1‑149 (MQ=255)
                aCCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTgc                                                                                                                                   <  2:4254623/149‑1 (MQ=255)
                  cATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTgcgc                                                                                                                                 >  2:963170/1‑149 (MQ=255)
                        ctctGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACg                                                                                                                           >  1:2818936/1‑149 (MQ=255)
                            gTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTcc                                                                                                                       <  2:1873400/149‑1 (MQ=255)
                             tAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCt                                                                                                                      <  1:2258704/149‑1 (MQ=255)
                             tAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCt                                                                                                                      <  1:2789374/149‑1 (MQ=255)
                               aGCCTTATTTTATTGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTgg                                                                                                                    >  1:1037085/1‑149 (MQ=255)
                               aGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTgg                                                                                                                    >  1:1035950/1‑149 (MQ=255)
                               aGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTgg                                                                                                                    >  2:1943171/1‑149 (MQ=255)
                                gCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGa                                                                                                                   >  1:4296237/1‑149 (MQ=255)
                                  cTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATg                                                                                                                 <  2:21625/149‑1 (MQ=255)
                                  cTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATg                                                                                                                 >  1:1875316/1‑149 (MQ=255)
                                  cTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGAAATGGCCATGGAAAACGGTGAAGTCCTACCTGTTGCGATGCTTGCGCGATACGTTCCTGGATg                                                                                                                 >  1:4286353/1‑149 (MQ=255)
                                     attttattGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGtt                                                                                                              <  1:4543027/149‑1 (MQ=255)
                                     attttattGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGtt                                                                                                              <  1:993441/149‑1 (MQ=255)
                                      ttttattGAAGCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGttt                                                                                                             >  1:2446482/1‑149 (MQ=255)
                                      ttttattGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGttt                                                                                                             >  1:709697/1‑149 (MQ=255)
                                       tttattGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTTCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTg                                                                                                            >  2:3480234/1‑149 (MQ=255)
                                       tttattGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTg                                                                                                            >  1:4122204/1‑149 (MQ=255)
                                        ttattGAAGCACGCAGGATAGCTAACCCGTGTACGCTGAACTAGTCCGATCAGTTCGGAATAAACCAGTACAGCCAACGTCATGGCCATGGAAAACGGTGACGTACTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGc                                                                                                           <  2:1650716/149‑1 (MQ=255)
                                        ttataGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGc                                                                                                           <  2:3987105/149‑1 (MQ=255)
                                           ttGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATc                                                                                                        <  1:1487482/149‑1 (MQ=255)
                                            tGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTCATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCa                                                                                                       >  2:4130463/1‑149 (MQ=255)
                                             gAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGTATGGTTTGCATCAg                                                                                                      >  1:3880047/1‑149 (MQ=255)
                                              aaGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGc                                                                                                     <  2:3549463/149‑1 (MQ=255)
                                                gCACGCAGGATAGCTAACTCGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGcgt                                                                                                   <  1:1846355/149‑1 (MQ=255)
                                                   cGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGcgtcgt                                                                                                <  2:3850566/149‑1 (MQ=255)
                                                     cAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGttt                                                                                              <  2:3919950/149‑1 (MQ=255)
                                                     cAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGttt                                                                                              <  2:1463809/149‑1 (MQ=255)
                                                     cAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGttt                                                                                              <  1:734731/149‑1 (MQ=255)
                                                      aGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTg                                                                                             <  1:1364825/149‑1 (MQ=255)
                                                          tAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCt                                                                                         <  1:3413703/149‑1 (MQ=255)
                                                           aGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCtt                                                                                        <  1:673427/149‑1 (MQ=255)
                                                               aaCACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                                                                                    >  1:234796/1‑149 (MQ=255)
                                                               aaCACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                                                                                    <  2:15770/149‑1 (MQ=255)
                                                               aaCACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                                                                                    >  2:1112619/1‑149 (MQ=255)
                                                               aaCACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGcc                                                                                    >  1:2926124/1‑149 (MQ=255)
                                                                    gtgtACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTTCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATccc                                                                               <  1:1196039/149‑1 (MQ=255)
                                                                    gtgtACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGGCGTCGCCCTccc                                                                               >  1:91773/1‑149 (MQ=255)
                                                                    gtgtACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTGGCCGTCGCCAAccc                                                                               >  1:2081854/1‑149 (MQ=255)
                                                                       tACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCcaca                                                                            <  2:2701178/149‑1 (MQ=255)
                                                                            tGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGa                                                                       <  2:411934/149‑1 (MQ=255)
                                                                            tGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGa                                                                       >  2:4672129/1‑149 (MQ=255)
                                                                              aaCAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGAtt                                                                     <  1:2659967/149‑1 (MQ=255)
                                                                                 aaGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaa                                                                  <  1:3029613/149‑1 (MQ=255)
                                                                                  aGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTGGCCTTAGCCAACCCACACCCCAGTtaat                                                                 >  1:2739075/1‑149 (MQ=255)
                                                                                   gTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTGGGCTTCGCCCTCCCACCCCAGATTtaata                                                                >  2:3145342/1‑149 (MQ=255)
                                                                                    tCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataa                                                               >  1:4216199/1‑149 (MQ=255)
                                                                                    tCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTtaataa                                                               <  2:3271613/149‑1 (MQ=255)
                                                                                        aTCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAAc                                                           <  2:2899040/149‑1 (MQ=255)
                                                                                        aTCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAAc                                                           <  2:2688469/149‑1 (MQ=255)
                                                                                         tCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACg                                                          <  1:2432385/149‑1 (MQ=255)
                                                                                          cAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCCTCCCACACCAGATTTAATAATAAccg                                                         >  2:2015231/1‑147 (MQ=255)
                                                                                          cAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACgg                                                         <  2:1374721/149‑1 (MQ=255)
                                                                                           aGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGc                                                        <  1:3509739/149‑1 (MQ=255)
                                                                                           aGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAAAGGc                                                        >  2:2586742/1‑149 (MQ=255)
                                                                                            gTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCa                                                       <  1:2938864/149‑1 (MQ=255)
                                                                                            gTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCa                                                       <  1:848702/149‑1 (MQ=255)
                                                                                             ttCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAg                                                      >  2:718638/1‑149 (MQ=255)
                                                                                                   aTAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTGGCCTTCGCCATCCCACACCAGAATTAATAATAACGGCAGTGCgtc                                                >  2:3916936/1‑149 (MQ=255)
                                                                                                     aaaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCgtcag                                              <  1:2328706/149‑1 (MQ=255)
                                                                                                      aaCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCAACCCACACCAGATTTAATAAAAAAGGCAGTGCgtcagt                                             >  1:2552056/1‑149 (MQ=255)
                                                                                                       aCCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCgtcagtc                                            <  2:3739434/149‑1 (MQ=255)
                                                                                                       aCCAGTTCAGCCAACGTAATGGCCAGGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCgtcagtc                                            <  2:3664650/149‑1 (MQ=255)
                                                                                                       aCCAGTTCAGCCAACGTAAGGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCgtcagtc                                            <  2:1199354/149‑1 (MQ=255)
                                                                                                         cAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTcac                                          <  2:2926124/149‑1 (MQ=255)
                                                                                                          aGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTcaca                                         <  1:4087039/149‑1 (MQ=255)
                                                                                                          aGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTcaca                                         >  2:2521184/1‑149 (MQ=255)
                                                                                                            ttCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATc                                       <  2:2350497/149‑1 (MQ=255)
                                                                                                                gCCAACGTATTGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCGGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCt                                   <  1:3143437/149‑1 (MQ=255)
                                                                                                                gCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTTCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCAACCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCt                                   >  2:4591240/1‑149 (MQ=255)
                                                                                                                  cAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTAc                                 <  2:3640800/149‑1 (MQ=255)
                                                                                                                   aaCGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTTGCCATCGCCCACCCACACCCGATTTAAAAATAAACGCCGTGCGTCAGTAACATAGtaaaaa                                >  1:2955421/1‑144 (MQ=255)
                                                                                                                          tGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGa                         <  1:4578264/149‑1 (MQ=255)
                                                                                                                              cttGGAAACCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGAATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGt                     <  1:818717/147‑1 (MQ=255)
                                                                                                                                tGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACAGCCCAGATGGTgt                   >  2:1547367/1‑148 (MQ=255)
                                                                                                                                 ggAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTTGGCTTCGCCAAACCACACCAGATTTAATAATAAACGCACTGCGGAAGTCACAAAGTCTACTTCCCCGATTGTGaa                  >  1:3989917/1‑149 (MQ=255)
                                                                                                                                    aaaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAAttt               <  2:1323856/149‑1 (MQ=255)
                                                                                                                                     aaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTTGCCGTCGCCAACCCACACCAGATTTAATAAAAAACGCAGTGCGTCAGTCACAACGTCTACAGCCCACATGGTGAATTTg              >  1:2696031/1‑149 (MQ=255)
                                                                                                                                     aaCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTGTGGGCTTCGCCCACCCACACCCCAGATAATAATAACCGCCGGTCGGCAGTAACAAACCCAACAGCCCACCTGGTGAAtttt              >  1:3166611/1‑148 (MQ=255)
                                                                                                                                      aCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTTGCCGTCGCCAAACCACACCACATTTAATAAAAAAGGCAGTGCGTCAGTCACAACGTCTACAGACCACATGGTGTAtttta             >  1:1604772/1‑147 (MQ=255)
                                                                                                                                      aCGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTACTAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGc             <  2:3539964/149‑1 (MQ=255)
                                                                                                                                            aCGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTTGCCTTCGCCAACCCACACCACATTTAATAAAAAAGGCAGTGCGTCAGTCACAAACTCTAAAGCCCAGATGGTGAATTTTCcttctt       >  1:4622896/1‑149 (MQ=255)
                                                                                                                                            aCGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCcttctt       <  1:2200339/149‑1 (MQ=255)
                                                                                                                                            aCGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCcttctt       <  1:750745/149‑1 (MQ=255)
                                                                                                                                             cGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCCttcttc      <  2:1410490/149‑1 (MQ=255)
                                                                                                                                                 cttcctGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGGCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCCTTCTTCTAcc  >  2:74790/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                 
GGATTCGCGTTTATCTACCATGTTCTCTGTAAGCCTTATTTTATTGAAGCACGCAGGATAGCTAACACGTGTACGCTGAACAAGTCCGATCAGTTCGGAATAAACCAGTTCAGCCAACGTAATGGCCATGGAAAACGGTGACGTCCTTCCTGTTGCGATGCTTGCGCGATACGTTCCTGGATGGTTTGCATCAGCGTCGTTTGGCCTTCGCCATCCCACACCAGATTTAATAATAACGGCAGTGCGTCAGTCACATCGTCTACTGCCCAGATGGTGAATTTGCCTTCTTCTACC  >  CP009273/1011906‑1012199

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