Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,229,201:1	+C	6.3%	intergenic (‑143/+33)	PP_4606 ← / ← PP_4607	ferric siderophore receptor/transmembrane sensor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,229,201	1	.	C	6.3%	258.7 / 11.3	63	intergenic (‑143/+33)	PP_4606/PP_4607	ferric siderophore receptor/transmembrane sensor
Reads supporting (aligned to +/- strand):  ref base . (31/28);  new base C (2/2);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.98e-01

TCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGC  >  NC_002947/5229055‑5229349                                                                                                                                                    |                                                                                                                                                     tCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑gg                                                                                                                                                    >  7:18812/1‑149 (MQ=255)  cATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACggg‑ggg                                                                                                                                                   >  1:111250/1‑149 (MQ=255)   aTGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGt                                                                                                                                                  <  4:64775/149‑1 (MQ=255)      gACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCC‑CCCGTAGGGGCGGGttt                                                                                                                                                <  3:170551/148‑1 (MQ=39)        cGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTTTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAcc                                                                                                                                             <  7:46876/149‑1 (MQ=255)         gCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTAccc                                                                                                                                            >  1:89325/1‑149 (MQ=255)           aGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGt                                                                                                                                          >  5:2277/1‑149 (MQ=255)                  tgAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAaga                                                                                                                                     <  2:248481/147‑1 (MQ=255)                          tGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGcaca                                                                                                                           <  1:196746/149‑1 (MQ=255)                             ttattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGc                                                                                                                        <  6:211031/149‑1 (MQ=255)                              tattGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGc                                                                                                                        >  3:208735/1‑148 (MQ=255)                                    gAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAcc                                                                                                                 <  5:219705/149‑1 (MQ=255)                                     aatcaatcGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTAcc                                                                                                                 >  6:94196/1‑148 (MQ=255)                                             gCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAggg                                                                                                        >  6:52758/1‑149 (MQ=255)                                             gCACTTTCACTCATTACACGGACGCCCCGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                         <  2:233289/148‑1 (MQ=255)                                               aCTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGt                                                                                                       >  1:246616/1‑148 (MQ=255)                                                     aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTTGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCTCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                         <  1:32684/124‑1 (MQ=255)                                                     aCTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAgg                                                                                                                         >  2:32684/1‑124 (MQ=255)                                                      cTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGcc                                                                                               >  6:148631/1‑149 (MQ=255)                                                          ttACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCaa                                                                                           >  1:246473/1‑149 (MQ=255)                                                                ggATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCt                                                                                                                                 <  1:217566/105‑1 (MQ=255)                                                                ggATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCt                                                                                                                                 >  2:217566/1‑105 (MQ=255)                                                                 gATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTc                                                                                    >  1:208227/1‑149 (MQ=255)                                                                     ccACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGcc                                                                                <  2:208227/149‑1 (MQ=255)                                                                      cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCt                                                                                   <  4:135798/145‑1 (MQ=255)                                                                      cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCt                                                                                   >  3:135798/1‑145 (MQ=255)                                                                      cACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                               <  2:111250/149‑1 (MQ=255)                                                                       aCGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGt                                                                              <  2:246616/149‑1 (MQ=255)                                                                              gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                         >  2:78397/1‑115 (MQ=255)                                                                              gTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAgg                                                                                                         <  1:78397/115‑1 (MQ=255)                                                                               ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAcc                                                                      >  5:153805/1‑149 (MQ=255)                                                                               ttACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAc                                                                       <  2:246473/148‑1 (MQ=255)                                                                                 aCCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAAccc                                                                     <  6:181842/148‑1 (MQ=255)                                                                                     cACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCAc                                                                >  3:251915/1‑149 (MQ=255)                                                                                             ccTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGccc                                                         <  8:18812/148‑1 (MQ=255)                                                                                                 gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                      <  4:251915/147‑1 (MQ=255)                                                                                                 gAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGc                                                      <  5:94196/147‑1 (MQ=255)                                                                                                  aaaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa                                                   >  1:127859/1‑149 (MQ=255)                                                                                                  aaaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa                                                   <  5:52758/149‑1 (MQ=255)                                                                                                  aaaTATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGaa                                                   >  5:218241/1‑149 (MQ=255)                                                                                                         ttgccTGTGCCAGCCTCTTCACGGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTACCGAACCAGCGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcg                                            >  6:59771/5‑149 (MQ=18)                                                                                                          ttGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAgcgc                                           >  5:137656/1‑149 (MQ=255)                                                                                                                  gcggcCGCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCg                                   <  7:169584/149‑1 (MQ=255)                                                                                                                   cggcCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCgg                                  <  2:127859/149‑1 (MQ=255)                                                                                                                           tcacgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                               <  5:249024/93‑1 (MQ=21)                                                                                                                           tcacgGGTAAACCC‑CCCGTAGGGGCGGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCg                                                                               >  6:249024/4‑96 (MQ=21)                                                                                                                           tCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGaa                                                                                                              <  8:49362/65‑1 (MQ=255)                                                                                                                           tCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGaa                                                                                                              >  7:49362/1‑65 (MQ=255)                                                                                                                            cgcgGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGc                            <  2:89325/146‑1 (MQ=255)                                                                                                                                ggTAAACCCGCCCCTACGGG‑GGGTTTACGCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                      >  1:153629/1‑148 (MQ=255)                                                                                                                                ggTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCgg                      >  5:134419/1‑148 (MQ=255)                                                                                                                                 gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                     >  6:31907/1‑148 (MQ=255)                                                                                                                                 gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                     >  6:35870/1‑148 (MQ=255)                                                                                                                                 gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGa                    <  5:148631/149‑1 (MQ=255)                                                                                                                                 gTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTGGCCCAGCGCATCCAGCGCCCGCAGCGGCTCGGGCAGCggg                     >  8:240112/1‑148 (MQ=255)                                                                                                                                  tAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCggg                     >  2:11548/1‑147 (MQ=255)                                                                                                                                    aaCCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGaa                   <  6:260752/147‑1 (MQ=255)                                                                                                                                    aaCCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGaa                   >  5:260752/1‑147 (MQ=255)                                                                                                                                      cccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTg               >  3:152578/1‑149 (MQ=255)                                                                                                                                       ccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc              >  4:141539/1‑149 (MQ=255)                                                                                                                                       ccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc              >  5:121142/1‑149 (MQ=255)                                                                                                                                       ccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc              <  8:1152/149‑1 (MQ=255)                                                                                                                                       ccGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGc              >  8:250540/1‑149 (MQ=255)                                                                                                                                         gCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCg            >  4:244190/1‑149 (MQ=255)                                                                                                                                             cTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGa        <  2:153629/149‑1 (MQ=255)                                                                                                                                                   g‑gggTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCAAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTgc  <  6:153805/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TCATGGACGCAGACCTGTGAGGGGGTGATTATTGCGAATCAATCGCACTTTCACTCATTACACGGATGCCACGTGCTGTTACCTCACCTGTGCCTGGAAATATTTTTGCTGTGGCGGCCTCTTCGCGGGTAAACCCGCCCCTACGGG‑GGGTTTACCCGTGAAGAGGCTGGCACAGGCAACTACCGAACCAGGGTCAAACGCCCCAACACGGTCTGCCGTGAAAAACCCATCACCTTGCCCAGCGAATCCAGCGCCAGCAGCGGCTCGGCCACCGGGAAGCTGCCGCTGACCTTGC  >  NC_002947/5229055‑5229349

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