Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,145,433	C→T	11.6%	intergenic (‑100/+41)	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,433	0	C	T	11.6%	184.1 / 11.6	69	intergenic (‑100/+41)	hemO/PP_1006	heme oxygenase/heme receptor
Reads supporting (aligned to +/- strand):  ref base C (33/28);  new base T (4/4);  total (37/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.78e-01

GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145286‑1145577                                                                                                                                                    |                                                                                                                                                  ggTGATATGGTTCAGGAGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTc                                                                                                                                                   <  8:128915/148‑1 (MQ=255)    gATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCgg                                                                                                                                               >  3:14593/1‑149 (MQ=255)    gATCTGGTTCAGGCGCGGGGCGCGCAGGGCGGGGCGGTCGGGCATGGGGGAGTCCATGGATAGAGGGCTTCTAGATACAAGACGAACCAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCgg                                                                                                                                               <  6:98073/149‑1 (MQ=255)          gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGc                                                                                                                                             <  3:274786/145‑1 (MQ=255)          gTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑Gt                                                                                                                                          <  2:165723/147‑1 (MQ=255)             cAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGccc                                                                                                                                      <  5:160218/148‑1 (MQ=255)             cAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTAGCCAGTACCGGCCTCTTCGCGGGC‑GTGccc                                                                                                                                      <  1:36613/148‑1 (MQ=255)                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                                                                                                                            <  5:88776/87‑1 (MQ=255)                                                                                    gACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGc                                                                                                                            >  6:88776/1‑87 (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                                                                       <  5:128585/122‑1 (MQ=255)                                                                                                      agaCAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTgc                                                                       >  6:128585/1‑122 (MQ=255)                                                                                                         cAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCTCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                          <  4:76604/148‑1 (MQ=255)                                                                                                          aGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtggc                                        >  7:233078/1‑149 (MQ=255)                                                                                                           gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                          >  4:327829/1‑146 (MQ=255)                                                                                                           gTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTtg                                          <  3:327829/146‑1 (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)                                                                                                                                                    |                                                                                                                                                  GGTGATCTGGTTCAGGCGCTGGGAGCGCAGGGCTGGGCGGTCGGTCATGGGGGAGTCCTTGGAAAGAGGGCGTCTAGATACAAGACGAAACAGCGGGAGCGAGACAGTAAAAAAAGCCGGTATTTGCGTTGCCAGTACCGGCCTCTTCGCGGGC‑GTGCCCGCAAAGAGGCCGGTACAGACGTAACGCCTCAGATATCCCAAACCACATTGATACCGAAGTTGCGCCCCGGCATGGTCAGGCGGTCGATGTTGGCTGGCTGGGTCACAGCCGCCTCACCCTGGCCGTCGTAGC  >  NC_002947/1145286‑1145577

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