Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,576,047 T→C 11.8% H148H (CAT→CAC ybhP‑II → phosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,576,0470TC11.8% 146.8 / 10.7 59H148H (CAT→CACybhP‑IIphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base T (27/25);  new base C (2/5);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.92e-01

CGGTCTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCC  >  NC_002947/4575899‑4576195
                                                                                                                                                     |                                                                                                                                                    
cgGTCTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCAt                                                                                                                                                      <  5:52683/149‑1 (MQ=255)
        ccTGACGGCCACCACGGCAATGCGTTAGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCt                                                                                                                                               >  1:168976/1‑149 (MQ=255)
         cTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGc                                                                                                                                             <  1:152851/149‑1 (MQ=255)
                 caccacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                     >  2:214132/1‑149 (MQ=255)
                  accacGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCACCTg                                                                                                                                    <  6:28371/149‑1 (MQ=255)
                         cAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGc                                                                                                                             <  1:214132/149‑1 (MQ=255)
                         cAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGc                                                                                                                             <  4:123209/149‑1 (MQ=255)
                         cAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAACGACAGGTGCACGCAGACCGCATGCACCTGGCCTTGc                                                                                                                             <  3:92687/149‑1 (MQ=38)
                            tGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGACAGGTGCACGCAGACCGCATGCACCTGGCCTTGcccc                                                                                                                          <  4:182819/149‑4 (MQ=38)
                                     tGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                  >  4:303116/1‑149 (MQ=255)
                                       tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGTCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCACCt                                                                                                                                     <  4:210002/128‑1 (MQ=255)
                                       tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGTCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCACGCAGACCGCATGCACCt                                                                                                                                     >  3:210002/1‑128 (MQ=255)
                                       tCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCAc                                                                                                                  <  3:303116/147‑1 (MQ=255)
                                             tACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCaaaa                                                                                                          <  5:34153/149‑1 (MQ=255)
                                                ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                     <  8:291032/119‑1 (MQ=255)
                                                ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCt                                                                                                                                     >  7:291032/1‑119 (MQ=255)
                                                ccGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAag                                                                                                         >  3:200809/1‑147 (MQ=255)
                                                 cGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGc                                                                                                                                         >  2:175717/1‑114 (MQ=255)
                                                 cGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGc                                                                                                                                         <  1:175717/114‑1 (MQ=255)
                                                       tcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCCAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCAACCGGCAAAAGCAGCTGCAg                                                                                                >  6:21812/1‑149 (MQ=255)
                                                       tcGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAg                                                                                                >  6:26371/1‑149 (MQ=255)
                                                           gCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgc                                                                                            >  8:104553/1‑149 (MQ=255)
                                                            cACTGCAACCTGGACGTGTCGATCACCGGCCCCGAACTCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCTGTCTGCGTGCACCTGTCGTTGCTGTAAAGCCACCGGCAAAAGCAGCTGCAGTtgct                                                                                           <  7:178625/149‑1 (MQ=255)
                                                              ccGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTtgctgc                                                                                         <  8:248952/149‑1 (MQ=255)
                                                               cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCg                                                                                        >  3:319338/1‑149 (MQ=255)
                                                               cGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGAGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCg                                                                                        >  3:245590/1‑149 (MQ=255)
                                                                         aCGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActgctg                                                                              <  4:200809/149‑1 (MQ=255)
                                                                         aCGTGTCGAGCACCGGCCCCGAACGACGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActgctg                                                                              <  4:118180/149‑1 (MQ=255)
                                                                                      cGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActg                                                                                 <  8:149463/133‑1 (MQ=255)
                                                                                      cGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAActg                                                                                 >  7:149463/1‑133 (MQ=255)
                                                                                      cGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCgg                                                                 <  4:266164/149‑1 (MQ=255)
                                                                                           ccGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAAc                                                            >  1:130540/1‑149 (MQ=255)
                                                                                             gAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGc                                                          >  1:232431/1‑149 (MQ=255)
                                                                                                         cTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcat                                              <  6:188449/149‑1 (MQ=255)
                                                                                                          tGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGatcatc                                             >  6:54930/1‑149 (MQ=255)
                                                                                                           gTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCg                                            >  1:106952/1‑149 (MQ=255)
                                                                                                           gTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCg                                            >  4:119451/1‑149 (MQ=255)
                                                                                                              gCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCg                                         >  1:65830/1‑149 (MQ=255)
                                                                                                               cACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCgg                                        >  5:154935/1‑149 (MQ=255)
                                                                                                                aCTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGc                                       <  7:152939/149‑1 (MQ=255)
                                                                                                                      tgCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACtt                                  >  3:299438/1‑148 (MQ=255)
                                                                                                                        cTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATTATCGCCGGCGACTTCaa                               <  2:305733/149‑1 (MQ=255)
                                                                                                                         tGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAc                              <  2:168976/149‑1 (MQ=255)
                                                                                                                          ggAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAAc                              >  4:32822/1‑148 (MQ=255)
                                                                                                                              gTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTg                         >  8:22327/1‑149 (MQ=255)
                                                                                                                                gCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGa                       <  3:32822/149‑1 (MQ=255)
                                                                                                                                 ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGat                                                 <  6:169492/122‑1 (MQ=255)
                                                                                                                                 ccGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCAGGCCAACGCCCCGGTGat                                                 >  5:169492/1‑122 (MQ=255)
                                                                                                                                  cGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                      >  4:84759/1‑148 (MQ=255)
                                                                                                                                   ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  3:119451/149‑1 (MQ=255)
                                                                                                                                   ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaaaa                    <  2:232431/149‑1 (MQ=255)
                                                                                                                                   ggggCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGaa                      >  4:117775/1‑147 (MQ=255)
                                                                                                                                      gCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGcc                 <  4:245590/149‑1 (MQ=255)
                                                                                                                                        aacgacAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCa                >  3:273357/6‑148 (MQ=255)
                                                                                                                                        aacgacAGGTGCACGCAGACCGCATGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCa                <  4:273357/143‑1 (MQ=255)
                                                                                                                                        aaGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCAt               >  4:64080/1‑149 (MQ=255)
                                                                                                                                             cAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCaa          <  7:104553/149‑1 (MQ=255)
                                                                                                                                                   cATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAc    <  3:117775/149‑1 (MQ=255)
                                                                                                                                                    aTGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTg                                                                                                   >  2:185339/1‑53 (MQ=255)
                                                                                                                                                    aTGCGGTCTGCGTGCACCTGTCGTAGCTGGAAAGCCACCGGCAAAAGCAGCTg                                                                                                   <  1:185339/53‑1 (MQ=255)
                                                                                                                                                     tGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCAccc  <  5:54930/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CGGTCTACCCTGACGGCCACCACGGCAATGCGTTGCTGTCCAAATACCCGATCATCGAGCACCGCAACCTGGACGTGTCGATCACCGGCCCCGAACGCCGCGGCCTGTTGCACTGTGTGCTGGAGGTGCCGGGGCAAGGCCAGGTGCATGCGGTCTGCGTGCACCTGTCGTTGCTGGAAAGCCACCGGCAAAAGCAGCTGCAGTTGCTGCGTAAACTGCTGGACGCCCTGCCGGCCAACGCCCCGGTGATCATCGCCGGCGACTTCAACGACTGGAAAAGCCATGGCAACCGCACCC  >  NC_002947/4575899‑4576195

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: