Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,875,613 A→G 6.8% K113E (AAA→GAA)  cobS → adenosylcobinamide‑GDP ribazoletransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,875,6130AG 6.8% 167.1 / 2.5 59K113E (AAA→GAA) cobSadenosylcobinamide‑GDP ribazoletransferase
Reads supporting (aligned to +/- strand):  ref base A (28/27);  new base G (2/2);  total (30/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.74e-01

CAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATCGGGCGGGCGGCGATGTTGGGT  >  NC_002947/1875466‑1875759
                                                                                                                                                   |                                                                                                                                                  
cAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGaa                                                                                                                                                   <  3:57832/149‑1 (MQ=255)
 aGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGa                                                                                                                                                    >  6:46115/1‑147 (MQ=255)
   ggTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGa                                                                                                                                                >  6:118246/1‑149 (MQ=255)
            tCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAg                                                                                                                                       <  7:318811/149‑1 (MQ=255)
               ccgcTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTgg                                                                                                                                    >  3:196059/1‑149 (MQ=255)
               ccgcTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTgg                                                                                                                                    >  7:310067/1‑149 (MQ=255)
               ccgcTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTgg                                                                                                                                    >  2:317329/1‑149 (MQ=255)
                 gcTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGc                                                                                                                                  <  8:310067/149‑1 (MQ=255)
                   tGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCg                                                                                                                                <  4:196059/149‑1 (MQ=255)
                      aTGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGAtt                                                                                                                             <  2:122245/149‑1 (MQ=255)
                             tttGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCCCAGTGGGCCGATTGCCGtcg                                                                                                                      >  4:165746/1‑147 (MQ=255)
                                       aCGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGt                                                                                                            <  2:209749/149‑1 (MQ=255)
                                            gtgGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTgg                                                                                                       >  1:47012/1‑149 (MQ=255)
                                            gtgGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTg                                                                                                        <  7:72445/148‑1 (MQ=255)
                                               ggTACCGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGTTCACGCTgg                                                                                                             <  5:238856/140‑1 (MQ=255)
                                               ggTACCGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTgg                                                                                                             >  6:238856/1‑140 (MQ=255)
                                                   cTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGtgttgt                                                                                                 <  4:100015/148‑1 (MQ=255)
                                                      cTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGttgttg                                                                                               <  8:318972/147‑1 (MQ=255)
                                                      cTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCt                                                                                             >  4:125312/1‑149 (MQ=255)
                                                        cAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGGAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGt                                                                                                            <  2:282539/132‑1 (MQ=255)
                                                        cAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGGAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGt                                                                                                            <  8:219251/132‑1 (MQ=255)
                                                        cAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGGAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGt                                                                                                            >  7:219251/1‑132 (MQ=255)
                                                        cAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGGAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGt                                                                                                            >  1:282539/1‑132 (MQ=255)
                                                              tgCCTTGCACCTGGATGGCCTGGCGGACAGTGGAGACGCCTGGCTAGGTGGTTTCGGCTACCGTGAGCGCACCTTGCGTATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGtttt                                                                                     <  3:165746/149‑1 (MQ=255)
                                                                 cTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTgcg                                                                                  >  3:271569/1‑149 (MQ=255)
                                                                  ttGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTgttgtt                                                                                                >  1:105451/1‑134 (MQ=255)
                                                                     cACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTTCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACTCTGGTCTTGGTGTTGTTGCTGAAGTTTTgcg                                                                                  <  2:270202/145‑1 (MQ=255)
                                                                     cACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTgcg                                                                                  >  1:270202/1‑145 (MQ=255)
                                                                      aCCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATtg                                                                             >  6:321919/1‑149 (MQ=255)
                                                                          ggATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTggg                                                                         >  5:117370/1‑149 (MQ=255)
                                                                           gATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGTGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCA‑GCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTg                                                                       <  5:46115/149‑1 (MQ=255)
                                                                              gtgcTGGCGGACAGTGGCGACGCGTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAATGTCCCGCGGAGTGGTCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGc                                                                      <  1:74778/145‑1 (MQ=255)
                                                                                         cAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCgggcagg                                                          <  2:47012/149‑1 (MQ=255)
                                                                                          aGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCgggcaggg                                                         >  2:273727/1‑149 (MQ=255)
                                                                                            tGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCgggcagggca                                                       >  3:71136/1‑149 (MQ=255)
                                                                                              ccGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGttgttg                                                                                               >  3:229786/1‑107 (MQ=255)
                                                                                              ccGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGttgttg                                                                                               <  4:229786/107‑1 (MQ=255)
                                                                                               cGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCg                                                    <  3:170429/149‑1 (MQ=255)
                                                                                                 aCGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCggg                                                  <  6:117370/149‑1 (MQ=255)
                                                                                                   gCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGgc                                                <  5:118246/149‑1 (MQ=255)
                                                                                                        gCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAAt                                           >  1:133198/1‑149 (MQ=255)
                                                                                                           aGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGGGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATtgc                                        <  2:16603/149‑1 (MQ=255)
                                                                                                             gtggtTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATtgt                                                                            >  7:110928/1‑111 (MQ=255)
                                                                                                             gtggtTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATtgt                                                                            <  8:110928/111‑1 (MQ=255)
                                                                                                             gtggtTTCGGCGACCGTCAGCGCACCTTGCTGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTAACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGAATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATtgat                                       >  6:202044/1‑146 (MQ=255)
                                                                                                               ggtTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTgctgct                                    <  3:125312/149‑1 (MQ=255)
                                                                                                                     ggCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCg                                                    <  5:413803/127‑1 (MQ=255)
                                                                                                                     ggCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCg                                                    >  6:413803/1‑127 (MQ=255)
                                                                                                                       cGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGt                            <  1:317329/149‑1 (MQ=255)
                                                                                                                        gACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGtt                           <  1:273727/149‑1 (MQ=255)
                                                                                                                           cGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGtttt                                                                                     <  4:371192/88‑1 (MQ=255)
                                                                                                                           cGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGtttt                                                                                     >  3:371192/1‑88 (MQ=255)
                                                                                                                           cGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGAt                        >  6:134584/1‑149 (MQ=255)
                                                                                                                            gTGATCGAATCTTTCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTTGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTTCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATc                       <  4:271569/149‑1 (MQ=255)
                                                                                                                            gTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGT‑GGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATcg                      >  3:113413/1‑149 (MQ=255)
                                                                                                                              gAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTCGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGt                            <  8:46626/142‑1 (MQ=255)
                                                                                                                              gAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTCGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGt                            >  7:46626/1‑142 (MQ=255)
                                                                                                                                gcgcACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTACTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATcgggc                   >  4:362204/1‑149 (MQ=255)
                                                                                                                                     ccTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGTGCCGTTGATCGggcgggcg              >  8:266742/1‑149 (MQ=255)
                                                                                                                                               cATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATCGGGCGGGCGGCGATGTTgg    >  3:124736/1‑149 (MQ=255)
                                                                                                                                                 tGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATCGGGCGGGCGGCGATGTTGGGt  >  3:105092/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CAGGGTACGCCGTCGCCGCTGCATGCGGCTTTGCTGCTGACGTTGTGGGTACTGCTCAGTGGTGCCTTGCACCTGGATGGCCTGGCGGACAGTGCCGACGCCTGGCTAGGTGGTTTCGGCGACCGTGAGCGCACCTTGCGGATCATGAAAGACCCGCGCAGTGGGCCGATTGCCGTGGTCACGCTGGTCTTGGTGTTGTTGCTGAAGTTTTGCGCATTGTGGGTGCTGGTCGGGCAGGGCATCGGGGCGCAATTGCTGCTGGCGCCGTTGATCGGGCGGGCGGCGATGTTGGGT  >  NC_002947/1875466‑1875759

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