Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,101,837	C→G	6.9%	A29G (GCC→GGC)	PP_3610 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,101,837	0	C	G	6.9%	216.2 / 4.0	69	A29G (GCC→GGC)	PP_3610	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (31/32);  new base G (2/3);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

GGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGA  >  NC_002947/4101692‑4101984                                                                                                                                                    |                                                                                                                                                    ggccggcGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCtt                                                                                                                                                  <  5:193594/149‑1 (MQ=255)  gccggcGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGTCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTg                                                                                                                                                 <  3:168262/149‑1 (MQ=255)  gccggcGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTg                                                                                                                                                 <  5:259196/149‑1 (MQ=255)    cggcGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGt                                                                                                                                               <  5:358607/149‑1 (MQ=255)             cATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGC‑GTTGACCAAGGCGGCTTGGTCCAccgccag                                                                                                                                     <  2:30862/149‑7 (MQ=255)                  ccGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAAACCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGCCGGAGCCGCCTTGGTCAACGCCTTGACAg                                                                                                                                 >  8:341472/1‑149 (MQ=255)                     aCAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGc                                                                                                                              <  3:319876/149‑1 (MQ=255)                            aCTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGcc                                                                                                                       <  6:133088/149‑1 (MQ=255)                            aCTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGcc                                                                                                                       <  6:215861/149‑1 (MQ=255)                              tGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCcag                                                                                                                     >  5:179985/1‑149 (MQ=255)                                cACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCcag                                                                                                                     >  3:201443/1‑147 (MQ=255)                                cACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCATGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCcagc                                                                                                                    <  4:36566/148‑1 (MQ=255)                                 aCCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGG‑TGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCcagcag                                                                                                                  >  2:104761/1‑148 (MQ=255)                                    cggCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCCACCGTCGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAt                                                                                                               <  8:345679/147‑1 (MQ=255)                                     tGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                              <  4:201443/149‑1 (MQ=255)                                     tGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                              <  4:321839/149‑1 (MQ=255)                                       ccACTGCCATAGATCCTTTGCCATGAGCGACGAACGAATGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTgg                                                                                                            <  6:256693/149‑1 (MQ=255)                                        caatGCCATAGAGCATTACCCATCAGCCACTAACGAAAGTCCCTGATCGACCTTGGCCACCTCCCAACCCGCAA‑‑ACTGGACCAATGCCTGCAGGTCGACCAAGCCGACTACGTCAACGCCAAGACAGGGCAACTGCCAGCAGCATAgg                                                                                                            >  4:340813/4‑148 (MQ=255)                                        cACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGt                                                                                                           <  1:366150/149‑1 (MQ=255)                                        cACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGt                                                                                                           <  2:234270/149‑1 (MQ=255)                                         aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAATTAGTGGCCCAATGCCTGGCGGTGGACCAAACCCCCATGCTCAACGCCTATACAGGGAAACTGCCAGCAGAATTgg                                                                                                            >  1:315003/1‑149 (MQ=255)                                         aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                              <  1:187748/145‑1 (MQ=255)                                         aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                              >  2:187748/1‑145 (MQ=255)                                         aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGt                                                                                                           >  2:26657/1‑148 (MQ=255)                                           tGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTa                                                                                                          >  3:312376/1‑147 (MQ=255)                                              cATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTagag                                                                                                       >  3:34558/1‑147 (MQ=255)                                              cATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTagag                                                                                                       <  4:34558/147‑1 (MQ=255)                                              cATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCa                                                                                                     >  8:50002/1‑149 (MQ=255)                                               aTAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCAc                                                                                                    >  7:10303/1‑149 (MQ=255)                                                       ttAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGc                                                                                            >  8:1639/1‑149 (MQ=255)                                                         aGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGc                                                                                            <  1:82869/147‑1 (MQ=255)                                                           ccATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGt                                                                                        >  1:124523/1‑149 (MQ=255)                                                            cATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTg                                                                                       <  6:179985/149‑1 (MQ=255)                                                                            aaGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGc                                                                       >  2:143780/1‑149 (MQ=255)                                                                              gTCCCTGATCGACCGTAGCCACATGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCAc                                                                     >  6:239779/1‑149 (MQ=255)                                                                               tCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACt                                                                    >  7:147726/1‑149 (MQ=255)                                                                                cccTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAgtgcctg                                                                   >  6:298611/1‑143 (MQ=255)                                                                                 ccTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTATAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTgg                                                                  >  3:126143/1‑149 (MQ=255)                                                                                     agcgcccGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGCATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCt                                                              <  4:312376/144‑1 (MQ=255)                                                                                         aCCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGTGCCTGGGCCTgctc                                                          <  5:298611/149‑4 (MQ=255)                                                                                         aCCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGTGCCTGGGCCTg                                                             >  1:144688/1‑146 (MQ=255)                                                                                         aCCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGTGCCTGGGCCTg                                                             <  2:144688/146‑1 (MQ=255)                                                                                          ccGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCa                                                          >  5:241635/1‑148 (MQ=255)                                                                                                ccACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAg                                                   >  4:216112/1‑149 (MQ=255)                                                                                                 cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                  >  1:260307/1‑149 (MQ=255)                                                                                                 cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                  <  5:186113/149‑1 (MQ=255)                                                                                                 cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                  <  5:308345/149‑1 (MQ=255)                                                                                                  aCCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTg                                                 >  7:69817/1‑149 (MQ=255)                                                                                                   ccTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTgg                                                <  1:155791/149‑1 (MQ=255)                                                                                                         aaCCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCa                                          >  1:190037/1‑149 (MQ=255)                                                                                                         aaCCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCa                                          >  3:85512/1‑149 (MQ=255)                                                                                                           cccGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCTGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCAtt                                        <  1:26657/149‑1 (MQ=255)                                                                                                            ccGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTc                                                           <  4:103322/129‑1 (MQ=255)                                                                                                            ccGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTc                                                           >  3:103322/1‑129 (MQ=255)                                                                                                            ccGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTg                                       >  3:271966/1‑149 (MQ=255)                                                                                                                aa‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCt                                   <  4:85512/149‑1 (MQ=255)                                                                                                                    ccTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTgg                                 <  1:143780/149‑1 (MQ=255)                                                                                                                          cGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCNCTgg                           <  5:109891/149‑1 (MQ=255)                                                                                                                          cGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTgg                           <  2:124523/149‑1 (MQ=255)                                                                                                                           gAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACAggg                          >  7:256302/1‑149 (MQ=255)                                                                                                                              tGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCaa                       <  7:50002/149‑1 (MQ=255)                                                                                                                                  tGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTgg                   >  4:5503/1‑149 (MQ=255)                                                                                                                                       gttgACCAAGGCGGCTTGGTCCACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTg                >  1:254770/4‑147 (MQ=255)                                                                                                                                       gttgACCAAGGCGGCTTGGTCCACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTg                <  2:254770/144‑1 (MQ=255)                                                                                                                                       gttgACCAAGGCGGCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTg                                                 <  3:189823/111‑1 (MQ=255)                                                                                                                                       gttgACCAAGGCGGCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTg                                                 >  4:189823/4‑114 (MQ=255)                                                                                                                                            ccAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGc                                                                <  3:66925/94‑1 (MQ=255)                                                                                                                                            ccAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGc                                                                >  4:66925/1‑94 (MQ=255)                                                                                                                                            ccAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAg         >  7:114111/1‑149 (MQ=255)                                                                                                                                                 ccgccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAgccgcc    >  3:107958/1‑149 (MQ=255)                                                                                                                                                  cgccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGAACCTCGCTGCCAGTGAGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGAACTGGGAAAGTGGCTGGAGCAAgccgac    >  4:128924/1‑146 (MQ=255)                                                                                                                                                   gccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTc                                                           >  7:344360/1‑92 (MQ=255)                                                                                                                                                   gccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTc                                                           <  8:344360/92‑1 (MQ=255)                                                                                                                                                   gccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAgccgcc    >  3:139708/1‑147 (MQ=255)                                                                                                                                                   gccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGa  <  8:147726/149‑1 (MQ=255)                                                                                                                                                   gccTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGa  <  4:271966/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCCGA  >  NC_002947/4101692‑4101984

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