Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,660,182 G→T 7.2% intergenic (+35/‑194) metK → / → ligB methionine adenosyltransferase/DNA ligase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,660,1820GT 7.2% 205.6 / 2.8 65intergenic (+35/‑194)metK/ligBmethionine adenosyltransferase/DNA ligase B
Reads supporting (aligned to +/- strand):  ref base G (21/39);  new base T (2/3);  total (23/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCA  >  NC_002947/5660034‑5660311
                                                                                                                                                      |                                                                                                                                 
ggCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTTCTGCAGCGCCAGACGAAAGCCCCTGCCGAg                                                                                                                                   <  2:224247/149‑1 (MQ=255)
          gCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAg                                                                                                                         <  7:205217/149‑1 (MQ=255)
           ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        <  6:134665/149‑1 (MQ=255)
           ccGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        <  6:312691/149‑1 (MQ=255)
           ccGTACTCCGCAGCAGAAGACTGTCGACGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgg                                                                                                                        >  1:168823/1‑149 (MQ=255)
            cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCGCCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAggg                                                                                                                       <  1:296869/149‑1 (MQ=255)
            cGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAggg                                                                                                                       <  4:4119/149‑1 (MQ=255)
             atcgggtgcagcaGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAgggg                                                                                                                      <  8:166569/142‑1 (MQ=255)
                   cGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttt                                                                                                                >  3:101347/1‑149 (MQ=255)
                      agcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  2:133259/149‑1 (MQ=255)
                      agcagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  7:95544/149‑1 (MQ=255)
                        cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           >  6:309194/1‑149 (MQ=255)
                        cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           <  1:48306/149‑1 (MQ=255)
                        cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGc                                                                                                           <  3:241335/149‑1 (MQ=255)
                        cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                           <  5:325206/149‑7 (MQ=255)
                        cagAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACAAAAGCCCCTGCCGAGCGATCGGCAGGGGCttttgtctgg                                                                                                           <  3:7776/149‑7 (MQ=255)
                          gaagaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCTTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCt                                                                                                         <  2:81524/149‑1 (MQ=255)
                            agaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGTA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGt                                                                                                       <  2:205981/149‑1 (MQ=255)
                            agaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTg                                                                                                        <  2:321744/148‑1 (MQ=255)
                             gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                               <  7:269709/140‑1 (MQ=255)
                             gaCTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCtttttt                                                                                                               >  8:269709/1‑140 (MQ=255)
                             gaCTGGCGGCGACGTCACCTTCACCACCTTCACCTGGGATCACGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTg                                                                                                        <  1:272154/149‑1 (MQ=255)
                                tGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgtctg                                                                                                            >  4:8250/1‑134 (MQ=255)
                                tGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAGGGGCtttcgtctg                                                                                                            <  3:8250/140‑7 (MQ=255)
                                  tCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             >  7:86546/1‑137 (MQ=255)
                                  tCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTtgt                                                                                                             <  8:86546/137‑1 (MQ=255)
                                     gCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCAtt                                                                                              <  8:65050/149‑1 (MQ=255)
                                         cgacACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgggg                                                                                                                      >  7:97421/1‑121 (MQ=255)
                                         cgacACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGATCGCTCGGCAgggg                                                                                                                      <  8:97421/121‑1 (MQ=255)
                                          gacACCTTCACCACGTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCAt                                                                                         >  7:273500/1‑149 (MQ=255)
                                            caccttcaccaccttcacCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCg                                                                                       <  2:97211/149‑1 (MQ=255)
                                            caccttcaccaccttcacCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCg                                                                                       <  5:82173/149‑1 (MQ=255)
                                                    ccaccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGACCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                  >  1:248571/1‑146 (MQ=255)
                                                     caccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGt                                                                                  >  8:274478/1‑145 (MQ=255)
                                                      accttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGAt                                                                             <  5:45536/149‑1 (MQ=255)
                                                       ccttcaccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCCTCGCCCGTCAGATc                                                                            <  7:274478/149‑1 (MQ=255)
                                                           caccTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAgc                                                                        >  8:303958/1‑149 (MQ=255)
                                                                         cACCGACCGCGCCCGGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                            >  8:193408/1‑149 (MQ=255)
                                                                         cACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGtt                                                            >  5:57215/1‑149 (MQ=255)
                                                                            cgaccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCt                                                         <  6:282101/149‑1 (MQ=255)
                                                                               ccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGc                                                      <  2:168823/149‑1 (MQ=255)
                                                                               ccgCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGc                                                      <  3:203146/149‑1 (MQ=255)
                                                                                 gcgcCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTc                                                    >  2:276506/1‑149 (MQ=255)
                                                                                    cccAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTcgtc                                                 >  3:26336/1‑146 (MQ=255)
                                                                                           gTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCggg                                           <  3:226960/148‑1 (MQ=255)
                                                                                            ttGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCt                                         <  8:273500/149‑1 (MQ=255)
                                                                                               cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       <  7:193408/148‑1 (MQ=255)
                                                                                               cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       <  4:108139/148‑1 (MQ=255)
                                                                                               cGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTg                                       >  2:269817/1‑148 (MQ=255)
                                                                                                 tGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgag                                    <  1:276506/149‑1 (MQ=255)
                                                                                                  gACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGgaga                                   <  7:303958/149‑1 (MQ=255)
                                                                                                    cGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACg                                 <  6:57215/149‑1 (MQ=255)
                                                                                                      cTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAg                               <  2:38314/149‑1 (MQ=255)
                                                                                                         ccGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCAc                            >  7:298592/1‑149 (MQ=255)
                                                                                                            gCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAg                         >  7:204621/1‑149 (MQ=255)
                                                                                                                gTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGa                      >  5:83616/1‑148 (MQ=255)
                                                                                                                 tAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGaa                     <  8:298592/148‑1 (MQ=255)
                                                                                                                        ccTGCAGCTCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCAc              <  5:309194/148‑1 (MQ=255)
                                                                                                                        ccTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCa               <  8:204621/147‑1 (MQ=255)
                                                                                                                         cTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCt            <  1:269817/149‑1 (MQ=255)
                                                                                                                         cTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCt            >  6:190432/1‑149 (MQ=255)
                                                                                                                           gCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCg                                             >  6:320526/1‑114 (MQ=255)
                                                                                                                            cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa          <  4:101347/148‑1 (MQ=255)
                                                                                                                            cAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGa          >  2:292464/1‑148 (MQ=255)
                                                                                                                                  cAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgc   <  2:248571/149‑1 (MQ=255)
                                                                                                                                   aGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTgca  >  1:261539/1‑149 (MQ=255)
                                                                                                                                                      |                                                                                                                                 
GGCCACTTCGGCCGTACTCCGCAGCAGAAGACTGTCGGCGACGACACCTTCACCACCTTCACCTGGGA‑‑GCGCACCGACCGCGCCCAGTCGTTGCGTGACGCTGCCGGCCTGTAAGCTTCCTGCAGCGCCAGACGAAAGCCCCTGCCGAGCGATCGGCAGGGGCTTTTTTGTGCCTGTTACGGCATTTTCATCGCCCGTCAGATCGAGCGCCGCGCTCGTTTCTTGCTCCGTCTCGGGCTTGGAGACGAGCACCAGCGAAATATCACCTGACAACTGCA  >  NC_002947/5660034‑5660311

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: