Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,293,967	C→G	6.6%	D13H (GAC→CAC)	PP_3763 ←	precorrin‑6A synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,293,967	0	C	G	6.6%	190.4 / 3.5	61	D13H (GAC→CAC)	PP_3763	precorrin‑6A synthase
Reads supporting (aligned to +/- strand):  ref base C (22/35);  new base G (2/2);  total (24/37)
Fisher's exact test for biased strand distribution p-value = 6.43e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.94e-01

GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCA  >  NC_002947/4293821‑4294113                                                                                                                                                    |                                                                                                                                                   gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCAc                                                                                                                                                  >  2:340499/1‑149 (MQ=255) gCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCAAAGGTGACCTGGCGCGGG‑TCAc                                                                                                                                                  >  5:231040/1‑149 (MQ=255)       accaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCgcgc                                                                                                                                            <  7:25167/149‑1 (MQ=255)       accaccCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCgcgc                                                                                                                                            <  1:37740/149‑1 (MQ=255)           cccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TGACCCGCGccagg                                                                                                                                        <  1:20118/149‑4 (MQ=255)           cccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGAt                                                                                                                                        <  3:353999/149‑1 (MQ=255)            ccTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATa                                                                                                                                       <  2:273301/149‑1 (MQ=255)             cTCGGGGCGGTAGCGGTGCATAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑ACACCCGCGCCGATAc                                                                                                                                      <  5:223277/149‑1 (MQ=255)             cTCGGGGCGGTAGCGCTGCAGAATGGTCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATAc                                                                                                                                      <  4:94547/149‑1 (MQ=255)             cTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TGACCcg                                                                                                                                               <  6:24504/140‑1 (MQ=255)             cTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TGACCcg                                                                                                                                               >  5:24504/1‑140 (MQ=255)                       tAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGA‑TCACCCGCGCCGATACCGATcaacaa                                                                                                                            <  2:31736/149‑1 (MQ=255)                           gcTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAgg                                                                                                                        <  2:31930/149‑1 (MQ=255)                            cTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATca                                                                                                                                >  5:314545/1‑140 (MQ=255)                            cTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATca                                                                                                                                <  6:314545/140‑1 (MQ=255)                               cAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCCCCCGCGGCGCTACCGCTGAACCACAGGTCtt                                                                                                                    >  1:205904/1‑149 (MQ=255)                                  aaTGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCa                                                                                                                 >  3:40704/1‑149 (MQ=255)                                    tGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCAtt                                                                                                               <  4:173872/149‑1 (MQ=255)                                      gCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCGTTCATtgc                                                                                                             >  8:337527/1‑149 (MQ=255)                                       ccTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgcc                                                                                                            <  5:122358/149‑1 (MQ=255)                                       ccTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgcc                                                                                                            >  7:184661/1‑149 (MQ=255)                                       ccTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgcc                                                                                                            <  1:157295/149‑1 (MQ=255)                                            cgcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccctgcc                                                                                                       <  1:332888/149‑1 (MQ=255)                                             gcAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATtgccctgccc                                                                                                      >  4:337054/1‑149 (MQ=255)                                                  cgcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTCTGc                                                                                                 <  8:38407/149‑1 (MQ=255)                                                   gcACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGGCTCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGc                                                                                                 <  1:284656/149‑1 (MQ=255)                                                    cACCAGTTCATCTTTGTCGACGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCgg                                                                                               <  2:180013/149‑1 (MQ=255)                                                        aGTTCATCTTTGTCGTCGCGCTTGTCGATCACGAAGAACACGCTGGCACTGCGCAGCGCATCTACCGATTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGTTCTTTCATTGCCCTGCCCTTTGCGGGAAc                                                                                           <  4:279722/149‑1 (MQ=255)                                                         gTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAAcc                                                                                          <  3:337054/149‑1 (MQ=255)                                                             aTCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCcgc                                                                                       <  4:144382/148‑1 (MQ=255)                                                                  tGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     <  7:337527/145‑1 (MQ=255)                                                                  tGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     <  4:40704/145‑1 (MQ=255)                                                                           ccTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCg                                                                        >  5:34703/1‑149 (MQ=255)                                                                               gTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGt                                                                      <  1:340499/147‑1 (MQ=255)                                                                                  gAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCg                                                                 >  5:299749/1‑149 (MQ=255)                                                                                    gCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCgg                                                                >  2:281755/1‑148 (MQ=255)                                                                                    gCACGAAGAACAAGCTGGCACTGCGCAGCGCATCGACCGATTCATAGGTGACCTGGCGCGTG‑TCACCCGATCCGATAAAGATCAACAACAGGTATTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGCCGCTTAACgg                                                                >  6:259313/1‑148 (MQ=255)                                                                                        gaagaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTCCCCTTTGCGGGAACCCGCGACCATTATCTGGCGGTTACCGGTAttt                                                           >  6:180380/1‑149 (MQ=255)                                                                                           gaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTgg                                                         >  5:60354/1‑148 (MQ=255)                                                                                           gaaCACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTgg                                                         >  7:317708/1‑148 (MQ=255)                                                                                                cgcgtgCACTGCGCAGCGCATCGACCGCTTCATATGTGACCTGGCGCTGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAagc                                                   <  2:317166/144‑1 (MQ=255)                                                                                                 gCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     >  8:263571/1‑114 (MQ=255)                                                                                                 gCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGa                                                                                     <  7:263571/114‑1 (MQ=255)                                                                                                 gCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAag                                                    >  7:111813/1‑147 (MQ=255)                                                                                                  gtgcgccTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGCAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAagc                                                   <  1:367925/141‑1 (MQ=255)                                                                                                    ggCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgct                                               >  1:32975/1‑149 (MQ=255)                                                                                                     gCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAgctg                                              >  5:45178/1‑149 (MQ=255)                                                                                                           cgcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                         <  1:281755/148‑1 (MQ=255)                                                                                                            gcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACACCAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTagaa                                       <  6:231040/149‑1 (MQ=255)                                                                                                            gcagcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTagaa                                       <  8:111813/149‑1 (MQ=255)                                                                                                               gcgcaTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTag                                         >  3:236155/1‑144 (MQ=255)                                                                                                                  caTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCga                                 <  8:184661/149‑1 (MQ=255)                                                                                                                   aTCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCgag                                <  4:362880/149‑1 (MQ=255)                                                                                                                      gACCGCTTCATAGGTGAGCTGGCGCGTG‑GCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGTGTACCCTCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTTGAAGATCGAGaa                              <  5:180380/148‑1 (MQ=255)                                                                                                                        ccGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGtt                           <  6:34703/149‑1 (MQ=255)                                                                                                                         agCTTCATAGGTTACCTGGCGCGGG‑TCACCCGCTCCGATACCGATCAACATCAGGTGTTTCATTGCCCTGCCCTTTGATGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGATTAGAAGATCGAGAAGTTg                          <  5:259313/148‑1 (MQ=255)                                                                                                                                      aCCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCACGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGa             >  4:44089/1‑149 (MQ=255)                                                                                                                                        cTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGAcc           <  6:299749/149‑1 (MQ=255)                                                                                                                                           gcgcGGG‑TGACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCcgc        >  7:351/1‑149 (MQ=255)                                                                                                                                            cgcgGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCcg         >  2:301449/1‑147 (MQ=255)                                                                                                                                                 g‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCa  <  5:84832/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GCGATAACCACCCTCGGGGCGGTAGCGCTGCAGAATGGCCTTGCGCAGGCGCACCAGTTCATCTTTGTCGTCGCCCTTGTCGAGCACGAAGAACACGCTGGCACTGCGCAGCGCATCGACCGCTTCATAGGTGACCTGGCGCGGG‑TCACCCGCGCCGATACCGATCAACAACAGGTCTTTCATTGCCCTGCCCTTTGCGGGAACCCGCGACCATTGTCGGGCGGTTACCGGTATTTGGCAAAGCAGCTGCTTAGAAGATCGAGAAGTTGTAGCTGACAATGACCCGCGTCTCA  >  NC_002947/4293821‑4294113

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