Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,435,723 G→T 11.7% D599Y (GAT→TAT)  tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,435,7230GT11.7% 537.6 / 29.9 159D599Y (GAT→TAT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (55/85);  new base T (8/11);  total (63/96)
Fisher's exact test for biased strand distribution p-value = 8.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 159 total aligned reads displayed.

ACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAG  >  CP009273/4435580‑4435869
                                                                                                                                               |                                                                                                                                                  
aCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGatat                                                                                                                                                <  1:1239240/148‑1 (MQ=255)
aCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGatat                                                                                                                                                <  1:1239245/148‑1 (MQ=255)
aCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGatat                                                                                                                                                >  1:1359904/1‑148 (MQ=255)
aCGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGatat                                                                                                                                                >  2:1156966/1‑148 (MQ=255)
 cgcgCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATc                                                                                                                                               <  2:2959501/148‑1 (MQ=255)
   cgcTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCaa                                                                                                                                             >  2:2805859/1‑148 (MQ=255)
     cTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAAt                                                                                                                                           >  1:2140221/1‑148 (MQ=255)
     cTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAAt                                                                                                                                           >  1:2645304/1‑148 (MQ=255)
           gggCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCg                                                                                                                                     <  1:1956038/148‑1 (MQ=255)
           gggCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCg                                                                                                                                     >  1:870979/1‑148 (MQ=255)
           gggCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCg                                                                                                                                     >  1:831479/1‑148 (MQ=255)
             gCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGAt                                                                                                                                   >  2:497692/1‑148 (MQ=255)
             gCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGAt                                                                                                                                   <  2:353128/148‑1 (MQ=255)
             gCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGAt                                                                                                                                   >  2:497682/1‑148 (MQ=255)
              cTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATc                                                                                                                                  <  1:3068015/148‑1 (MQ=255)
              cTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATc                                                                                                                                  <  1:3067922/148‑1 (MQ=255)
              cTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATc                                                                                                                                  <  1:204037/148‑1 (MQ=255)
                tGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAg                                                                                                                                <  1:2683062/148‑1 (MQ=255)
                 ggCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGa                                                                                                                               <  1:1108213/148‑1 (MQ=255)
                     gTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGc                                                                                                                           <  2:152103/148‑1 (MQ=255)
                         aGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCggg                                                                                                                       <  1:1570302/148‑1 (MQ=255)
                            gAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGaaa                                                                                                                    <  1:2616443/148‑1 (MQ=255)
                             aaaGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAAc                                                                                                                   <  1:497692/148‑1 (MQ=255)
                             aaaGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAAc                                                                                                                   >  1:2917194/1‑148 (MQ=255)
                              aaGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACt                                                                                                                  >  2:1863194/1‑148 (MQ=255)
                               aGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTc                                                                                                                 >  2:3230157/1‑148 (MQ=255)
                                gggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCg                                                                                                                <  2:920919/148‑1 (MQ=255)
                                gggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCg                                                                                                                >  1:2991427/1‑148 (MQ=255)
                                 ggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGa                                                                                                               >  2:1789765/1‑148 (MQ=255)
                                 ggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGa                                                                                                               >  2:2264879/1‑148 (MQ=255)
                                 ggCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGa                                                                                                               <  1:1792683/148‑1 (MQ=255)
                                    tGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGt                                                                                                            <  1:1987791/148‑1 (MQ=255)
                                     ggTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTa                                                                                                           >  1:398268/1‑148 (MQ=255)
                                         aaaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTAcgcg                                                                                                       <  1:2264879/148‑1 (MQ=255)
                                         aaaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTAcgcg                                                                                                       >  2:2106466/1‑148 (MQ=255)
                                         aaaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTAcgcg                                                                                                       >  2:2106488/1‑148 (MQ=255)
                                          aaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                      <  1:2565160/148‑1 (MQ=255)
                                          aaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                      <  2:1110019/148‑1 (MQ=255)
                                          aaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                      <  1:2564402/148‑1 (MQ=255)
                                          aaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                      <  1:132925/148‑1 (MQ=255)
                                          aaGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGt                                                                                                      >  2:2830035/1‑148 (MQ=255)
                                           aGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGtt                                                                                                     <  2:398268/148‑1 (MQ=255)
                                            gTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTg                                                                                                    >  1:2740532/1‑148 (MQ=255)
                                                 cgGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCaa                                                                                               <  1:2604456/148‑1 (MQ=255)
                                                     aCGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAAttt                                                                                           <  1:2432424/148‑1 (MQ=255)
                                                     aCGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAAttt                                                                                           <  1:2432431/148‑1 (MQ=255)
                                                     aCGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAAttt                                                                                           <  1:1859332/148‑1 (MQ=255)
                                                        gTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGc                                                                                        >  1:109957/1‑148 (MQ=255)
                                                        gTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGc                                                                                        >  2:1698570/1‑148 (MQ=255)
                                                           gAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGAACGCGTTGAGCAAATTTCGCAGc                                                                                     >  2:2268950/1‑148 (MQ=255)
                                                             aGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCg                                                                                   >  1:1609169/1‑148 (MQ=255)
                                                               cgcCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGa                                                                                 >  1:2084482/1‑148 (MQ=255)
                                                                 ccACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATg                                                                               <  2:2451506/148‑1 (MQ=255)
                                                                 ccACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATg                                                                               <  2:1637524/148‑1 (MQ=255)
                                                                     aaCTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTaaa                                                                           <  1:1598348/148‑1 (MQ=255)
                                                                      actactGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAat                                                                          >  2:3243867/1‑148 (MQ=255)
                                                                       ctactGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGCTCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAata                                                                         <  1:23447/148‑1 (MQ=255)
                                                                       ctactGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAata                                                                         <  2:2122370/148‑1 (MQ=255)
                                                                          ctGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCa                                                                      <  2:1722537/148‑1 (MQ=255)
                                                                          ctGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCa                                                                      <  2:1922495/148‑1 (MQ=255)
                                                                           tGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCaa                                                                     <  2:1274545/148‑1 (MQ=255)
                                                                           tGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCaa                                                                     <  1:353711/148‑1 (MQ=255)
                                                                                 atatCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAACAAATTTCGCAGCCGGATGTAAATATCAACCTCGt                                                               >  2:594677/1‑148 (MQ=255)
                                                                                   atCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCa                                                             <  1:2512778/148‑1 (MQ=255)
                                                                                    tCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAc                                                            >  2:2611526/1‑148 (MQ=255)
                                                                                     cACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAcc                                                           <  1:2865892/148‑1 (MQ=255)
                                                                                     cACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAcc                                                           <  1:3189466/148‑1 (MQ=255)
                                                                                      aCCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCAccc                                                          >  1:709160/1‑148 (MQ=255)
                                                                                         gcgcgcgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGa                                                       <  2:109957/148‑1 (MQ=255)
                                                                                         gcgcgcgGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGa                                                       <  1:1366739/148‑1 (MQ=255)
                                                                                                gCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaa                                                <  1:1127453/148‑1 (MQ=255)
                                                                                                gCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaa                                                <  1:1127457/148‑1 (MQ=255)
                                                                                                 ccTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCaaa                                               <  2:2575112/148‑1 (MQ=255)
                                                                                                   tGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAgg                                             >  2:1359253/1‑148 (MQ=255)
                                                                                                    gcgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGc                                            <  2:2923381/148‑1 (MQ=255)
                                                                                                    gcgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGc                                            >  2:698654/1‑148 (MQ=255)
                                                                                                    gcgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGc                                            <  2:355972/148‑1 (MQ=255)
                                                                                                     cgcTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGAGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCa                                           <  1:2366771/148‑1 (MQ=255)
                                                                                                           cAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGaaa                                     >  1:3219659/1‑148 (MQ=255)
                                                                                                            aGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGaaaa                                    <  1:1758009/148‑1 (MQ=255)
                                                                                                              gAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagc                                  >  1:784483/1‑148 (MQ=255)
                                                                                                              gAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagc                                  >  1:784518/1‑148 (MQ=255)
                                                                                                              gAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGACTGCCAAAGGCAGCGAAAagc                                  >  2:2430138/1‑148 (MQ=255)
                                                                                                              gAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagc                                  >  2:77175/1‑148 (MQ=255)
                                                                                                                aCTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagcag                                <  2:1864720/148‑1 (MQ=255)
                                                                                                                 cTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAagcagc                               <  2:1609169/148‑1 (MQ=255)
                                                                                                                  tttCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAgcagca                              <  2:2461664/148‑1 (MQ=255)
                                                                                                                    tCCGTGGCGCAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACg                            <  1:1789765/148‑1 (MQ=255)
                                                                                                                    tCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACg                            <  2:2170534/148‑1 (MQ=255)
                                                                                                                     ccGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGa                           <  1:1772998/148‑1 (MQ=255)
                                                                                                                             cAGGTTCGCGTGGAATGCTATATCAAATCCACCGATCAGCTCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGt                   <  2:784622/148‑1 (MQ=255)
                                                                                                                             cAGGTTCGCGTGGAAGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGt                   <  2:784518/148‑1 (MQ=255)
                                                                                                                                   cgcgTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGa             <  2:2140221/148‑1 (MQ=255)
                                                                                                                                    gcgTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGAt            <  2:1330833/148‑1 (MQ=255)
                                                                                                                                         gAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAgg       >  1:552524/1‑148 (MQ=255)
                                                                                                                                          aaGGCTATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAggg      >  1:996915/1‑148 (MQ=255)
                                                                                                                                           aGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAgggg     <  1:1968464/148‑1 (MQ=255)
                                                                                                                                              cGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAg  >  2:3125159/1‑148 (MQ=255)
                                                                                                                                              cGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAg  >  2:3124399/1‑148 (MQ=255)
                                                                                                                                              cGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAg  <  2:2740532/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                  
ACGCGCTGCCGGGGCTTGGCGGTACAGCGAAAGGGCTGGTGAAAGTACGCGGCACGGTGGAAGCGCCACAACTACTGGCAGATATCACCGCGCGCGGCCTGCGCTGGCAGGAACTTTCCGTGGCGCAGGTTCGCGTGGAAGGCGATATCAAATCCACCGATCAGATCGCCGGGAAACTCGACGTACGCGTTGAGCAAATTTCGCAGCCGGATGTAAATATCAACCTCGTCACCCTGAATGCCAAAGGCAGCGAAAAGCAGCACGAGCTACAGTTGCGGATTCAGGGGGAG  >  CP009273/4435580‑4435869

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: ‑