Predicted mutation
evidence seq id position mutation annotation gene description
RA CP009974 2,680,048 G→A A189T (GCC→ACC)  RPPX_11900 → nitrate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0099742,680,0480GA93.2% 699.1 / 24.8 221A189T (GCC→ACC) RPPX_11900nitrate reductase
Reads supporting (aligned to +/- strand):  ref base G (6/9);  new base A (97/108);  total (103/118)
Fisher's exact test for biased strand distribution p-value = 7.90e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.75e-02

Only 100 of 221 total aligned reads displayed.

CATGCTCTGGACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCATGTGTTCG  >  CP009974/2679900‑2680195
                                                                                                                                                    |                                                                                                                                                   
cATGCTCTGGACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGcc                                                                                                                                                   >  2:1543405/1‑151 (MQ=255)
      ctGGACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTg                                                                                                                                             >  1:82669/1‑151 (MQ=255)
         gACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCg                                                                                                                                          <  1:1653248/151‑1 (MQ=255)
         gACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCg                                                                                                                                          <  2:1468554/151‑1 (MQ=255)
           cGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGa                                                                                                                                        <  2:1504871/151‑1 (MQ=255)
               gtgCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCtg                                                                                                                                    <  2:1994128/151‑1 (MQ=255)
                tgCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCtgt                                                                                                                                   <  1:890159/151‑1 (MQ=255)
                ggcgtCGGTGCCCTGGGCACGCTGCTGGCTGTCGGCGGGCCACTGGCGTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCtgt                                                                                                                                   <  1:2596586/147‑1 (MQ=255)
                   ttCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGa                                                                                                                                <  1:1317244/151‑1 (MQ=255)
                    tCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGAGCGCTTGACTGTGGCTGACCTGGACCCTGTTAGCCAAGACCGACCTGGCGGAGCTGTGGaa                                                                                                                               <  2:1192106/151‑1 (MQ=255)
                     gggTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAAc                                                                                                                              <  2:3575199/150‑1 (MQ=255)
                      ggTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAAcc                                                                                                                             >  1:78106/1‑151 (MQ=255)
                       gTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCg                                                                                                                            <  1:4115531/151‑1 (MQ=255)
                        tGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGt                                                                                                                           <  2:3825500/151‑1 (MQ=255)
                          cccTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGc                                                                                                                         <  2:3996637/151‑1 (MQ=255)
                          cccTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGc                                                                                                                         <  1:5180764/151‑1 (MQ=255)
                          cccTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGc                                                                                                                         <  1:3282931/151‑1 (MQ=255)
                           ccTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCGGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGcc                                                                                                                        >  2:4902171/1‑151 (MQ=255)
                              ggCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCggg                                                                                                                     <  1:2094768/151‑1 (MQ=255)
                              ggCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCggg                                                                                                                     <  1:24816/151‑1 (MQ=255)
                               gCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCgggg                                                                                                                    >  2:378024/1‑151 (MQ=255)
                                ccACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGa                                                                                                                   >  1:1916766/1‑151 (MQ=255)
                                ccACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGa                                                                                                                   >  2:2322084/1‑151 (MQ=255)
                                  aCCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATg                                                                                                                 <  1:4011329/151‑1 (MQ=255)
                                  aCCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCCTGCCGGTGGCTGACCTGGAACCTGTTGGCCAAGACCGACCTGGCGGAGCTGTGGAACCGGGCCGGGGATg                                                                                                                 >  2:2638287/1‑151 (MQ=255)
                                       gctggctgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCg                                                                                                            <  1:3818064/151‑1 (MQ=255)
                                         tggctgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGTTGGTTCGAt                                                                                                          >  1:2172651/1‑151 (MQ=255)
                                         tggctgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGag                                                                                                          >  2:4417684/1‑150 (MQ=255)
                                          ggctgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATg                                                                                                         >  2:2254650/1‑151 (MQ=255)
                                            ctgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTc                                                                                                       <  1:1222894/151‑1 (MQ=255)
                                             tgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCg                                                                                                      >  2:3018250/1‑151 (MQ=255)
                                             tgTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCg                                                                                                      <  1:2958972/151‑1 (MQ=255)
                                               tCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCt                                                                                                    >  1:312409/1‑151 (MQ=255)
                                                cAGCGGGCCACTGGCCTTCGTACGCAAGGTCCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGCGCTGTGGAACCGTGCCGGGGCTGGTTCGATGTCGCTg                                                                                                   >  2:4758817/1‑151 (MQ=255)
                                                 ggCGGGCCACTGGCCGTCGTACGCGCGGTGCTGCGCGCCTGGGGCACCTGGTTGTTGCTGGGCGCGTGCCTGTGGCTGCCCTGGGCCCTGTTCGCCAAGACCGTCCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTgg                                                                                                  <  2:129325/150‑1 (MQ=255)
                                                    gggCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCgg                                                                                               >  1:614615/1‑151 (MQ=255)
                                                      gCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGtt                                                                                             <  1:2399949/151‑1 (MQ=255)
                                                        cACTGGCCTTCGTACGCAAGGTGCTGCGCGGCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTgg                                                                                           <  2:4873053/151‑1 (MQ=255)
                                                        cACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTgg                                                                                           <  2:5232616/151‑1 (MQ=255)
                                                        cACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTgg                                                                                           <  2:3883178/151‑1 (MQ=255)
                                                           tGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCtt                                                                                        >  1:2282079/1‑151 (MQ=255)
                                                           tGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCtt                                                                                        <  1:5035916/151‑1 (MQ=255)
                                                             gCCGTCGTACGGAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCGAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCg                                                                                      <  2:1653633/151‑1 (MQ=255)
                                                              ccTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGa                                                                                     >  2:1525564/1‑151 (MQ=255)
                                                               cTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGat                                                                                    >  2:737084/1‑151 (MQ=255)
                                                               cTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGat                                                                                    <  2:951399/151‑1 (MQ=255)
                                                                   gTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCg                                                                                >  1:855251/1‑151 (MQ=255)
                                                                   gTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCg                                                                                >  2:5060918/1‑151 (MQ=255)
                                                                     aCGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTc                                                                              <  1:1083778/151‑1 (MQ=255)
                                                                      cGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCa                                                                             <  2:4575907/151‑1 (MQ=255)
                                                                      cGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCa                                                                             >  2:312930/1‑151 (MQ=255)
                                                                      cGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCAATGTCGCTGGCGGTTGGCTTCGATATCGTcc                                                                             >  1:4229700/1‑150 (MQ=255)
                                                                        cAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATc                                                                           <  2:3172521/151‑1 (MQ=255)
                                                                        cAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATc                                                                           <  2:1032339/151‑1 (MQ=255)
                                                                        cAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATc                                                                           <  1:1747195/151‑1 (MQ=255)
                                                                          aGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGc                                                                         <  1:3384939/151‑1 (MQ=255)
                                                                           ggTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGAGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGcc                                                                        >  2:4020608/1‑151 (MQ=255)
                                                                            gtggtgCGCGCCTCGGGCATCTGGGTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTGCGCCAAGACCGAGCTGGCGGCGCTGTGGATCCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCa                                                                       <  2:5007889/151‑1 (MQ=255)
                                                                              gctgcGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATg                                                                     >  1:2470991/1‑151 (MQ=255)
                                                                               ctgcGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGc                                                                    >  2:3316595/1‑151 (MQ=255)
                                                                                tgcGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGACCCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGcc                                                                   <  1:3137425/151‑1 (MQ=255)
                                                                                  cgcgcCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGc                                                                 <  2:863691/151‑1 (MQ=255)
                                                                                    cgcCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTg                                                               >  1:2732008/1‑151 (MQ=255)
                                                                                    cgcCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTg                                                               <  1:4641023/151‑1 (MQ=255)
                                                                                      ccTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTc                                                             <  1:3004443/151‑1 (MQ=255)
                                                                                       cTGGGGCAGCTGGTTGGTGCTGGGCGCTTGGCTGTGGCTGACCTGGACCCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTcc                                                            <  2:3783176/151‑1 (MQ=255)
                                                                                        tGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCt                                                           <  2:2202421/151‑1 (MQ=255)
                                                                                        tGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCt                                                           >  1:2818272/1‑151 (MQ=255)
                                                                                            gCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCt                                                       <  1:1105021/151‑1 (MQ=255)
                                                                                             cATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTg                                                      <  1:3598869/151‑1 (MQ=255)
                                                                                              aTCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCCTGCCGCTGTCCTGGCTGc                                                     >  1:2250734/1‑151 (MQ=255)
                                                                                              aTCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGc                                                     >  1:4846086/1‑151 (MQ=255)
                                                                                               tCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGcc                                                    <  1:2213756/151‑1 (MQ=255)
                                                                                                  ggTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCt                                                 <  2:88063/151‑1 (MQ=255)
                                                                                                     tgtggCTGGGCGCTTGCCTGTGGCTGACCTGGAACCGGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGTTGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGAt                                              <  2:1896845/151‑1 (MQ=255)
                                                                                                      gtggCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATa                                             <  1:4292841/148‑1 (MQ=255)
                                                                                                             ggCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGact                                      >  2:2902133/1‑151 (MQ=255)
                                                                                                              gcgcTGGCCTGTGGCTGACCTGGACCCTGTTCGCCAAGACCGACCTGGGGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGacta                                     <  2:1098558/151‑1 (MQ=255)
                                                                                                               cgcTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTTATCGCCATGCCGCTGTCCTGGGTGCCGCTGGTAGCTGactac                                    >  1:1304345/1‑151 (MQ=255)
                                                                                                               cgcTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGactac                                    <  2:1852853/151‑1 (MQ=255)
                                                                                                                 cTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTc                                  >  2:47109/1‑151 (MQ=255)
                                                                                                                 cTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTc                                  >  2:2975492/1‑151 (MQ=255)
                                                                                                                  ttGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCa                                 >  1:3423433/1‑151 (MQ=255)
                                                                                                                   tGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCAc                                <  2:306363/151‑1 (MQ=255)
                                                                                                                    gCCTGTGGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTCCTCACg                               >  2:1812696/1‑151 (MQ=255)
                                                                                                                       tgtgGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCtt                            <  2:3309273/151‑1 (MQ=255)
                                                                                                                         tgGCTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTg                          >  2:1795430/1‑151 (MQ=255)
                                                                                                                            cTGACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGccc                       >  1:3495336/1‑151 (MQ=255)
                                                                                                                              gACCTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGc                     >  1:277565/1‑151 (MQ=255)
                                                                                                                                 cTGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCCAt                  >  2:2755757/1‑151 (MQ=255)
                                                                                                                                  tGGAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATg                 >  1:3050315/1‑151 (MQ=255)
                                                                                                                                   ggAACCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATgg                >  1:1827870/1‑151 (MQ=255)
                                                                                                                                      aCCTGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATggcag             <  2:4656715/151‑1 (MQ=255)
                                                                                                                                         tGTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATggcaggca          <  1:1522043/151‑1 (MQ=255)
                                                                                                                                          gTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCAt         <  2:11402/151‑1 (MQ=255)
                                                                                                                                          gTTCGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCGGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCAt         <  2:3701062/151‑1 (MQ=255)
                                                                                                                                             cGCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCAtgtg      <  2:2078973/151‑1 (MQ=255)
                                                                                                                                              gCCAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGCCTACTCACGCTTTGCCCGCAATGGCAGGCATgtgt     >  2:4172514/1‑151 (MQ=255)
                                                                                                                                                cAAGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCATGTGTTc   <  1:2301470/151‑1 (MQ=255)
                                                                                                                                                 aaGACCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTTGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCATGTGTTcg  >  1:1968761/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CATGCTCTGGACGCTGTGCTTCGGTGCCCTGGCCACCCTGCTGGCTGTCAGCGGGCCACTGGCCTTCGTACGCAAGGTGCTGCGCGCCTGGGGCATCTGGTTGTTGCTGGGCGCTTGCCTGTGGCTGACCTGGAACCTGTTCGCCAAGGCCGACCTGGCGGAGCTGTGGAACCGTGCCGGGGATGGTTCGATGTCGCTGGCGGTTGGCTTCGATATCGTCATCGCCATGCCGCTGTCCTGGCTGCCGCTGATAGCGGACTACTCACGCTTTGCCCGCAATGGCAGGCATGTGTTCG  >  CP009974/2679900‑2680195

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