Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,404,288 Δ1 bp 6.2% coding (872/1209 nt) PP_2106 → ammonium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,404,2880T. 6.2% 219.8 / 10.1 80coding (872/1209 nt)PP_2106ammonium transporter
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base . (3/2);  total (41/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.24e-01

GATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAGGTGGCGTTGGGT  >  NC_002947/2404142‑2404436
                                                                                                                                                   |                                                                                                                                                     
gATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCtggt                                                                                                                                                    <  1:671232/150‑1 (MQ=255)
 aTGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCtggtg                                                                                                                                                   >  2:1675967/1‑150 (MQ=255)
     ccATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGAGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCAGGCCGGTGGCga                                                                                                                                              <  2:388567/150‑2 (MQ=255)
       aTGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCgggg                                                                                                                                            >  2:234265/1‑150 (MQ=255)
         ggtgggtgGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGgc                                                                                                                                          >  1:677076/1‑150 (MQ=255)
              gtggCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGt                                                                                                                                     >  1:47287/1‑150 (MQ=255)
              gtggCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGt                                                                                                                                     <  2:1433513/150‑1 (MQ=255)
               tggCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTATGCATCAAGGCCCGCTGGCCGGGCTGGTTTCGATTAGTTCATGGGCCTATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGtt                                                                                                                                    >  2:1161559/1‑150 (MQ=255)
                 gCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTg                                                                                                                                  >  2:1378432/1‑150 (MQ=255)
                 gCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGATCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTg                                                                                                                                  >  1:1303562/1‑150 (MQ=255)
                  cACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCtggtg                                                                                                                                                   >  1:768294/1‑133 (MQ=255)
                  cACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCtggtg                                                                                                                                                   <  2:768294/133‑1 (MQ=255)
                  cACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGt                                                                                                                                 >  1:1685174/1‑150 (MQ=255)
                  cACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGt                                                                                                                                 <  1:234265/150‑1 (MQ=255)
                   ccccTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTc                                                                                                                                <  1:1856237/149‑1 (MQ=255)
                    cccTTTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCt                                                                                                                               >  2:1268940/1‑150 (MQ=255)
                    cccTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGCGGTCTGTTTGTCt                                                                                                                               <  1:921378/150‑1 (MQ=255)
                        gTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGAc                                                                                                                           >  2:1240423/1‑150 (MQ=255)
                         tCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGGGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGAcc                                                                                                                          >  2:279273/1‑150 (MQ=255)
                              gTTGCTGGCCGGGCGCAACGACCGGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCGGGTGGCGGTTTGCGCAGGTTCCGATCTGATTCGCCCGGGGGGTGCGCTCGCGACCGGCCTGGTGGCGGGGGCGGTGTTTGTCTGGaggttaac                                                                                                                     >  1:185865/1‑143 (MQ=255)
                                       cGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCa                                                                                                            >  1:847587/1‑150 (MQ=255)
                                          gcgcAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTTGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCaaaa                                                                                                         <  1:1445262/150‑1 (MQ=255)
                                          gcgcAACGACCCAGGCTTTCTCCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCaaaa                                                                                                         <  1:883389/150‑1 (MQ=255)
                                              aaCGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCgg                                                                                                                                              >  2:173305/1‑109 (MQ=255)
                                              aaCGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCgg                                                                                                                                              <  1:173305/109‑1 (MQ=255)
                                              aaCGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGc                                                                                                     >  2:203792/1‑150 (MQ=255)
                                               aCGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCt                                                                                                    <  2:23503/150‑1 (MQ=255)
                                                 gACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTgg                                                                                                  <  2:1923574/150‑1 (MQ=255)
                                                    ccGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGaaa                                                                                               >  1:29785/1‑150 (MQ=255)
                                                     cGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCCCCCGCCACCAGGCCGGTGGCgag                                                                                                                                             <  1:27929/103‑3 (MQ=38)
                                                     cGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCCCCCGCCACCAGGCCGGTGGCgag                                                                                                                                             >  2:27929/1‑101 (MQ=38)
                                                          tttCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCAGGCCGGTGGCgag                                                                                                                                             >  1:437521/1‑96 (MQ=255)
                                                          tttCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATcgac                                                                                         <  1:1240423/150‑1 (MQ=255)
                                                          tttCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCAGGCCGGTGGCgag                                                                                                                                             <  2:437521/98‑3 (MQ=255)
                                                              tGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATcgacga                                                                                       >  2:1710045/1‑148 (MQ=255)
                                                              tGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATcgacga                                                                                       <  1:1710045/148‑1 (MQ=255)
                                                              tGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGacgacg                                                                                     >  1:1730438/1‑150 (MQ=255)
                                                              tGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGacgacg                                                                                     >  1:519633/1‑150 (MQ=255)
                                                               gCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGt                                                                                    >  1:733612/1‑150 (MQ=255)
                                                                 aTAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTcc                                                                                  >  1:1362688/1‑150 (MQ=255)
                                                                  tAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCt                                                                                 <  2:579971/150‑1 (MQ=255)
                                                                   aaCGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTg                                                                                <  2:47287/150‑1 (MQ=255)
                                                                   aaCGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTg                                                                                >  2:146359/1‑150 (MQ=255)
                                                                      ggCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGc                                                                             <  2:909226/150‑1 (MQ=255)
                                                                      ggCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGc                                                                             <  1:203792/150‑1 (MQ=255)
                                                                         ccGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTa                                                                          >  2:1874960/1‑150 (MQ=255)
                                                                         ccGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTa                                                                          >  1:1132361/1‑150 (MQ=255)
                                                                           gCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATg                                                                        <  1:1378432/150‑1 (MQ=255)
                                                                               gCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCg                                                                    <  2:677076/150‑1 (MQ=255)
                                                                                 cGGGCTGGTGGCGATCTGCGCACTTTCCTATCTGATGCACCCGGTGGGTGAGCTCGCCTCCGGCCTGTTGGCGGGCGCCTTGTTTGTCTGCATCTTCATCGCCGCTCAAAATTGCTCGAAAATCCACGAAGTCCTGGGCTTATGGCCGtt                                                                  <  2:1303562/150‑1 (MQ=255)
                                                                                     cTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCAt                                                              <  2:1767911/150‑1 (MQ=255)
                                                                                     cTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCAt                                                              <  2:1685174/150‑1 (MQ=255)
                                                                                         tggCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCGACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTc                                                          >  2:1708554/1‑150 (MQ=255)
                                                                                          ggCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAAt                                                                                             <  2:208186/114‑1 (MQ=255)
                                                                                          ggCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAAt                                                                                             >  1:208186/1‑114 (MQ=255)
                                                                                          ggCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCt                                                         <  2:1825431/150‑1 (MQ=255)
                                                                                                  gcgcAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCg                                                 >  1:295188/1‑150 (MQ=255)
                                                                                                  gcgcAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCg                                                 >  1:868016/1‑150 (MQ=255)
                                                                                                   cgcAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGcat                                                <  2:836196/150‑3 (MQ=255)
                                                                                                    gcAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGAGGGTGGACGAGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCgtg                                               >  2:267904/1‑150 (MQ=255)
                                                                                                        gTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTggg                                           <  1:1206022/150‑1 (MQ=255)
                                                                                                         ttCCGATCTGATGCACCCGGTGGGTGCGCCCGCCACCAGGCCGGTGGCga                                                                                                                                              <  1:1088656/50‑2 (MQ=18)
                                                                                                         ttCCGATCTGATGCACCCGGTGGGTGCGCCCGCCACCAGGCCGGTGGCga                                                                                                                                              >  2:1088656/1‑49 (MQ=18)
                                                                                                          tCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGt                                         <  2:519633/150‑1 (MQ=255)
                                                                                                             gATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGc                                                                                                     >  2:237017/1‑87 (MQ=255)
                                                                                                             gATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGc                                                                                                     <  1:237017/87‑1 (MQ=255)
                                                                                                               tCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCAt                                    >  1:1168840/1‑150 (MQ=255)
                                                                                                                cTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATc                                   <  1:403995/150‑1 (MQ=255)
                                                                                                                 tGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATcgac                                                                                         <  2:72017/95‑1 (MQ=255)
                                                                                                                 tGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATcgac                                                                                         >  1:72017/1‑95 (MQ=255)
                                                                                                                       aCCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCtgtg                            >  2:1826254/1‑150 (MQ=255)
                                                                                                                       aCCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCtgtg                            >  2:1664722/1‑150 (MQ=255)
                                                                                                                       aCCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCtgtg                            <  1:1874960/150‑1 (MQ=255)
                                                                                                                       aCCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCtgtg                            <  1:1675967/150‑1 (MQ=255)
                                                                                                                         ccGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGAGGGGTGGCATCGCCTTTGGc                          <  2:1168840/150‑1 (MQ=255)
                                                                                                                         ccGGTGGGTGCGATCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGc                          <  1:1340851/150‑1 (MQ=255)
                                                                                                                               ggtgCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTc                    >  2:102645/1‑150 (MQ=255)
                                                                                                                                 tgCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCgg                  >  1:1028647/1‑150 (MQ=255)
                                                                                                                                    gcTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCa               >  1:1530613/1‑150 (MQ=255)
                                                                                                                                      tCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAgg             >  2:422776/1‑150 (MQ=255)
                                                                                                                                      tCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAgg             >  2:1729640/1‑150 (MQ=255)
                                                                                                                                       cGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAGGt            <  2:1028647/150‑1 (MQ=255)
                                                                                                                                             cGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAGGTGGCGtt      <  1:1161559/150‑1 (MQ=255)
                                                                                                                                                  cTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAGGTGGCGTTGGGt  <  2:1362688/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCGCCACCGGCCTGGTGGCGGGGGCGCTGTTTGTCTGGACCTTCACCGCCGCGCAAAATCGCTGGAAAATCGACGACGTCCTGGGCGTATGGCCGTTGCATGGTCTGTGTGGCGTGTGGGGTGGCATCGCCTGTGGCGTGTTCGGCCAGGTGGCGTTGGGT  >  NC_002947/2404142‑2404436

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