Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_002947 1,647,927 G→A Y265Y (TAC→TAT)  oprB‑II ← carbohydrate‑selective porin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,647,9270GA98.4% 207.9 / ‑3.7 63coding (795/1335 nt)oprB‑IIcarbohydrate‑selective porin
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (30/32);  minor base . (1/0);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTC  >  NC_002947/1647783‑1648067
                                                                                                                                                |                                                                                                                                            
gTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaata                                                                                                                                            <  5:115592/147‑1 (MQ=255)
 tCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTgtaatagta                                                                                                                                         <  3:98395/149‑1 (MQ=255)
          aGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCAc                                                                                                                                 >  4:286714/1‑148 (MQ=255)
                gTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTc                                                                                                                          >  4:110990/1‑149 (MQ=255)
                   aCCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTcccc                                                                                                                       >  7:239109/1‑149 (MQ=255)
                    ccTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCt                                                                                                                      <  8:255688/149‑1 (MQ=255)
                               gcgcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTg                                                                                                           <  5:297714/149‑1 (MQ=255)
                                 gcCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGa                                                                                                         <  5:350195/149‑1 (MQ=255)
                                    aaTCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCACGCCTTGGATAc                                                                                                      >  8:330157/1‑149 (MQ=255)
                                        cacaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTg                                                                                                   <  4:357046/148‑1 (MQ=255)
                                         acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGat                                                                                                                              >  1:134410/1‑120 (MQ=255)
                                         acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGat                                                                                                                              <  2:134410/120‑1 (MQ=255)
                                         acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                 >  3:62996/1‑149 (MQ=255)
                                         acaAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACAGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTggg                                                                                                 >  5:288701/1‑149 (MQ=255)
                                          caAGCCGTGCTTTCTCGAACTGCTGCGGGAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTgggg                                                                                                <  7:311516/149‑1 (MQ=255)
                                           aaGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCGCCTTTCAAGCCTTGGATACGTGGGGt                                                                                               <  2:119695/149‑1 (MQ=255)
                                                 tgcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGCTACGTGGGGTCCAtcc                                                                                         >  7:334798/1‑147 (MQ=255)
                                                 tgcttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAc                                                                                         >  2:211202/1‑149 (MQ=255)
                                                   cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTg                                                                                                           <  8:156693/129‑1 (MQ=255)
                                                   cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTg                                                                                                           >  7:156693/1‑129 (MQ=255)
                                                   cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCa                                                                                       <  8:45982/149‑1 (MQ=255)
                                                   cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCa                                                                                       <  4:192132/149‑1 (MQ=255)
                                                   cttgctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCa                                                                                       <  3:286714/149‑1 (MQ=255)
                                                   cttgctCGAAATGCTGCGGTAGGCGGCGCCGCTGAGGTCGTCCGTCTGAGCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATAcaa                                                                                       <  2:278120/149‑3 (MQ=255)
                                                   cttgatCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGAGGCCGGCTGACCGTTGATGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATTGTAGCCGGAACGATATTCCCATTTCAAGACTTGGATAAGTGGGGTCCTTag                                                                                         >  8:147692/1‑146 (MQ=255)
                                                       ctCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTc                                                                                   >  2:89928/1‑149 (MQ=255)
                                                        tCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTc                                                                                   >  1:185886/1‑148 (MQ=255)
                                                            aCTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACg                                                                              <  7:19327/149‑1 (MQ=255)
                                                             ctgctgCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACgg                                                                             <  6:160259/149‑1 (MQ=255)
                                                                ctgcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTA‑TAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCAt                                                                          >  2:331516/1‑148 (MQ=255)
                                                                  gcGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTa                                                                        >  3:134484/1‑149 (MQ=255)
                                                                         gcggcgCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCAcc                                                                 <  4:24849/149‑1 (MQ=255)
                                                                             cgccgcTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTgg                                                             <  4:62996/149‑1 (MQ=255)
                                                                                 gcTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgc                                                         <  3:110990/149‑1 (MQ=255)
                                                                                  cTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGtgcc                                                        <  7:99520/149‑1 (MQ=255)
                                                                                      gggCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCCGGGtgccgctg                                                    >  5:183223/1‑149 (MQ=255)
                                                                                               gcTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGc                                           >  5:41882/1‑149 (MQ=255)
                                                                                                 tGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCtt                                         <  6:288701/149‑1 (MQ=255)
                                                                                                   aCCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCtt                                         <  4:134484/147‑1 (MQ=255)
                                                                                                        tGCTGTACTTGAGCACATCTTCTGGCATGGCATTACTGTAATAGTAGCCGTCACGATATTCCCCTTTCAAGCCTTGGATACCTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGATCTTGAAACCa                                  <  7:147692/149‑1 (MQ=255)
                                                                                                           tGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCGGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTg                               >  1:259540/1‑149 (MQ=255)
                                                                                                             tCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAATCCTTGGATACGTGGGGACAATACCACTTCGAACGGCATTACAGCACCCTGGGTGCCGCTGCCGATGAGATTGAAACCAtttac                             >  3:315739/1‑146 (MQ=255)
                                                                                                               cTTGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCg                            >  2:104139/1‑148 (MQ=255)
                                                                                                                ttGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                          >  5:310020/1‑149 (MQ=255)
                                                                                                                ttGAGAACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGa                          >  5:73442/1‑149 (MQ=255)
                                                                                                                  gagaACATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGAtt                        >  5:362570/1‑149 (MQ=255)
                                                                                                                     aaCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTc                       >  1:371458/1‑147 (MQ=255)
                                                                                                                      aCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc                             >  6:142358/1‑140 (MQ=255)
                                                                                                                      aCATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGcc                             <  5:142358/140‑1 (MQ=255)
                                                                                                                       cATCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAgg                   >  7:329472/1‑149 (MQ=255)
                                                                                                                        aTCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCACACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGt                  <  8:239109/149‑1 (MQ=255)
                                                                                                                         tCTTGTGCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGtt                 >  3:209519/1‑149 (MQ=255)
                                                                                                                             gtgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGtt                 <  6:41882/145‑1 (MQ=255)
                                                                                                                              tgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGaa             <  4:209519/148‑1 (MQ=255)
                                                                                                                              tgCCTTGGCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAg            <  6:362570/149‑1 (MQ=255)
                                                                                                                                 catGGAATTACTGTAAAAGTAGACGGCACGAAATTCCCATTTAAAGCCTTGGATACCTGGGGTCCATACAAGTTCGAACGGAATTACCGCACCCTCGGTGCCCCTGCCGCTGAGCTTGAAACCATAGCCGGATTCGAGGAAGGAAGGGt         >  4:168468/3‑149 (MQ=255)
                                                                                                                                    ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt        <  2:185886/147‑1 (MQ=255)
                                                                                                                                    ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt        <  4:271602/147‑1 (MQ=255)
                                                                                                                                    ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGtt        <  1:211202/147‑1 (MQ=255)
                                                                                                                                    ggCATTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTc       <  1:104139/148‑1 (MQ=255)
                                                                                                                                       aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCt   >  2:337588/1‑149 (MQ=255)
                                                                                                                                       aTTACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCt   <  2:259540/149‑1 (MQ=255)
                                                                                                                                        ttACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc  >  3:271263/1‑149 (MQ=255)
                                                                                                                                        ttACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc  <  1:337588/149‑1 (MQ=255)
                                                                                                                                         tACTGTAATAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTc  >  7:169065/1‑148 (MQ=255)
                                                                                                                                                |                                                                                                                                            
GTCGGACGCCAGCGAGGTGACCTGCTGCTGGGCGCCAATCCACAAGCCGTGCTTGCTCGAACTGCTGCGGTAGGCGGCGCCGCTGAGGGCGGCCGGCTGACCGTTGCTGTCCTTGAGAACATCTTGTGCCTTGGCATTACTGTAGTAGTAGCCGGCACGATATTCCCCTTTCAAGCCTTGGATACGTGGGGTCCATACCAGTTCGAACGGCATTACCGCACCCTGGGTGCCGCTGCCGCTGAGCTTGAAACCATTGCCGGATTCGAGGTTGGAAGGGTTCTGCTC  >  NC_002947/1647783‑1648067

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: ‑