Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,842,904:1	+G	19.0%	intergenic (+48/‑249)	ccoP‑II → / → PP_4259	cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,842,904	1	.	G	19.0%	145.8 / 44.3	58	intergenic (+48/‑249)	ccoP‑II/PP_4259	cbb3‑type cytochrome c oxidase subunit/iron‑sulfur cluster‑binding protein
Reads supporting (aligned to +/- strand):  ref base . (23/24);  new base G (3/8);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.36e-01

CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCATACTGCC  >  NC_002947/4842758‑4843051                                                                                                                                                      |                                                                                                                                                     ccGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCt‑t                                                                                                                                                     <  7:337678/148‑1 (MQ=255)   ggCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTa                                                                                                                                                  <  1:160328/149‑1 (MQ=255)    gCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAg                                                                                                                                                 <  6:20143/149‑1 (MQ=255)      ccAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGa                                                                                                                                               <  7:298102/149‑1 (MQ=255)        aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGc                                                                                                                                             <  5:190544/149‑1 (MQ=255)        aGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGc                                                                                                                                             <  2:156573/149‑1 (MQ=255)          gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                           <  1:144539/149‑1 (MQ=255)          gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                           <  6:44254/149‑1 (MQ=255)          gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGCgg                                                                                                                                           <  6:385619/149‑1 (MQ=255)          gAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGGCGCTCCATGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGGGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCgg                                                                                                                                           <  2:52462/149‑1 (MQ=255)             caacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGtt                                                                                                                                        >  8:238045/1‑149 (MQ=255)             caacTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGtt                                                                                                                                        <  2:83132/149‑1 (MQ=255)                cTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACCAAGGATCCT‑TGTAGGAGCGGGTTTAc                                                                                                                                     >  4:76054/1‑149 (MQ=255)                   cAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTAccc                                                                                                                                   <  2:43677/148‑1 (MQ=255)                         aaCGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCGCCTACCCAGGAGACG‑TGTTGGAGCGGGTGTAGCcccgga                                                                                                                              >  7:151428/1‑145 (MQ=255)                               aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTa                                                                                                                                                  <  2:371233/121‑1 (MQ=255)                               aaGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTa                                                                                                                                                  >  1:371233/1‑121 (MQ=255)                                      ccTTGCTGGCGGCTTCCGGCTACAGCCTGTCGCCCCCGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAC‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGc                                                                                                               <  8:285851/148‑1 (MQ=37)                                      aCTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTgg                                                                                                                >  7:285851/1‑148 (MQ=255)                                          gCTGGCGGCTTACGACTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCACCTACACAGCTTCCT‑TGTAGGATCGGGTTTAACCCCGAAGAGGCCGGTACTGGCACAg                                                                                                           >  8:347871/1‑149 (MQ=38)                                              gcggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCa                                                                                                         <  5:308881/147‑1 (MQ=255)                                               cggcTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCAttt                                                                                                      <  7:238045/149‑1 (MQ=255)                                                   ttACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCg                                                                                                  >  3:237310/1‑149 (MQ=255)                                                      cttCTACCGCCGGTCGCACCAGGCTGAGCCGGCCACGGCTTAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGCATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCa                                                                                                <  8:12288/146‑1 (MQ=37)                                                          tACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCCGGTAAA‑‑ACCGCTCCTAAACAGGATCCG‑TGTAGGAGCGGGTTTTACCGCGAAGAGGCCCGTTCTGGCACAGCATTTCTCGCAGGc                                                                                             >  8:22575/1‑147 (MQ=255)                                                                cTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑ACCCCAACTACACAACATGCA‑TCTAGGAGCCGGATTACCCCCGAAAAAGCCCGTACAGGCAAAGCCATTCACGCACCCAAGTGaa                                                                                      >  6:342829/1‑147 (MQ=255)                                                                  gTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAACCCCCGCTACTACACACGATCCG‑TGTAGGAGCGGGTTTACCCGCGCATCGGCCGGTACTGGCACAGCATTTCTAGCAGGCAAGTGAcc                                                                                     >  1:107720/1‑149 (MQ=38)                                                                   tCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTgg                                                                                  <  4:205665/149‑1 (MQ=255)                                                                         cAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAt                                                                            <  2:112423/149‑1 (MQ=255)                                                                            ggtGAGCCGGCCATGGCTGAGTATGGCAGGCGTGGCCTGTTCGCGGTTAAA‑‑CCCGCTCCTACACTGGTTCCT‑TTTAGGTGGGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAAtt                                                                           <  5:342829/145‑1 (MQ=37)                                                                                       cAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTGCTACCCAGGATGCT‑TGGTGGGGCGGGTTTACCCGCGGAGAGGCCGGGTCGGGCCCAGCATTTGTGGCGGGGAAGTGGCCTGGATGAATTGAGACCCGCCAt                                                               >  4:82282/1‑148 (MQ=21)                                                                                             tGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTTAA‑‑ACCCCGGCGTCCCAGGAGCAG‑TGGAGGGGCGGGTTTACCCCGGAAACGGCCCGGACAGGGGGAGGAATTGTCGCAGGGAAGTGAAGTGGGTCGAGTGGCGCCCGCCGTaaaagct                                                        >  8:23546/1‑144 (MQ=255)                                                                                               aGTACGGCAGGCCTGGCCTCTTCGCGGGTAAC‑‑CCCGCTCCTACACCGGATCCT‑TGTAGGAGCGGTTTAACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGtt                                                       <  1:292073/148‑1 (MQ=38)                                                                                                  aaGGCAGGCCTGGCCTCTTCGCGGGACAC‑‑CCCGCTCCTACACAGCATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCa                                                    <  6:330736/148‑1 (MQ=255)                                                                                                   aggcaggcCTGGCCTCTTCGCGGGTAAA‑‑CCCCCTACTACACAGGATCCT‑TGTATGAGCGGGGTTAACCGCGAAGAGGCCGGTACTGGGACAGCATTTGTGCCAGGCAAGTGACCAGGAAAAATTTAAACCCCCCAAAAAAAGTTTCATa                                                  >  6:151925/1‑149 (MQ=21)                                                                                                          cctggcctGTTCGCGGGTAAA‑‑CCCGCTCCTACACCGGATCCT‑TGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                           <  7:22575/149‑1 (MQ=255)                                                                                                          cctggcctCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGaa                                           <  4:237310/149‑1 (MQ=255)                                                                                                            tggcctCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAAcc                                         >  2:255612/1‑149 (MQ=255)                                                                                                                   ttCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGTAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑cg                                   <  3:76054/147‑1 (MQ=255)                                                                                                                        gggTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTa                             >  1:242485/1‑148 (MQ=255)                                                                                                                           tAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAgg                         >  2:193047/1‑149 (MQ=255)                                                                                                                           tAAA‑‑CCCGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAgg                         >  2:84828/1‑149 (MQ=255)                                                                                                                            aaa‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGt                        >  6:287374/1‑149 (MQ=255)                                                                                                                              a‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTTGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGa                                <  7:76847/139‑1 (MQ=255)                                                                                                                              a‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGa                                >  8:76847/1‑139 (MQ=255)                                                                                                                                  ccGCTCCTACA‑AGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCg                      >  5:45514/1‑147 (MQ=255)                                                                                                                                       ccTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCg               <  2:242485/149‑1 (MQ=255)                                                                                                                                        cTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCg               >  2:41876/1‑148 (MQ=255)                                                                                                                                         tACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGaa             >  5:356291/1‑149 (MQ=255)                                                                                                                                           acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                               <  6:190345/112‑1 (MQ=255)                                                                                                                                           acaaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCa                                               >  5:190345/4‑115 (MQ=255)                                                                                                                                            caaGGATCCTGTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGaa             <  1:84828/145‑1 (MQ=255)                                                                                                                                             caGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAACCGCGACTAAAGGTCGCAGTGCGAAccc          <  1:193047/149‑1 (MQ=255)                                                                                                                                             caGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAAccc          >  3:100542/1‑148 (MQ=255)                                                                                                                                              aGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCa         >  7:343987/1‑148 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTg                 >  2:48320/1‑138 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTg                 <  1:48320/138‑1 (MQ=255)                                                                                                                                                gATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCATAc      <  7:347871/149‑1 (MQ=255)                                                                                                                                                    cT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCATACTGc   <  1:41876/148‑1 (MQ=255)                                                                                                                                                     t‑tGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCATACTGcc  >  3:48729/1‑148 (MQ=255)                                                                                                                                                      |                                                                                                                                                     CCGGCCCAGGAACAACTGCAGGGCAACGACAAGGTGCACTTGCTGGCGGCTTACGTCTACAGCCTGTCGCACCAGGCTGAGCCGGCCAAGGCTGAGTAAGGCAGGCCTGGCCTCTTCGCGGGTAAA‑‑CCCGCTCCTACACAGGATCCT‑TGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGTACTGGCACAGCATTTCTCGCAGGCAAGTGACCTGGATCAATTGACACCCGCCATAACACGTTTCATACCACGAACCCAA‑CGCGACTAAAGGTCGCAGTGCGAACCCATACTGCC  >  NC_002947/4842758‑4843051

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