Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,761,591	C→G	6.2%	P84A (CCG→GCG)	pvdN →	pyoverdine biosynthesis‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,761,591	0	C	G	6.2%	205.8 / 3.5	64	P84A (CCG→GCG)	pvdN	pyoverdine biosynthesis‑like protein
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (2/2);  total (32/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 = 8.45e-01

CCGGCCACGGATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAGGCCTCGGGATGATCTACGGCGG  >  NC_002947/4761447‑4761739                                                                                                                                                 |                                                                                                                                                      ccGGCCACGGATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGc                                                                                                                                                   <  8:239581/149‑1 (MQ=255)    gCCACGGATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACg                                                                                                                                                >  4:185821/1‑149 (MQ=255)          gATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGt                                                                                                                                          <  1:215018/149‑1 (MQ=255)          gATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGt                                                                                                                                          <  5:64141/149‑1 (MQ=255)               cTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGt                                                                                                                                     >  7:66370/1‑149 (MQ=255)                  gAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTggg                                                                                                                                  <  5:167615/149‑1 (MQ=255)                  gAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTg                                                                                                                                    <  3:185821/147‑1 (MQ=255)                    gCGGTTTTCGCCAACTGTTCGAGCTCGCCCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGa                                                                                                                                 >  8:74821/1‑148 (MQ=255)                     cGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGac‑a                                                                                                                              <  1:75383/149‑1 (MQ=255)                        ttttCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑Acc                                                                                                                            >  1:231870/1‑148 (MQ=255)                          ttCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCa                                                                                                                          <  8:66370/148‑1 (MQ=255)                            cGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCaga                                                                                                                        <  8:289801/148‑1 (MQ=255)                              ccAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGc                                                                                                                      >  6:89975/1‑148 (MQ=255)                                 aCTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGt                                                                                                                  <  3:221886/149‑1 (MQ=255)                                    gTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTgg                                                                                                                >  1:263969/1‑148 (MQ=255)                                     ttCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTgg                                                                                                                >  7:222189/1‑147 (MQ=255)                                     ttCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGaa                                                                                                              <  5:298468/149‑1 (MQ=255)                                       cGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAg                                                                                                             <  4:116682/148‑1 (MQ=255)                                          gCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCAc                                                                                                          >  6:96948/1‑148 (MQ=255)                                            tCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACGCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACg                                                                                                         >  7:156249/1‑147 (MQ=255)                                            tCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACg                                                                                                         >  7:213601/1‑147 (MQ=255)                                                  cacACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAg                                                                                                 <  2:231870/147‑1 (MQ=255)                                                   ccACTACGCGCACTTCGCCAACTTCCTCATCACCTCGTACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAgg                                                                                                >  5:95923/1‑149 (MQ=255)                                                   ccACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGaa                                                                                                              >  2:269768/1‑135 (MQ=255)                                                   ccACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGaa                                                                                                              <  1:269768/135‑1 (MQ=255)                                                   ccACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAgg                                                                                                >  3:77944/1‑149 (MQ=255)                                                        aCGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGt                                                                                           >  2:80037/1‑149 (MQ=255)                                                         cgcgCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCg                                                                                            <  7:282394/147‑1 (MQ=255)                                                         cgcgCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGt                                                                                           <  7:74821/148‑1 (MQ=255)                                                         cgcgCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCGGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGTAGCACGAAGCCGAGGTGCGt                                                                                           <  8:233788/148‑1 (MQ=255)                                                           cgcACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGaa                                                                                        >  3:44343/1‑149 (MQ=255)                                                            gcACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAAt                                                                                       <  8:222189/149‑1 (MQ=255)                                                            gcACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAAt                                                                                       <  4:77944/149‑1 (MQ=255)                                                               cTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATggg                                                                                    <  7:233291/149‑1 (MQ=255)                                                                ttCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGt                                                                                                                  >  8:190391/1‑118 (MQ=255)                                                                  cGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCAACCAACACCCGGCACGCTGCGTCGAGTGGGac‑a                                                                                                                              <  5:194228/104‑1 (MQ=255)                                                                  cGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCAACCAACACCCGGCACGCTGCGTCGAGTGGGac‑a                                                                                                                              >  6:194228/1‑104 (MQ=255)                                                                    ccAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATgggcgggc                                                                               >  5:182731/1‑149 (MQ=255)                                                                      aaCTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGACGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCg                                                                             >  5:35058/1‑149 (MQ=255)                                                                        cTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACGCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCt                                                                           <  8:156249/149‑1 (MQ=255)                                                                         ttCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGACAACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCt                                                                           <  1:80037/149‑1 (MQ=255)                                                                              catcaCCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCtt                                                                      >  8:14879/1‑148 (MQ=255)                                                                                tcaCCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGaa                                                                   <  2:263969/149‑1 (MQ=255)                                                                                tcaCCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGaa                                                                   <  5:89975/149‑1 (MQ=255)                                                                                tcaCCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGGCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGaa                                                                   <  5:289148/149‑1 (MQ=255)                                                                                 caCCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGGAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAg                                                                  >  7:282802/1‑149 (MQ=255)                                                                                 caCCTCGCACCCGCGCCCCGTGCGCGAAGCCATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGGACCCAGCGACGACTGGGAC‑ACCCa                                                                                                                          >  1:28331/1‑93 (MQ=255)                                                                                    cTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTc                                                                >  2:244083/1‑148 (MQ=255)                                                                                       gCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGCAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAg                                                             <  8:282802/148‑1 (MQ=255)                                                                                       gCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGc                                                            >  3:148000/1‑149 (MQ=255)                                                                                         aCCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGa                                                                                                  <  6:70188/109‑1 (MQ=255)                                                                                         aCCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGa                                                                                                  >  5:70188/1‑109 (MQ=255)                                                                                          cccGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCAc                                                         <  7:14879/149‑1 (MQ=255)                                                                                                ccccGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGc                                                       >  5:135061/1‑145 (MQ=255)                                                                                                           aaGCGATCGCAAGCCTTCGCGCACAACTCGACCACCACCCTGCACGCTGCGTCGAGTGGGAG‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCAcc                                        <  5:52759/149‑1 (MQ=255)                                                                                                           aaGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCa                                          >  6:11604/1‑147 (MQ=255)                                                                                                           aaGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCa                                          <  5:11604/147‑1 (MQ=255)                                                                                                                   gCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGACTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGcac                                >  5:214679/1‑149 (MQ=255)                                                                                                                     aaGCCTTCGCGCACAACTCGACCAACACGCGGCACGCTGCATCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGcacca                              >  1:172353/1‑149 (MQ=255)                                                                                                                         cTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGa                          <  4:44343/149‑1 (MQ=255)                                                                                                                            cgcgCACAACTCGACCAACACGCGGCACGCTGCATCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAg                        <  2:172353/148‑1 (MQ=255)                                                                                                                                   aaCTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAGGCCTCgg                 >  6:302491/1‑148 (MQ=255)                                                                                                                                     cTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAGGCCTCGGgat              >  7:243272/1‑149 (MQ=255)                                                                                                                                             acacCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTCCGTGAATGGGCGGGCCGCGATCTTGAAGTCAAGCCACGCCAGATCGCACTCACCGGCCGCACCACCGAAGGCCTCGGGATGATCTAcg      >  6:58221/1‑149 (MQ=255)                                                                                                                                                cccGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAGGCCTCGGGATGATCTAcggcg   <  4:148000/149‑1 (MQ=255)                                                                                                                                                 ccGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCAACGGCAGCACCACAGAAGGCCTCGGGATGATATAcggcgg  >  6:240687/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                      CCGGCCACGGATGCCTGGAGCGGTTTTCGCCAACTGTTCGAGCTCGACCCCCACTACGCGCACTTCGCCAACTTCCTCATCACCTCGCACCCGCGCCCCGTGCGCGAAGCGATCGCAAGCCTTCGCGCACAACTCGACCAACACCCGGCACGCTGCGTCGAGTGGGAC‑ACCCAGAGCGAGTGGAAGCACGAAGCCGAGGTGCGTGAATGGGCGGGCCGCTATCTTGAAGTCAAGCCACGCCAGATCGCCCTCACCGGCAGCACCACCGAAGGCCTCGGGATGATCTACGGCGG  >  NC_002947/4761447‑4761739

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