Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,101,834	C→G	6.9%	A28G (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,834	0	C	G	6.9%	220.4 / 4.0	69	A28G (GCC→GGC)	PP_3610	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (32/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.17e-01

CTTCCGGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCC  >  NC_002947/4101687‑4101982                                                                                                                                                      |                                                                                                                                                     catccGCCCGGCGGTCGCAGCGCCGACAATCGACCGCACCCTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAgcc                                                                                                                                                     <  8:139665/147‑1 (MQ=255)  ttccGGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAgccg                                                                                                                                                    >  6:358607/1‑149 (MQ=255)  ttccGGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAgccg                                                                                                                                                    >  6:342570/1‑149 (MQ=255)      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:215861/149‑1 (MQ=255)                                 aCTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGcc                                                                                                                     <  6:133088/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                                                                                                         <  2:234270/149‑1 (MQ=255)                                             cACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGt                                                                                                         <  1:366150/149‑1 (MQ=255)                                              aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAATTAGTGGCCCAATGCCTGGCGGTGGACCAAACCCCCATGCTCAACGCCTATACAGGGAAACTGCCAGCAGAATTgg                                                                                                          >  1:315003/1‑149 (MQ=255)                                              aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                            >  2:187748/1‑145 (MQ=255)                                              aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCAtt                                                                                                            <  1:187748/145‑1 (MQ=255)                                              aCTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGt                                                                                                         >  2:26657/1‑148 (MQ=255)                                                tGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTa                                                                                                        >  3:312376/1‑147 (MQ=255)                                                   cATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTagag                                                                                                     <  4:34558/147‑1 (MQ=255)                                                   cATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTagag                                                                                                     >  3:34558/1‑147 (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                                                           <  2:144688/146‑1 (MQ=255)                                                                                              aCCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGTGCCTGGGCCTg                                                           >  1:144688/1‑146 (MQ=255)                                                                                               ccGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCa                                                        >  5:241635/1‑148 (MQ=255)                                                                                                     ccACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAg                                                 >  4:216112/1‑149 (MQ=255)                                                                                                      cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                <  5:186113/149‑1 (MQ=255)                                                                                                      cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                <  5:308345/149‑1 (MQ=255)                                                                                                      cACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGt                                                >  1:260307/1‑149 (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              <  2:254770/144‑1 (MQ=255)                                                                                                                                            gttgACCAAGGCGGCTTGGTCCACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTg              >  1:254770/4‑147 (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)                                                                                                                                                      |                                                                                                                                                     CTTCCGGCCGGCGTGCGCATGGCCGACAATCGACTGCACCGTGCCACTGCCATAGAGCCTTAGCCATGAGCGACGAACGAAAGTCCCTGATCGACCGTGGCCACCTGCCAACCCGCAA‑‑CCTGGCCGAATGCCTGGCGGTGGACCAAGCCGCCTTGGTCAACGCCTTGACAGGGCAACTGCCAGCAGCATTGGTAGAGCACCTCGCTGCCAGTGCGGAGCAGGCCCAGGCACTGGGCCTGTCAAAAAAAGTGGCGGCCATTGGCCTGGCACTGGGCAAGTGGCTGGAGCAAGCCGCC  >  NC_002947/4101687‑4101982

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