Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,199,963:1	+C	6.5%	coding (829/1278 nt)	PP_1944 →	aminomethyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,199,963	1	.	C	6.5%	245.2 / 10.5	62	coding (829/1278 nt)	PP_1944	aminomethyltransferase
Reads supporting (aligned to +/- strand):  ref base . (31/27);  new base C (2/2);  total (33/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 = 1.00e+00

GTTACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199816‑2200111                                                                                                                                                      |                                                                                                                                                     gTTACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑t                                                                                                                                                     >  5:94872/1‑149 (MQ=255)        aCTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctg                                                                                                                                              >  1:286926/1‑149 (MQ=255)         cTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                             <  8:242989/149‑1 (MQ=255)            cATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAGATCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgg                                                                                                                                             <  1:260540/147‑1 (MQ=255)             aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgggctg                                                                                                                                         >  6:376483/1‑149 (MQ=255)             aTGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAgctgggctg                                                                                                                                         >  8:246029/1‑149 (MQ=255)              tGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAGCTCGa‑ctgggctgc                                                                                                                                        <  4:151904/149‑2 (MQ=255)                 ccGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCtgggctgggc                                                                                                                                      >  7:308328/1‑148 (MQ=255)                     gATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCaaaa                                                                                                                                 >  8:270338/1‑149 (MQ=255)                       ttttCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCTAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAAcc                                                                                                                               <  8:308328/149‑1 (MQ=255)                           cGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATc                                                                                                                            >  3:91044/1‑148 (MQ=255)                                cAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTc                                                                                                                      >  8:21146/1‑149 (MQ=255)                                 agagATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCg                                                                                                                     <  5:279786/149‑1 (MQ=255)                                  gagaTGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTTAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGa                                                                                                                    <  2:80517/149‑1 (MQ=255)                                      tGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGTCCAAAACCATCAAGTACGATCAc                                                                                                                >  2:212850/1‑149 (MQ=255)                                         cGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCTATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGAc                                                                                                             <  4:280349/149‑1 (MQ=255)                                          gggTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACt                                                                                                            >  2:352751/1‑149 (MQ=255)                                             ttATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTa                                                                                                         >  2:168678/1‑149 (MQ=255)                                               aTCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATc                                                                                                       >  3:362394/1‑149 (MQ=255)                                                   cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACTACTTTATCGGTc                                                                                                   >  1:370771/1‑149 (MQ=255)                                                   cAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTc                                                                                                   >  1:366519/1‑149 (MQ=255)                                                    aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                              >  4:197189/1‑137 (MQ=255)                                                    aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGa                                                                                                              <  3:197189/137‑1 (MQ=255)                                                    aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCgg                                                                                                     <  5:376483/146‑1 (MQ=255)                                                    aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCgg                                                                                                     <  1:168678/146‑1 (MQ=255)                                                    aaGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGGTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGATGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCg                                                                                                  >  1:83064/1‑149 (MQ=255)                                                             cgggcgGGGATGGCAATGTTTGCCTCGACTTTCAATGTTGGGGGGAGCTTTGAG‑TCGAATTATGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCGAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         <  5:197297/146‑1 (MQ=255)                                                             cgcgcgGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         <  6:94872/149‑1 (MQ=255)                                                                  gggTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGaaa                                                                                    <  7:246029/149‑1 (MQ=255)                                                                                    cTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCt                                                                   <  1:352751/148‑1 (MQ=255)                                                                                    cTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGACTc                                                                  >  1:248735/1‑149 (MQ=255)                                                                                     tCGACTTTCAATGTTGCGGGGAGCTTTTAG‑TAGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAAACATCAAGTTCGATAACGACTTTATCGGTCGTGCCGCGTTGGAAAAACAGGTTGCAGAGACTc                                                                  >  5:254747/1‑148 (MQ=255)                                                                                     tCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCTATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGcc                                                                    >  5:375441/1‑146 (MQ=255)                                                                                       gtcTTTCACTGTTGGGGTGATCTTTGTG‑TGGATTAAGGTATGGGTTTGGTTTCTCAGCCCTG‑TCGTTGTGGGCTGGGCCATAAGCATCAAGTTCGATCACGACTTTATCGGTCGTGCGGCGTTGGAAAAAGAGGTTGCAGAGCCTCg                                                                 <  6:254747/145‑1 (MQ=255)                                                                                           ttCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGt                                                           >  1:222649/1‑149 (MQ=255)                                                                                              aaTGTTGCGGGGCGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGAt                                                        <  8:287491/149‑1 (MQ=255)                                                                                               aTGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATg                                                       >  4:338005/1‑149 (MQ=255)                                                                                                tGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATgc                                                      >  6:109968/1‑149 (MQ=255)                                                                                                  ttGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAgg                                                                              >  8:133979/1‑123 (MQ=255)                                                                                                  ttGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAgg                                                                              <  7:133979/123‑1 (MQ=255)                                                                                                    gCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         >  8:231252/1‑110 (MQ=255)                                                                                                    gCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGtt                                                                                         <  7:231252/110‑1 (MQ=255)                                                                                                      ggggAGCTTTGAT‑TCGAATAAGGTAAGCGATTGGTATGGCACCCCAG‑TCGAGCGGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCt                                                <  2:366519/149‑1 (MQ=255)                                                                                                        ggAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                              >  3:238850/1‑149 (MQ=255)                                                                                                        ggAGCTTTGAG‑TATAATAATGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTTGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTg                                              <  6:375441/149‑1 (MQ=255)                                                                                                         taGCTTTGAG‑TCTAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCTATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTgg                                             <  1:212850/148‑1 (MQ=255)                                                                                                            cTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGttt                                          <  5:109968/149‑1 (MQ=255)                                                                                                              ttGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTAAAGGTCGTGCCGCGTTGGATAAAGAGGTTGCAGAGCCTCGACGGGTGATGGGCACCGTGGTTTg                                         >  2:164764/1‑148 (MQ=255)                                                                                                               tGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTgg                                        <  4:238850/148‑1 (MQ=255)                                                                                                               tGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTg                                         <  7:270338/147‑1 (MQ=255)                                                                                                                 aG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCa                                                                         <  3:167784/113‑1 (MQ=255)                                                                                                                 aG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCa                                                                         >  4:167784/1‑113 (MQ=255)                                                                                                                       aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtga                                >  5:293357/1‑149 (MQ=255)                                                                                                                       aaTAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtga                                <  7:21146/149‑1 (MQ=255)                                                                                                                          aaGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGtgatga                             >  3:222587/1‑149 (MQ=255)                                                                                                                                     ttGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATg                  <  2:286926/149‑1 (MQ=255)                                                                                                                                        gTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGtt                <  4:222587/148‑1 (MQ=255)                                                                                                                                            cccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGc           >  3:125649/3‑149 (MQ=255)                                                                                                                                             ccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGt                                            >  8:184060/2‑115 (MQ=255)                                                                                                                                             ccAGCCCAGCTCGA‑CTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGt                                            <  7:184060/114‑1 (MQ=255)                                                                                                                                                gCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCg       >  5:178335/1‑149 (MQ=255)                                                                                                                                                 cccAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGc      >  5:49505/1‑149 (MQ=255)                                                                                                                                                  ccAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCt     >  5:209688/1‑149 (MQ=255)                                                                                                                                                     g‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGtt  >  2:56192/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                     GTTACTGACTACATGCCCGCGATTTTCGACCCAGAGATGGCGGGTTATCTCAAGGAATTCCGCGCGGGTATGCCAAGTTTTGCCTCGACTTTCAATGTTGCGGGGAGCTTTGAG‑TCGAATAAGGTAAGCGATTGGTATCGCAGCCCAG‑TCGAGCTGGGCTGGGCCAAAACCATCAAGTTCGATCACGACTTTATCGGTCGTGCCGCGTTGGAAAAAGAGGTTGCAGAGCCTCGACGGGTGATGCGCACCTTGGTTTGGAATGGTGATGACGTTGCAGATGTTCATGCCTCGCTGTT  >  NC_002947/2199816‑2200111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑