Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,246,242	A→G	18.3%	intergenic (+31/+60)	ibpA → / ← PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,246,242	0	A	G	18.3%	123.8 / 21.7	60	intergenic (+31/+60)	ibpA/PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (24/25);  new base G (6/5);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.35e-01

AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAT  >  NC_002947/2246098‑2246387                                                                                                                                                 |                                                                                                                                                    aaGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                  <  5:323363/148‑1 (MQ=255)       gccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑t                                                                                                                                         <  4:188807/149‑1 (MQ=255)          ggCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGGTGAACTGATCGAGCTGCAACAAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttt                                                                                                                                      <  4:84139/149‑1 (MQ=255)           gCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACG‑AAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                    <  1:302804/149‑1 (MQ=255)             cTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCg                                                                                                                                   >  1:227665/1‑149 (MQ=255)                gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttc                                                                                                                                  <  7:376402/147‑4 (MQ=255)                gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGc                                                                                                                                <  4:356047/149‑1 (MQ=255)                gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGc                                                                                                                                <  1:207486/149‑1 (MQ=255)                gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGc                                                                                                                                >  2:376743/1‑149 (MQ=255)                gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGc                                                                                                                                >  2:207486/1‑149 (MQ=255)                 ccAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCCATGGCGAAAAACCAGCGCAGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCACAGGTGCCC‑‑TTTTTCGTGCt                                                                                                                               >  1:290038/1‑149 (MQ=255)                  cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCt                                                                                                                               >  8:277876/1‑148 (MQ=255)                           gCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAg                                                                                                                     >  4:153079/1‑149 (MQ=255)                           gCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAg                                                                                                                     >  8:59507/1‑149 (MQ=255)                                 cATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCcgc                                                                                                               >  6:74048/1‑149 (MQ=255)                                   tCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                             <  7:143321/149‑1 (MQ=255)                                     gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGc                                                                                                                                              <  6:200834/116‑1 (MQ=255)                                     gACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGc                                                                                                                                              >  5:200834/1‑116 (MQ=255)                                      aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAgg                                                                                                          >  7:214492/1‑149 (MQ=255)                                          gctgcGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTgcg                                                                                                      >  2:391241/1‑149 (MQ=255)                                               gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGcc                                                                                                 <  1:301958/149‑1 (MQ=255)                                                     tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATg                                                                                           <  6:171192/149‑1 (MQ=255)                                                       ccTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGc                                                                                          >  6:394753/1‑148 (MQ=255)                                                        cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATTGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGa                                                                                        <  7:59507/149‑1 (MQ=255)                                                           aagaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGAt                                                                                       <  3:358566/147‑1 (MQ=255)                                                            agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAg                                                                                    <  1:376743/149‑1 (MQ=255)                                                             gaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAg                                                                                    >  5:194388/1‑148 (MQ=255)                                                                gccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAAt                                                                                <  1:391241/149‑1 (MQ=255)                                                                  caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCt                                                                              >  5:287013/1‑149 (MQ=255)                                                                   aaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTg                                                                             >  6:182424/1‑149 (MQ=255)                                                                         aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                       >  7:215700/1‑149 (MQ=255)                                                                           gcgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGcc                                                                     >  1:261671/1‑149 (MQ=255)                                                                            cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                             >  3:134590/1‑108 (MQ=255)                                                                            cgcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGt                                                                                                             <  4:134590/108‑1 (MQ=255)                                                                             gcATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTc                                                                   <  6:206512/149‑1 (MQ=255)                                                                                   cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                         <  3:275098/105‑1 (MQ=255)                                                                                   cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGt                                                                                                         >  4:275098/1‑105 (MQ=255)                                                                                   cGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                             <  5:394753/149‑1 (MQ=255)                                                                                          tGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGc                                                      <  8:214492/149‑1 (MQ=255)                                                                                              gAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGa                                                   >  6:106245/1‑148 (MQ=255)                                                                                               aaaaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                                 <  7:277876/149‑1 (MQ=255)                                                                                                      aGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaa                                                  >  2:50306/1‑141 (MQ=255)                                                                                                      aGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaa                                                  <  1:50306/141‑1 (MQ=255)                                                                                                       gcgcTGAACTGATCGAGCTGCAACGAAAAACAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggta                                         <  2:290038/149‑1 (MQ=255)                                                                                                              aCTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacggtaa                                  <  2:227665/149‑1 (MQ=255)                                                                                                               cTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAAt                                 <  8:215700/149‑1 (MQ=255)                                                                                                                    cGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAAc                                                            <  2:64929/117‑1 (MQ=255)                                                                                                                    cGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAAc                                                            >  1:64929/1‑117 (MQ=255)                                                                                                                       gCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCt                         <  5:106245/149‑1 (MQ=255)                                                                                                                             gcacgaaaaaGGGCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                      >  3:107087/6‑148 (MQ=38)                                                                                                                             aCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                      >  8:91642/1‑146 (MQ=255)                                                                                                                              cGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                      >  4:30461/1‑129 (MQ=255)                                                                                                                              cGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacg                                      <  3:30461/129‑1 (MQ=255)                                                                                                                                   aaaaGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAcc             >  2:100085/1‑149 (MQ=255)                                                                                                                                      aGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg          >  5:370109/1‑149 (MQ=255)                                                                                                                                      aGGGCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTg          >  1:297600/1‑149 (MQ=255)                                                                                                                                       gggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGc                     >  6:325375/1‑137 (MQ=255)                                                                                                                                       gggCACCTCTGGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGc                     <  5:325375/137‑1 (MQ=255)                                                                                                                                             cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGa   <  7:91642/149‑1 (MQ=255)                                                                                                                                             cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCAGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGa   >  6:167461/1‑149 (MQ=255)                                                                                                                                                cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAt  <  3:153079/147‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                    AAGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGGTGCGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGAT  >  NC_002947/2246098‑2246387

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