Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,547,815 G→A 6.1% L419L (CTG→CTA)  gntT → gluconate transporter, high‑affinity GNT I system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,547,8150GA 6.1% 203.7 / 2.4 66L419L (CTG→CTA) gntTgluconate transporter, high‑affinity GNT I system
Reads supporting (aligned to +/- strand):  ref base G (39/22);  new base A (2/2);  total (42/24)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GGTGTTAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAG  >  NC_000913/3547675‑3547960
                                                                                                                                            |                                                                                                                                                 
ggTGTTAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCa                                                                                                                                            <  8:41613/148‑1 (MQ=255)
 gtgtTAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCAt                                                                                                                                           >  2:383257/1‑148 (MQ=255)
     tAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAtc                                                                                                                                       >  7:269953/1‑148 (MQ=255)
     tAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACAATCAAATCCTGGTCGATGCTGGAAACCATTAt                                                                                                                                        >  4:303558/1‑147 (MQ=255)
         ccGGAGCTGATGGTTATTGCGGGTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCgg                                                                                                                                   <  8:96799/148‑1 (MQ=255)
                  aTGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAtc                                                                                                                                       >  3:97810/1‑135 (MQ=255)
                  aTGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTAtc                                                                                                                                       <  4:97810/135‑1 (MQ=255)
                     gTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTgg                                                                                                                       >  1:6325/1‑148 (MQ=255)
                      ttATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGt                                                                                                                      <  6:52560/148‑1 (MQ=255)
                      ttATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGt                                                                                                                      <  2:6325/148‑1 (MQ=255)
                       tATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTgg                                                                                                                       >  6:332332/1‑146 (MQ=255)
                         ttGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAgg                                                                                                                   >  4:16389/1‑148 (MQ=255)
                         ttGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAg                                                                                                                    <  5:332332/147‑1 (MQ=255)
                           gCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTAGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCt                                                                                                                 >  3:19668/1‑148 (MQ=255)
                               ttggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAATAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCGTGGTCGATGCTGGAAACCATTTTCTCGGTGTGCGGTGTGGTAGGgtgtat                                                                                                             >  8:379233/1‑146 (MQ=255)
                                tggttCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGTGCTTCTGTCTTTTCAAAGAGTACTTTAACCTGAGTATCGGCGAGACCATCAAATCCTTGTCTATGCTTGGAAACATTATCTCGGTGTTCGGTCTGGTAGGATGTctg                                                                                                            >  6:264067/1‑148 (MQ=255)
                                     ccGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTAGAAACCATTATCTCGGTGTGCGGTCTGGt                                                                                                                      >  5:114418/1‑133 (MQ=255)
                                     ccGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTAGAAACCATTATCTCGGTGTGCGGTCTGGt                                                                                                                      <  6:114418/133‑1 (MQ=255)
                                      cGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTctgctgct                                                                                                       >  3:209893/1‑147 (MQ=255)
                                       ggTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGa                                                                                                     <  2:50208/148‑1 (MQ=255)
                                         tAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATTCTGGAAACCATTATATCGGTGTGCGGTCTTGTAGGCTGTCTGCTGCTGAat                                                                                                   >  7:219324/1‑148 (MQ=255)
                                            tgtgATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATg                                                                                                <  8:269953/148‑1 (MQ=255)
                                             gtgATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATgg                                                                                               >  7:341472/1‑148 (MQ=255)
                                             gtgATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATgg                                                                                               >  1:32709/1‑148 (MQ=255)
                                             gtgATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATgg                                                                                               >  5:102572/1‑148 (MQ=255)
                                                   ttCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGAttt                                                                                         >  8:10033/1‑148 (MQ=255)
                                                          aTGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAaca                                                                                  >  8:250383/1‑148 (MQ=255)
                                                           tgtgAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATGGGCGAGACCATCATATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAacac                                                                                 <  7:346363/148‑1 (MQ=255)
                                                             tgAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACactg                                                                               <  3:16389/148‑1 (MQ=255)
                                                                    tCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCg                                                                        >  3:124815/1‑148 (MQ=255)
                                                                     ccGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGCCCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCg                                                                        >  6:370889/1‑147 (MQ=255)
                                                                       gggCTTCCGGCGGTTCACAGATTACTTTAACCTGACTATCGGCGAGACCCTCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTTGTAGGCTGTCTGCTGCTGAATATTGTGATTTGATAACACTGACTGCCGGAt                                                                     <  4:209893/148‑1 (MQ=255)
                                                                               ggCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGt                                                              >  2:199448/1‑147 (MQ=255)
                                                                                gtttttCAAAGAGTAGTTTTTCCTGACTTTCGGCTTGACGATCAAATCCTTTTAGATGCGGGAAACGATTATCTCGGTGTGCGGTCTGGTAGTCTGTTTGATGCTGAATATGGTGATTTGATAACACTGACTGCAGGATGATGCGTg                                                             <  2:179185/143‑1 (MQ=255)
                                                                                     tCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATTCGGCGTGACCGcc                                                       >  3:156882/1‑148 (MQ=255)
                                                                                     tCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTAGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATa                                                                                     <  4:21442/118‑1 (MQ=255)
                                                                                     tCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTAGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATa                                                                                     >  3:21442/1‑118 (MQ=255)
                                                                                       aaaGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCtt                                                     <  7:10033/148‑1 (MQ=255)
                                                                                        aaGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTa                                                    >  5:160981/1‑148 (MQ=255)
                                                                                        aaGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTa                                                    >  1:193661/1‑148 (MQ=255)
                                                                                         agagTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTAt                                                   >  7:224383/1‑148 (MQ=255)
                                                                                             tACTTTAACCTGACTATCGGCGAGAACATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCgg                                               >  8:4463/1‑148 (MQ=255)
                                                                                                tttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGcc                                             >  6:236098/1‑147 (MQ=255)
                                                                                                tttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGcc                                             >  2:297711/1‑147 (MQ=255)
                                                                                                 ttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGa                                                            <  2:334739/131‑1 (MQ=255)
                                                                                                 ttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGa                                                            >  1:334739/1‑131 (MQ=255)
                                                                                                 ttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTa                                                    <  1:102349/139‑1 (MQ=255)
                                                                                                 ttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTa                                                    >  2:102349/1‑139 (MQ=255)
                                                                                                 ttAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTa                                           >  8:135912/1‑148 (MQ=255)
                                                                                                        gACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCg                                    >  5:374211/1‑148 (MQ=255)
                                                                                                         aCTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCg                                    <  1:297711/147‑1 (MQ=255)
                                                                                                         aCTATCGGCCAGACCATCAAATCCTGGTCGAAGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTTTCCGGCCTACGATTc                                     <  2:51528/146‑1 (MQ=255)
                                                                                                                 cGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATg                                                                                                >  6:328635/1‑79 (MQ=255)
                                                                                                                 cGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATg                                                                                                <  5:328635/79‑1 (MQ=255)
                                                                                                                    gaacATCAATTCCTGGTCGATGCTGGTAACCATTTTCTCGGTGTGGGGTCTGGTAGGCTGTCTGCTGCTGTATATGGTGTTTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTTCGATTCGGGATGAATtagt                        <  8:341472/145‑1 (MQ=38)
                                                                                                                    gaCCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGGCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATtagt                        <  7:379233/148‑1 (MQ=255)
                                                                                                                        aTCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCAGTCTGATGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGTCGTGACCTCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGc                    >  6:167679/1‑148 (MQ=255)
                                                                                                                          cAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCg                  >  7:133361/1‑148 (MQ=255)
                                                                                                                           aaaTCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCg                  >  6:274423/1‑147 (MQ=255)
                                                                                                                            aaTCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGa                <  5:274423/148‑1 (MQ=255)
                                                                                                                                cTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAg            <  4:124815/148‑1 (MQ=255)
                                                                                                                                    tCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGAcgc        >  5:218510/1‑148 (MQ=255)
                                                                                                                                     cGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGAcgcg       >  5:115681/1‑148 (MQ=255)
                                                                                                                                       aTGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGt      >  3:328564/1‑147 (MQ=255)
                                                                                                                                        tGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTc     >  7:376477/1‑147 (MQ=255)
                                                                                                                                           tGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGATGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGTCCGGATAAGACGCGTCAAg  >  4:290516/1‑147 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
GGTGTTAGCCCGGAGCTGATGGTTATTGCGGTTGGTTCCGGTAGTGTGATTTTCTCTCATGTGAACGATCCGGGCTTCTGGCTGTTCAAAGAGTACTTTAACCTGACTATCGGCGAGACCATCAAATCCTGGTCGATGCTGGAAACCATTATCTCGGTGTGCGGTCTGGTAGGCTGTCTGCTGCTGAATATGGTGATTTGATAACACTGACTGCCGGATGCGGCGTGACCGCCTTATCCGGCCTACGATTCGGGATGAATTAGTAGGCCGGATAAGACGCGTCAAG  >  NC_000913/3547675‑3547960

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