Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,620 A→T 8.0% intergenic (‑24/+33) PP_5073 ← / ← hemE hypothetical protein/uroporphyrinogen decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,790,6200AT 8.0% 175.0 / 4.1 62intergenic (‑24/+33)PP_5073/hemEhypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (32/25);  new base T (2/3);  total (34/28)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

GCCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCC  >  NC_002947/5790475‑5790766
                                                                                                                                                 |                                                                                                                                                   
gccGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑c                                                                                                                                                  >  2:31384/1‑148 (MQ=255)
  cGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑Cg                                                                                                                                                 >  3:275453/1‑147 (MQ=255)
   gCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAg                                                                                                                                              >  5:10354/1‑149 (MQ=255)
     aaaTGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGATATAACGGATTTTCATGATAATAGGGTCCTGTTTCGGTTTCGATGGGACCACATTAGAGACGATCGCTGAATTCAAACATAACTTACTAAAAACCCGGCAT‑CGAATcc                                                                                                                                            >  7:104649/1‑149 (MQ=255)
          cgcCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTT‑CGATGCCGGGtt                                                                                                                                       <  7:36450/149‑1 (MQ=255)
             cAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCTT‑CGATGCCGGGttttt                                                                                                                                    <  6:291248/149‑3 (MQ=255)
              aGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACt                                                                                                                                   <  1:266920/149‑1 (MQ=255)
                  cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACtgtg                                                                                                                                >  6:123880/1‑148 (MQ=255)
                  cTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACtgtg                                                                                                                                <  3:279145/148‑1 (MQ=255)
                   tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGATCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCgg                                                                                                                                          >  2:195067/1‑137 (MQ=255)
                   tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGATCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCgg                                                                                                                                          <  1:195067/137‑1 (MQ=255)
                     aaCAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTgtgcagt                                                                                                                            <  6:256742/149‑1 (MQ=255)
                          gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                        >  4:102193/1‑148 (MQ=255)
                          gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                        >  8:101849/1‑148 (MQ=255)
                          gggggCCAGGATGAGGGCGAGATAACGGTTTTTAATGTTTATCGGTTCCTGTTTAGGTTAAGATGGGACCAGATTAGCGACGTACGATGAATTCAACGTTAACTTACAAAAAACCC‑GCATGGGAAGCCGGTTTTAATgtgcagtgatt                                                                                                                        >  2:162459/1‑145 (MQ=255)
                                gAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCa                                                                                                                 >  6:239547/1‑149 (MQ=255)
                                   gATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtg                                                                                                              <  1:31384/149‑1 (MQ=255)
                                       aGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                           <  6:10354/148‑1 (MQ=255)
                                            gagaTAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTg                                                                                                     >  6:65009/1‑149 (MQ=255)
                                                taccGGATTTTCATGATAATCGGGTCCTGTTGCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGcc                                                                                                 <  3:266600/146‑1 (MQ=255)
                                                tAACGGATTTTCATTATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGcc                                                                                                 <  3:262363/149‑1 (MQ=255)
                                                 aaCGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGcc                                                                                                 >  1:13682/1‑148 (MQ=255)
                                                    ggATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGa                                                                                               >  7:22628/1‑147 (MQ=255)
                                                        tttCATGATAATCGGGTTCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTc                                                                                         >  6:14018/1‑149 (MQ=255)
                                                           cATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGt                                                                                       <  2:13682/148‑1 (MQ=255)
                                                                   tCGCGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGt                                                                              >  2:148463/1‑149 (MQ=255)
                                                                      ggTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTAACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtga                                                                           >  5:229002/1‑149 (MQ=255)
                                                                       gTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgat                                                                          <  5:123880/149‑1 (MQ=255)
                                                                        tCCTGTTTCGGTTTCGATGGGACCAGATTAGAGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTtgat                                                                          >  5:146172/1‑148 (MQ=255)
                                                                          cTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGaa                                                                       <  7:114615/149‑1 (MQ=255)
                                                                            gtttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATg                                                                     >  2:170033/1‑149 (MQ=255)
                                                                             tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                      <  3:273453/147‑1 (MQ=255)
                                                                             tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                      <  7:101849/147‑1 (MQ=255)
                                                                                   tttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCg                                                               >  6:311695/1‑148 (MQ=255)
                                                                                       gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                           >  7:112508/1‑100 (MQ=255)
                                                                                       gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACtggtg                                                                                                           <  8:112508/100‑1 (MQ=255)
                                                                                       gATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGAGCATCAAGTCACTGGTGGTACTGCGCCGACAGATCGTGCACGGCGTTTCTGAATGCACCCGCATGt                                                          >  4:242984/1‑149 (MQ=255)
                                                                                             aCCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCggg                                                    <  7:25437/149‑1 (MQ=255)
                                                                                               cAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTc                                                  <  2:16521/149‑1 (MQ=255)
                                                                                                 gATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCaa                                                <  6:229002/149‑1 (MQ=255)
                                                                                                  aTTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAAc                                               >  8:51221/1‑149 (MQ=255)
                                                                                                      gCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTcc                                           <  5:65009/149‑1 (MQ=255)
                                                                                                        gacgaTCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCgg                                         >  6:20606/1‑149 (MQ=255)
                                                                                                                 tGAATTCAAACTTAACGTACAAAAAACGCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                >  8:175552/1‑149 (MQ=255)
                                                                                                                 tGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGcc                                >  7:114548/1‑149 (MQ=255)
                                                                                                                    aTTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                             <  5:239547/149‑1 (MQ=255)
                                                                                                                    aTTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                             >  6:303119/1‑149 (MQ=255)
                                                                                                                    aTTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTg                             <  8:22628/149‑1 (MQ=255)
                                                                                                                     ttCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                            >  3:266841/1‑149 (MQ=255)
                                                                                                                     ttCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                            >  7:148180/1‑149 (MQ=255)
                                                                                                                       cAAACTTAACGTACAAAAAACCCGGAAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGc                           <  3:102193/148‑1 (MQ=255)
                                                                                                                        aaaCTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCa                         <  5:311695/149‑1 (MQ=255)
                                                                                                                               aCGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTg                   >  5:33850/1‑148 (MQ=255)
                                                                                                                                cGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCTACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGaa                 <  6:146172/149‑1 (MQ=255)
                                                                                                                                 gTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGaa                 >  5:277727/1‑148 (MQ=255)
                                                                                                                                        aaaCCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc         <  8:114548/149‑1 (MQ=255)
                                                                                                                                         aaCCCGGCTT‑CGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc         >  5:182101/1‑148 (MQ=255)
                                                                                                                                         aaCCCGGCTT‑CGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGc         >  8:275955/1‑148 (MQ=255)
                                                                                                                                          aCCCGGCTT‑CGATGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCa       <  6:182101/149‑1 (MQ=255)
                                                                                                                                            ccGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGGATGTTCCGGGTCATCTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGt     <  1:148463/149‑1 (MQ=255)
                                                                                                                                             cGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTg    >  3:236295/1‑149 (MQ=255)
                                                                                                                                               gCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGcc  >  2:19814/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
GCCGCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCC  >  NC_002947/5790475‑5790766

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