Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,404,087 G→C 6.5% A238A (GCG→GCC)  PP_2106 → ammonium transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,404,0870GC 6.5% 199.2 / 3.8 61A238A (GCG→GCC) PP_2106ammonium transporter
Reads supporting (aligned to +/- strand):  ref base G (31/26);  new base C (2/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TGGGCGGTACCGCGATGGGCGTTTGGTGGCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCG  >  NC_002947/2403940‑2404235
                                                                                                                                                   |                                                                                                                                                    
tGGGCGGTACCGCGATGGGCGTTTGGTGGCTTTTGCGCTATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGa                                                                                                                                                     <  4:167912/149‑1 (MQ=255)
       tACCGCGATGGGCGTTTGGTGGCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCCATCAAc                                                                                                                                                >  3:95285/1‑147 (MQ=255)
               tGGGCGTTTGGTGGCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATgg                                                                                                                                      <  5:138911/149‑1 (MQ=255)
                    gTTTGGTGGCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCAtg                                                                                                                                 >  3:63648/1‑149 (MQ=255)
                         gtggCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCCATCAACGAGTTGATCGCCaccaggc                                                                                                                             <  4:95285/148‑7 (MQ=255)
                           ggCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCGCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCCGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCAtgg                                                                                                                                <  7:236200/143‑1 (MQ=255)
                           ggCTTTTGCGCCATCGAGCAGTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGCCCTTGCAGGGCGTCAGCGGCCTGGTGGCGAGCACCTCGTTGATGGCCAggg                                                                                                                                >  8:236200/1‑143 (MQ=255)
                              ttttGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCa                                                                                                                       >  4:186519/1‑149 (MQ=255)
                                    gcCATCGAGCATTCCCTTTCTGGCGCTGTGCTCGTGGATCCTCATCATCTGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCCGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATTGTGGGTGGCACCCTGt                                                                                                                 <  4:187993/149‑1 (MQ=255)
                                     ccATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTTATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGt                                                                                                                 >  2:36890/1‑148 (MQ=255)
                                        tCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGc                                                                                                             <  7:78071/149‑1 (MQ=255)
                                          gAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCCCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGt                                                                                                           <  1:167677/149‑1 (MQ=255)
                                           aGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGtt                                                                                                          >  2:2434/1‑149 (MQ=255)
                                              aTTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCt                                                                                                       <  6:130469/149‑1 (MQ=255)
                                              aTTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCt                                                                                                       >  6:55878/1‑149 (MQ=255)
                                                 cccTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTgg                                                                                                     <  7:262462/148‑1 (MQ=255)
                                                  ccTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAg                                                                                                              >  5:260627/1‑138 (MQ=255)
                                                  ccTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAg                                                                                                              <  6:260627/138‑1 (MQ=255)
                                                    tttCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCACAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCgg                                                                                                 >  8:157967/1‑149 (MQ=255)
                                                        tGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTAATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGgcg                                                                                              >  6:35042/1‑148 (MQ=255)
                                                          gcgcTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGgcg                                                                                              <  4:73129/146‑1 (MQ=255)
                                                          gcgcTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCaa                                                                                           <  3:186519/149‑1 (MQ=255)
                                                                ggCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGAccc                                                                                     >  2:171359/1‑149 (MQ=255)
                                                                     gTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCt                                                                                >  5:319274/1‑149 (MQ=255)
                                                                       ggATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCtt                                                                               >  3:285369/1‑148 (MQ=255)
                                                                          tCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTg                                                                           <  4:63648/149‑1 (MQ=255)
                                                                             tcatcaTCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCa                                                                         >  4:149002/1‑148 (MQ=255)
                                                                             tcatcaTCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCa                                                                         >  5:171981/1‑148 (MQ=255)
                                                                              catcatCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTc                                                                                                                >  1:11600/1‑108 (MQ=255)
                                                                              catcatCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTc                                                                                                                <  2:11600/108‑1 (MQ=255)
                                                                               atcatcGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATaa                                                                      <  1:171359/149‑1 (MQ=255)
                                                                                     ggCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGccc                                                                >  2:200935/1‑149 (MQ=255)
                                                                                     ggCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGccc                                                                >  6:253644/1‑149 (MQ=255)
                                                                                              ggCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCgg                                                       <  6:319274/149‑1 (MQ=255)
                                                                                              ggCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCg                                                        <  5:253644/148‑1 (MQ=255)
                                                                                                   cAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCtg                                                   <  1:2434/148‑1 (MQ=255)
                                                                                                     aCGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCtggtg                                                >  6:95081/1‑149 (MQ=255)
                                                                                                       gTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCAGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGc                                              >  2:188894/1‑149 (MQ=255)
                                                                                                          aTGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGAt                                           >  8:92135/1‑149 (MQ=255)
                                                                                                          aTGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGAt                                           <  5:159239/149‑1 (MQ=255)
                                                                                                          aTGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCAGCATAACGGCCCGCTGGCCGGGCTGGTGGCGAt                                           >  6:453/1‑149 (MQ=255)
                                                                                                                 cccAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTgcg                                      >  2:258189/1‑147 (MQ=255)
                                                                                                                 cccAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCa                                    <  5:95081/149‑1 (MQ=255)
                                                                                                                  ccAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAg                                   >  5:252954/1‑149 (MQ=255)
                                                                                                                     gACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGCCAGGTGGCGATTTGCGCAGGtt                                >  8:249759/1‑149 (MQ=255)
                                                                                                                         ttGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCg                             >  2:55542/1‑148 (MQ=255)
                                                                                                                         ttGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGa                            >  8:243701/1‑149 (MQ=255)
                                                                                                                           gCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGAt                           <  6:252954/148‑1 (MQ=255)
                                                                                                                             aGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCt                         >  3:130332/1‑148 (MQ=255)
                                                                                                                             aGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCt                         >  3:198077/1‑148 (MQ=255)
                                                                                                                                  gTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCa                   >  2:264145/1‑149 (MQ=255)
                                                                                                                                        ggCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCggt             <  6:171981/149‑1 (MQ=255)
                                                                                                                                         gCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCggtg            >  4:125242/1‑149 (MQ=255)
                                                                                                                                          ccTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCAGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCggtgg           <  5:453/149‑1 (MQ=255)
                                                                                                                                            ccatggCCATCAACGAGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCggtgggt         >  6:321727/4‑149 (MQ=255)
                                                                                                                                            ccatggCCATCAACGAGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCggtgggt         <  5:321727/146‑1 (MQ=255)
                                                                                                                                              gtggCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGAt                      >  5:2507/1‑134 (MQ=255)
                                                                                                                                              gtggCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGAt                      <  6:2507/134‑1 (MQ=255)
                                                                                                                                                 gCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCACCGGGTgcatc    <  5:55878/149‑4 (MQ=255)
                                                                                                                                                  cGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTc   <  7:157967/149‑1 (MQ=255)
                                                                                                                                                  cGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTc   <  1:200935/149‑1 (MQ=255)
                                                                                                                                                   gATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGa                       >  7:197666/1‑128 (MQ=255)
                                                                                                                                                   gATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGa                       <  8:197666/128‑1 (MQ=255)
                                                                                                                                                   gATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGCGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCg  <  1:264145/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
TGGGCGGTACCGCGATGGGCGTTTGGTGGCTTTTGCGCCATCGAGCATTCCCTTTCTGGCGCTGGGCTCGTGGATCCTCATCATCGGCTGGTTCGGCTTCAACGTCATGAGTGCCCAGACCTTGCAGGGCGTCAGCGGCCTGGTGGCGATCAACTCGTTGATGGCCATGGTGGGTGGCACCCTGTCAGCGTTGCTGGCCGGGCGCAACGACCCGGGCTTTCTGCATAACGGCCCGCTGGCCGGGCTGGTGGCGATTTGCGCAGGTTCCGATCTGATGCACCCGGTGGGTGCGCTCG  >  NC_002947/2403940‑2404235

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