Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,790,625 A→T 7.9% intergenic (‑29/+28) 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,6250AT 7.9% 157.6 / 3.0 60intergenic (‑29/+28)PP_5073/hemEhypothetical protein/uroporphyrinogen decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (30/25);  new base T (2/3);  total (32/28)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGC  >  NC_002947/5790478‑5790771
                                                                                                                                                    |                                                                                                                                                  
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                                                                                                                                               <  1:195067/137‑1 (MQ=255)
                tGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGATCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCgg                                                                                                                                               >  2:195067/1‑137 (MQ=255)
                  aaCAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTgtgcagt                                                                                                                                 <  6:256742/149‑1 (MQ=255)
                       gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                             >  8:101849/1‑148 (MQ=255)
                       gggggCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCAt                                                                                                                             >  4:102193/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                                                                           <  7:101849/147‑1 (MQ=255)
                                                                          tttcggtttcgATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAAt                                                                           <  3:273453/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                                 >  7:148180/1‑149 (MQ=255)
                                                                                                                  ttCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTgg                                 >  3:266841/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)
                                                                                                                                                   aaGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTgg   >  6:295191/1‑147 (MQ=255)
                                                                                                                                                   aaGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTgg   <  5:295191/147‑1 (MQ=255)
                                                                                                                                                   aaGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGGGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGc  >  5:258714/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
GCAAATGCGCCAGTGCTGAACAAGGGGGCGAGGATGAGGGCGAGATAACGGATTTTCATGATAATCGGGTCCTGTTTCGGTTTCGATGGGACCAGATTAGCGACGATCGCTGAATTCAAACTTAACGTACAAAAAACCCGGCAT‑CGAAGCCGGGTTTACTGTGCAGTGCATCAAGTCACTGGTGGTACTGCGCCGACAGTTCGTGCACGGCGTTGATGAATGCACCCGCATGTTCCGGGTCAACTTCCGGCGTGATGCCGTGGCCAAGGTTGAACACATGGCCAGTGCCGTGGC  >  NC_002947/5790478‑5790771

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