Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,448	A→G	12.3%	intergenic (‑115/+26)	hemO ← / ← PP_1006	heme oxygenase/heme receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,145,448	0	A	G	12.3%	163.1 / 13.0	65	intergenic (‑115/+26)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base A (34/23);  new base G (4/4);  total (38/27)
Fisher's exact test for biased strand distribution p-value = 7.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCC  >  NC_002947/1145301‑1145591                                                                                                                                                     |                                                                                                                                                gcgcTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCa                                                                                                                                                 >  3:148427/1‑148 (MQ=255)    cTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAaga                                                                                                                                             <  7:209547/149‑1 (MQ=255)       ggAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGC‑ACGCCCGCGAagag                                                                                                                                            <  3:182662/147‑1 (MQ=255)        gAGCGCAGGGCGGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGGCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGcc                                                                                                                                         <  5:344583/149‑1 (MQ=255)         aGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGGTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCg                                                                                                                                        <  1:332687/149‑1 (MQ=255)           ctcAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGC‑ACGCCCGCGAAGAGGCCgg                                                                                                                                       <  8:192492/146‑1 (MQ=255)              aGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTa                                                                                                                                     >  2:252644/1‑147 (MQ=255)                        ggtcgttcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGc                                                                                                                         <  4:148427/149‑1 (MQ=255)                         gtcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGcc                                                                                                                        >  1:217401/1‑149 (MQ=255)                          tcggtcATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCt                                                                                                                       >  6:204761/1‑149 (MQ=255)                               cATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGC‑ACGCCCGCGAAGAGGCCg                                                                                                                                        <  3:229962/127‑1 (MQ=255)                               cATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTTGCGGGC‑ACGCCCGCGAAGAGGCCg                                                                                                                                        >  4:229962/1‑127 (MQ=255)                                       aGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCaaa                                                                                                          >  1:159671/1‑149 (MQ=255)                                               gAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAtt                                                                                                   >  3:76604/1‑148 (MQ=255)                                                    agGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTTCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCACAGAGGCCGGTACCGAACTAACGCCTCAGCTATACCAACCCACCTTGAAAc                                                                                              >  6:37633/1‑148 (MQ=255)                                                            gaGACACAAGACGAAACACCCGGAGCCACCCCGAAAAAAAAGCCCCTACCTCCGATGCCAGAACCGGCCTCATCGCGCGC‑GTGCCCCCAAAGAGGCCGGTACAGACGTAACCCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                     <  6:104540/148‑1 (MQ=255)                                                             aGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgcg                                                                                    >  8:363577/1‑149 (MQ=255)                                                                 aCAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGcccc                                                                                >  4:56902/1‑149 (MQ=255)                                                                     gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGc                                                                                                                                          >  6:88776/1‑87 (MQ=255)                                                                     gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGc                                                                                                                                          <  5:88776/87‑1 (MQ=255)                                                                      aCGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GGGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCAt                                                                           <  4:14593/149‑1 (MQ=255)                                                                        gacaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGCCGTAACGCCTCAGATATCCCAACCCACATTGCTACCGAAGTTGCGCCCCGGCATgg                                                                         >  6:87275/4‑149 (MQ=255)                                                                         aaaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGt                                                                        <  5:204761/149‑1 (MQ=255)                                                                          aaCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGATGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTc                                                                       <  1:107396/149‑1 (MQ=255)                                                                           aCAGCGCGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCACATCGCCGGCGGCGCCCCCAAAACGGCCCGT‑CACACGGAACGCCGCTGCAATACCAAAACACAAAGATACCGGAGATGGGCCCCGGCAAGGt                                                                        >  2:193633/1‑147 (MQ=255)                                                                            cAGCGGGAGCGAGACAGTAAAAAAAGCCGGTCTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAg                                                                     <  5:37633/149‑1 (MQ=255)                                                                             aGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAgg                                                                    <  2:111300/149‑1 (MQ=255)                                                                               cGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCg                                                                  >  3:171902/1‑149 (MQ=255)                                                                                     cGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGAt                                                            <  1:252644/149‑1 (MQ=255)                                                                                      gagaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATg                                                           >  2:35433/1‑149 (MQ=255)                                                                                       agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                     >  6:128585/1‑122 (MQ=255)                                                                                       agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                                     <  5:128585/122‑1 (MQ=255)                                                                                          cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCTCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                        <  4:76604/148‑1 (MQ=255)                                                                                           aGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggc                                                      >  7:233078/1‑149 (MQ=255)                                                                                            gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                        <  3:327829/146‑1 (MQ=255)                                                                                            gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                                        >  4:327829/1‑146 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                    >  5:144725/1‑101 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACAt                                                                                                    <  6:144725/101‑1 (MQ=255)                                                                                             tAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGTTCGATGTtggctg                                                    >  5:38822/1‑149 (MQ=255)                                                                                               aaaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggctggc                                                  >  4:115361/1‑149 (MQ=255)                                                                                                aaaaaGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTggctggct                                                 >  6:150464/1‑149 (MQ=255)                                                                                                  ataGCCGGTATTTGCGGTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCGCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGctggctgg                                               <  1:193633/147‑1 (MQ=255)                                                                                                    aGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGt                                             <  3:56902/149‑1 (MQ=255)                                                                                                        ggTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTcaca                                         <  7:12885/149‑1 (MQ=255)                                                                                                           aTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGc                                                                                                                                          >  5:271253/1‑49 (MQ=255)                                                                                                           aTTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGc                                                                                                                                          <  6:271253/49‑1 (MQ=255)                                                                                                              tGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAgccgc                                    <  3:188961/148‑1 (MQ=255)                                                                                                                     ccAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCt                            >  3:231120/1‑149 (MQ=255)                                                                                                                      cAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTg                           >  1:49461/1‑149 (MQ=255)                                                                                                                            cGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                     >  5:391887/1‑149 (MQ=255)                                                                                                                            cGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                     >  2:77173/1‑149 (MQ=255)                                                                                                                            cGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGACCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCAGGCagtc                     >  8:56667/1‑149 (MQ=255)                                                                                                                            cGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATCGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtc                     >  5:179267/1‑149 (MQ=255)                                                                                                                              gCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcgt                   <  5:150464/149‑1 (MQ=255)                                                                                                                              gCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                    >  6:373828/1‑148 (MQ=255)                                                                                                                              gCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGAGACCCGGCATGGTCAGGAGGTCGATGTAGGCTGGCAGGGTCACAGCAGCCTCTCCCTGGCcgtcg                    >  6:226920/1‑148 (MQ=255)                                                                                                                               ccTCTTTGCGGGC‑ACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                    >  2:197366/1‑147 (MQ=255)                                                                                                                               ccTCTTTGCGGGC‑ACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa                  <  1:197366/149‑1 (MQ=255)                                                                                                                               ccTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCcgtcg                    >  7:92227/1‑147 (MQ=255)                                                                                                                               ccTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTa                  <  7:363577/149‑1 (MQ=255)                                                                                                                                ctctTTGCGGGC‑ACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTCATACAGAAGTTGCGCCACGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACCGCCGCCTCACCCTGTCCGTCGTac                 >  6:169714/1‑148 (MQ=255)                                                                                                                                 tctTTGCGGGC‑ACGCCCGCGAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGc                >  4:80898/1‑149 (MQ=255)                                                                                                                                 tctACGAGGCC‑GTCCCCGCAAAGCGGCCGGTACACACGTCACGCCTCAGATATCCCAAACCACAAACATAACCAAGTTGCGCCACGGCATCCACACCAGCTCCATCATGGCAGGCTCCCTCACAGCCACCTCACACTGGCCGACGTac                 >  4:113622/1‑147 (MQ=255)                                                                                                                                             gTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGaa     >  4:7605/1‑149 (MQ=255)                                                                                                                                             acgCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGa      >  2:29700/3‑148 (MQ=255)                                                                                                                                                 ccGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATcc  >  6:369228/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                GCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGCTGCGCACCGAATCC  >  NC_002947/1145301‑1145591

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