Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,637,336	G→A	10.7%	intergenic (+45/+21)	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,336	0	G	A	10.7%	155.3 / 7.3	56	intergenic (+45/+21)	pdxJ/PP_1437	pyridoxine 5'‑phosphate synthase/heavy metal sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (27/23);  new base A (3/3);  total (30/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCT  >  NC_002947/1637188‑1637481                                                                                                                                                      |                                                                                                                                                  aTGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAag                                                                                                                                                   <  5:73485/149‑1 (MQ=255)  tGAGTTGAACCTTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTAAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAaga                                                                                                                                                  <  5:248602/149‑1 (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                                                                      >  5:203766/1‑149 (MQ=255)                                                                              gcgcTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGAc                                                                      <  4:256957/149‑1 (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)                                                                                                                                            aGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTc         <  5:183018/149‑1 (MQ=255)                                                                                                                                              cccGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                                <  4:115773/92‑1 (MQ=255)                                                                                                                                              cccGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGAc                                                                >  3:115773/1‑92 (MQ=255)                                                                                                                                              cccGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCCACCCGGTccc       >  8:215169/1‑149 (MQ=255)                                                                                                                                               ccGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCg      >  4:291089/1‑149 (MQ=255)                                                                                                                                                   gAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCt  <  4:269715/149‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  ATGAGTTGAACATTGGCCATGCGCTGGTGGCGCATGCCTTGTTTGTCGGCTTCAAGGCTGCTGTGGCCGAGATGAAGGCGCTGATCGTGGCGGCTTCGCGCTGATTG‑TTGGTCGCCTGTACCGGCCTTTTCGCGGGCAAGCCCGCGAAGAGGCCAGCACAGGCGGCCGATCACTCGGCAGGAAACACCAGGGTAAACCGGGTAGGCCCCAGGGGGCAGCTGCTGACCTCGACCCGCCCGCCATGCAACTGCATGATCGACTTCACGATCGCCAACCCCAACCCGGTCCCGCCCT  >  NC_002947/1637188‑1637481

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