Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	895,210	G→C	11.3%	intergenic (+47/+17)	PP_0777 → / ← PP_0778	glutathione peroxidase/group 1 family glycosyl transferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	895,210	0	G	C	11.3%	181.0 / 13.6	62	intergenic (+47/+17)	PP_0777/PP_0778	glutathione peroxidase/group 1 family glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (24/31);  new base C (2/5);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 6.89e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

TTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGC  >  NC_002947/895062‑895355                                                                                                                                                      |                                                                                                                                                  tttCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCcgcgc                                                                                                                                                  >  2:1073622/1‑150 (MQ=255)    cTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCgcggg                                                                                                                                               >  1:1279382/1‑150 (MQ=255)      ggtggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCgcgggcg                                                                                                                                             <  1:714165/150‑1 (MQ=255)        tggGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGGAGGCCATCGAGAAAGCCCTAGCCTGAGGAGTGCCTGACGTGATGTTTTTGCACCTGTG‑AGATCGAGTGCCGCgcgggcgac                                                                                                                                           >  2:324405/1‑148 (MQ=255)             aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTgg                                                                                                                                      >  2:410702/1‑150 (MQ=255)             aaGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCCAGCGCCGCCCGcgcggcactcg                                                                                                                                      <  1:94316/150‑6 (MQ=255)               ggATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCCCGCGCGGCACTCGAt                                                                                                                                    <  2:237891/150‑1 (MQ=255)               ggATGGCCGGGTGCTGGCGCGTTTCTCGACGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCCAGCGCCGCCCGCGCGGCACTCGAt                                                                                                                                    <  2:1114090/150‑1 (MQ=38)                gATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGAtc                                                                                                                                   <  1:123169/150‑1 (MQ=255)                   ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                <  1:147952/150‑1 (MQ=255)                   ggCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCATGAAGAGTGCCTGACGTGTTTTTGTTGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTca                                                                                                                                <  1:666220/150‑1 (MQ=255)                     ccGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcaca                                                                                                                              <  1:308008/150‑1 (MQ=255)                          tGCTGGCGCGGTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTTGCATCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTTGATCTCACAGgcgc                                                                                                                         <  1:176512/150‑1 (MQ=255)                            cTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGATGCAGGAAATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTg                                                                                                                       >  1:1055597/1‑150 (MQ=255)                             tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                      >  1:745349/1‑150 (MQ=255)                             tGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGc                                                                                                                      >  1:650997/1‑150 (MQ=255)                                 gcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCCCGcgcggcactcga                                                                                                                                     >  1:1217153/1‑125 (MQ=255)                                 gcgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCCCGcgcggcactcga                                                                                                                                     <  2:1217153/131‑7 (MQ=255)                                  cgTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACa                                                                                                                 <  1:1132022/150‑1 (MQ=255)                                        tcGCCGCGCACGGCGCCGGATGACCCAGAGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcgc                                                                                                           >  1:954492/1‑150 (MQ=255)                                                    gcgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCa                                                                                                 >  2:432780/1‑148 (MQ=255)                                                    gcgcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCa                                                                                                 <  1:432780/148‑1 (MQ=255)                                                      gcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAg                                                                                                                             <  1:356710/118‑1 (MQ=255)                                                      gcCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAg                                                                                                                             >  2:356710/1‑118 (MQ=255)                                                       ccGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCCCCTGTG‑AGATCGAGTGCGCCGCGGGCGGAGCTGAATATCGTTGTCTAGGGAAAGATTGTGTAGAACGCCCAGCCGt                                                                                            >  2:654676/1‑150 (MQ=255)                                                         ggATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                          >  1:275325/1‑150 (MQ=255)                                                          gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGc                                                                                         >  2:1157850/1‑150 (MQ=255)                                                          gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGc                                                                                         >  1:552952/1‑150 (MQ=255)                                                          gATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGc                                                                                         <  2:1123724/150‑1 (MQ=255)                                                              aCCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCa                                                                                     <  1:445529/150‑1 (MQ=255)                                                               cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                    <  2:745349/150‑1 (MQ=255)                                                               cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                    >  1:840935/1‑150 (MQ=255)                                                               cccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAg                                                                                    <  1:1157850/150‑1 (MQ=255)                                                                ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTTGCCTGAAGAGTGCCTGACGTGATGTTTTGGCGCCTGTG‑AGATCCAGCGCCGCCCGCGCGGCACACGATCTCACAGgtgcc                                                                                                                        >  1:1225256/1‑112 (MQ=18)                                                                ccAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCGCCTGTG‑AGATCCAGCGCCGCCCGCGCGGCACTCGATCTCACAGgtgcc                                                                                                                        <  2:1225256/113‑2 (MQ=18)                                                                     ggtggtGCAGGGCATCCAGAAAGCCCTAGCCTTAATATTGCCTTACGTGATGTTTTGGCACCTCTAAACATCG‑TTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                              <  1:1185780/150‑1 (MQ=255)                                                                     ggtggtGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                              <  2:552952/150‑1 (MQ=255)                                                                                 cATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                          <  1:1149046/126‑1 (MQ=255)                                                                                 cATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAg                                                                                          >  2:1149046/1‑126 (MQ=255)                                                                                  aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                 <  1:1073622/150‑1 (MQ=255)                                                                                  aTCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGtt                                                                 <  1:34490/150‑1 (MQ=255)                                                                                    cGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCa                                                               <  2:183338/150‑1 (MQ=255)                                                                                        aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcac                                                                                                                               <  2:750707/82‑1 (MQ=255)                                                                                        aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTcac                                                                                                                               >  1:750707/1‑82 (MQ=255)                                                                                        aaaGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGTCATGCGCTGTTCAGATa                                                           >  1:65800/1‑150 (MQ=255)                                                                                           gCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATAccc                                                        <  2:897841/150‑1 (MQ=255)                                                                                                   cTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGAGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACt                                                >  1:28545/1‑150 (MQ=255)                                                                                                    tGAAGAGTGGCTGACGTGATGTTTTGGCACCTGTG‑AGACCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTg                                               >  2:857470/1‑150 (MQ=255)                                                                                                       agagTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTc                                            <  1:410702/150‑1 (MQ=255)                                                                                                        gagTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCt                                                                    >  2:1137177/1‑125 (MQ=255)                                                                                                        gagTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCt                                                                    <  1:1137177/125‑1 (MQ=255)                                                                                                                gACGTGCTGTTTTGGCACCTGTG‑AGATCGAGTGCCGTGCGGGCGGCGCTGGCTCTCAAAAGCGGTGCAAACCTTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGCTTCCCCTCGAACTGTTCGACAATCGc                                   <  1:324405/150‑1 (MQ=255)                                                                                                                  cGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCg                                 <  2:28545/150‑1 (MQ=255)                                                                                                                   gTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCgg                                >  2:837825/1‑150 (MQ=255)                                                                                                                     gATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                            <  1:1033606/72‑1 (MQ=255)                                                                                                                     gATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTgcg                                                                                                            >  2:1033606/1‑72 (MQ=255)                                                                                                                         ttttGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGcc                          >  1:796396/1‑150 (MQ=255)                                                                                                                              gCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCGGACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggc                     <  1:12280/150‑1 (MQ=255)                                                                                                                                aCCTGTG‑AGAGCGAGGGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGGAGGACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACGCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTggcgg                   <  2:167955/150‑1 (MQ=255)                                                                                                                                     tg‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                              >  2:126358/1‑86 (MQ=255)                                                                                                                                     tg‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGcc                                                                              <  1:126358/86‑1 (MQ=255)                                                                                                                                     tg‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGc              <  2:650997/150‑1 (MQ=255)                                                                                                                                        agaTCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCAt            >  2:802217/1‑150 (MQ=255)                                                                                                                                           tCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCt         <  2:1279382/150‑1 (MQ=255)                                                                                                                                           tCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCt         >  1:227967/1‑150 (MQ=255)                                                                                                                                                  ccGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATCCACCGCAAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCgcgc  <  2:954492/150‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                  TTTCTGGTGGGCAAGGATGGCCGGGTGCTGGCGCGTTTCTCGCCGCGCACGGCGCCGGATGACCCAGCGGTGGTGCAGGCCATCGAGAAAGCCCTGGCCTGAAGAGTGCCTGACGTGATGTTTTGGCACCTGTG‑AGATCGAGTGCCGCGCGGGCGGCGCTGGATCTCACAGGCGCTGCAAACATTGCGCAGAACGCCCAGCCGTAGCACCAGGTTGCCGGCATGCGCTGTTCAGATACCCCTCGAACTGTTCGACAATCGCCGGCCAGCCCTGGCGGCTGGCATGCTGTCGCGC  >  NC_002947/895062‑895355

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑