Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,637,325	A→T	19.1%	intergenic (+34/+32)	pdxJ → / ← PP_1437	pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,637,325	0	A	T	19.1%	123.4 / 18.3	57	intergenic (+34/+32)	pdxJ/PP_1437	pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (25/21);  new base T (3/8);  total (28/29)
Fisher's exact test for biased strand distribution p-value = 1.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCC  >  NC_002947/1637183‑1637470                                                                                                                                                |                                                                                                                                                  tATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGATTCAATGCTGCTGTGCCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATAG‑TTGGTCGCCTGTACCGGACTTTTCGCGGGCAAGCCcg                                                                                                                                             >  5:113991/1‑149 (MQ=255)      aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAag                                                                                                                                        <  5:73485/149‑1 (MQ=255)       tGAGTTGAACCTTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTAAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                       <  5:248602/149‑1 (MQ=255)       tGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCTTGCCCGCGaaaa                                                                                                                                       <  7:93395/149‑3 (MQ=255)       tGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTCTTCGCGGGCTTGCCCGCGaaaa                                                                                                                                       <  8:314554/149‑3 (MQ=255)         aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCTTGCCCGCGAaaag                                                                                                                                      <  3:140429/148‑1 (MQ=255)         aGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCCGAGTGATCGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAaaag                                                                                                                                      <  5:317772/149‑1 (MQ=25)             gAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCa                                                                                                                                  <  1:201394/148‑1 (MQ=255)             gAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAg                                                                                                                                 >  1:183594/1‑149 (MQ=255)              aaCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGccgg                                                                                                                                 >  6:73722/1‑146 (MQ=38)               aCATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGca                                                                                                                               >  3:188429/1‑149 (MQ=255)                    ggCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TCGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAggc                                                                                                                          >  2:180775/1‑149 (MQ=25)                       cATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTTTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGTTTG‑TTGGTCGCCTGTAGCGGCCTTATCGCGGGCTAGCCCGCGAAGAGGCCAGCACAggcggc                                                                                                                       <  6:113991/149‑1 (MQ=255)                        aTGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGcc                                                                                                                      >  4:211487/1‑149 (MQ=255)                           cgcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGAccaa                                                                                                                    <  5:73722/148‑3 (MQ=25)                           cgcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGa                                                                                                                    <  2:183594/148‑1 (MQ=255)                            gcTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TCGGCCGCCTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGAccaaca                                                                                                                  <  1:180775/149‑5 (MQ=21)                                     gcATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAg                                                                                                         >  4:268027/1‑149 (MQ=255)                                                 tgtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAgg                                                                                              <  3:309026/148‑1 (MQ=255)                                                  gtCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAggg                                                                                             <  8:150716/148‑1 (MQ=255)                                                          cAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                                    >  1:242889/1‑149 (MQ=255)                                                           aaGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCggg                                                                                    >  4:109841/1‑148 (MQ=255)                                                            tggCTGGTGTGGCCGAGATGAAGGCTCTGATCGTGGCGGCTTCGCGGTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCTGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTa                                                                                  <  3:192643/148‑1 (MQ=255)                                                              gctgctGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAgg                                                                                >  3:256957/1‑149 (MQ=255)                                                                     tgGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCa                                                                           <  3:109841/147‑1 (MQ=255)                                                                       gCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggg                                                                        <  2:263160/148‑1 (MQ=255)                                                                        ccGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                                      >  1:265025/1‑149 (MQ=255)                                                                         cGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                                      >  4:159733/1‑148 (MQ=255)                                                                           agaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                       <  5:271273/81‑1 (MQ=255)                                                                           agaTGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                       >  6:271273/1‑81 (MQ=255)                                                                               gAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctgc                                                               >  5:177783/1‑149 (MQ=255)                                                                                   gcgcTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                           <  4:256957/149‑1 (MQ=255)                                                                                   gcgcTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                           >  5:203766/1‑149 (MQ=255)                                                                                     gcTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCt                                                         <  3:159733/149‑1 (MQ=255)                                                                                         aTCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                     <  4:188429/149‑1 (MQ=255)                                                                                          tCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAcc                                                    <  3:268027/149‑1 (MQ=255)                                                                                                 ggcTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgggg                                                                       <  1:309265/123‑1 (MQ=255)                                                                                                 ggcTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAgggg                                                                       >  2:309265/1‑123 (MQ=255)                                                                                                    ttCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCAt                                          >  4:328879/1‑149 (MQ=255)                                                                                                       gcgcTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGc                                        >  1:326217/1‑148 (MQ=255)                                                                                                        cgcTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCaa                                      <  3:211487/149‑1 (MQ=255)                                                                                                          cTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACt                                    <  2:242889/149‑1 (MQ=255)                                                                                                             aTTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCa                                 >  3:314413/1‑149 (MQ=255)                                                                                                              ttg‑ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCAt                                >  6:183018/1‑149 (MQ=255)                                                                                                                g‑ttgGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATg                               <  6:177783/148‑1 (MQ=255)                                                                                                                     gTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGa                          >  3:269715/1‑149 (MQ=255)                                                                                                                          cTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                                      >  2:265617/1‑100 (MQ=17)                                                                                                                          cTGTGCTGGCCTCTTCGCGGGCTTGCCCGCGAAAAGGCCGGTACAGGCGACCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAggggg                                                                      <  1:265617/100‑1 (MQ=17)                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGCGGCCAGCACAGGCGGCCGATCGCTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGCCCTCGACCCGCCCGCCATGCACCTGCATGATCGCCTGCa                     >  5:261211/1‑149 (MQ=255)                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                                >  5:347330/1‑106 (MQ=255)                                                                                                                          cTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAgctg                                                                <  6:347330/106‑1 (MQ=255)                                                                                                                                cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc                      >  2:214085/1‑142 (MQ=255)                                                                                                                                cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTc                      <  1:214085/142‑1 (MQ=255)                                                                                                                                cGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCg               >  8:85827/1‑149 (MQ=255)                                                                                                                                  gCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGc              >  6:87027/1‑148 (MQ=255)                                                                                                                                     ttttCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGccaac          >  7:184066/1‑149 (MQ=255)                                                                                                                                             ggCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCcaaccccaaccc  <  2:326217/149‑1 (MQ=255)                                                                                                                                                aaGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCg               >  4:283127/1‑133 (MQ=255)                                                                                                                                                aaGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCg               <  3:283127/133‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                  TATCAATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCC  >  NC_002947/1637183‑1637470

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