Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,525,375	G→A	5.7%	G336D (GGC→GAC)	murC →	UDP‑N‑acetylmuramate‑‑L‑alanine ligase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,525,375	0	G	A	5.7%	212.6 / 2.3	70	G336D (GGC→GAC)	murC	UDP‑N‑acetylmuramate‑‑L‑alanine ligase
Reads supporting (aligned to +/- strand):  ref base G (34/32);  new base A (2/2);  total (36/34)
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.41e-01

AACCACAACGTGCTCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGCACCCGCGACTTGTAC  >  NC_002947/1525230‑1525517                                                                                                                                                   |                                                                                                                                               aaCCACAACGTGCTCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAg                                                                                                                                             >  5:233632/1‑149 (MQ=255)     acaACGTGCTCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCAGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAgtgt                                                                                                                                          >  7:266884/1‑148 (MQ=255)           tGCTCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTg                                                                                                                                   <  3:37477/149‑1 (MQ=255)              tCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGACAGTGTGATGCTGGTc                                                                                                                                >  7:178046/1‑149 (MQ=255)              tCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGACAGTGTGATGCTGGTc                                                                                                                                >  7:172632/1‑149 (MQ=255)               cAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTc                                                                                                                                >  3:164778/1‑148 (MQ=255)                aaCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGa                                                                                                                                        <  5:97574/139‑1 (MQ=255)                aaCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGa                                                                                                                                        >  6:97574/1‑139 (MQ=255)                aaCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcgc                                                                                                                              >  5:95476/1‑148 (MQ=255)                aaCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcga                                                                                                                              >  3:210234/1‑149 (MQ=255)                aaCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcga                                                                                                                              >  8:375525/1‑149 (MQ=255)                 aCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcgac                                                                                                                             >  1:174053/1‑149 (MQ=255)                 aCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcgac                                                                                                                             <  6:382709/149‑1 (MQ=255)                 aCGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTcgac                                                                                                                             >  1:47950/1‑149 (MQ=255)                         gCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACg                                                                                                                      <  5:204131/148‑1 (MQ=255)                         gCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACg                                                                                                                      <  4:210234/148‑1 (MQ=255)                                   ccatcgccaCTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGGACTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCAccc                                                                                                            >  1:318736/1‑149 (MQ=255)                                          caCTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAg                                                                                                    >  7:45928/1‑149 (MQ=255)                                             tGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCAACGAAGTgg                                                                                                 >  1:141416/1‑149 (MQ=255)                                                   atgCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCGTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGcggcgg                                                                                           <  2:329192/147‑1 (MQ=255)                                                   aGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGcggcgg                                                                                           <  2:174053/149‑1 (MQ=255)                                                   aGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGcggcgg                                                                                           <  3:395505/149‑1 (MQ=255)                                                         tACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATc                                                                                      >  7:93092/1‑148 (MQ=255)                                                         tACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGACGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATc                                                                                      >  3:44159/1‑148 (MQ=255)                                                              aCGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCt                                                                                >  5:151363/1‑149 (MQ=255)                                                                gAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCtgt                                                                              <  2:215528/149‑1 (MQ=255)                                                                    ccATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTgcgc                                                                          <  2:47950/149‑1 (MQ=255)                                                                          tCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGtggc                                                                    <  4:44159/149‑1 (MQ=255)                                                                          tCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGtggc                                                                    <  8:93092/149‑1 (MQ=255)                                                                            cAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGtggctg                                                                  <  6:95476/149‑1 (MQ=255)                                                                               ggCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGtggctggc                                                                >  5:325118/1‑148 (MQ=255)                                                                               ggCCTGTCCGG‑CTTCCAGGGGGTCGGTCGTCGCTTCCAGGTATACTGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGcc                                                               <  8:371358/149‑1 (MQ=255)                                                                                  cTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTGCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAg                                                            <  1:80148/149‑1 (MQ=255)                                                                                   tGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTggcgg                                                                                              >  1:187030/1‑114 (MQ=255)                                                                                   tGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTggcgg                                                                                              <  2:187030/114‑1 (MQ=255)                                                                                    gTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAgcc                                                          >  5:369702/1‑149 (MQ=255)                                                                                     tCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAgc                                                           <  2:141416/147‑1 (MQ=255)                                                                                       cGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAgccgc                                                        <  8:266884/148‑1 (MQ=255)                                                                                            ttCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCt                                                   >  1:282879/1‑149 (MQ=255)                                                                                               cAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGt                                                <  6:325118/149‑1 (MQ=255)                                                                                                 gggTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGa                                              <  2:282879/149‑1 (MQ=255)                                                                                                   gtgtCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTggcgg                                                                                              <  8:185882/99‑1 (MQ=255)                                                                                                   gtgtCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTggcgg                                                                                              >  7:185882/1‑99 (MQ=255)                                                                                                        ggtcgTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGTGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTa                                       >  5:26571/1‑149 (MQ=255)                                                                                                         gtcgtcGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTAc                                      >  1:329617/1‑149 (MQ=255)                                                                                                          tcgtcgCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTAcc                                     <  6:151363/149‑1 (MQ=255)                                                                                                             tcgCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCt                                                   <  2:64441/132‑1 (MQ=255)                                                                                                             tcgCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCt                                                   >  1:64441/1‑132 (MQ=255)                                                                                                             tcgCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGc                                  >  4:257890/1‑149 (MQ=255)                                                                                                               gagtCCAGGTCTACGGCGAGCTGCGGGTCGACGTCGGCAGTGTGATTCTGGCCGACGACTACTGTCAGCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGGGGTTGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCg                                <  8:170929/146‑1 (MQ=255)                                                                                                                 ttCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCa                              >  5:337037/1‑149 (MQ=255)                                                                                                                   ccAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCAc                             <  7:375525/148‑1 (MQ=255)                                                                                                                       gTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGt                                                                                          >  3:292145/1‑83 (MQ=255)                                                                                                                       gTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGt                                                                                          <  4:292145/83‑1 (MQ=255)                                                                                                                         ataCGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTa                        >  3:32342/1‑147 (MQ=255)                                                                                                                            cGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGcc                   <  7:200555/149‑1 (MQ=255)                                                                                                                            cGGCGAGCTGCCGGTCGAAGGCGACAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGcc                   <  8:178046/149‑1 (MQ=255)                                                                                                                            cGGCGAGCTGCCGGTCGAAGGCGACAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGcc                   <  8:172632/149‑1 (MQ=255)                                                                                                                              gCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCg                  >  3:129530/1‑148 (MQ=255)                                                                                                                                  gCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGccca                                                             >  1:338297/1‑99 (MQ=255)                                                                                                                                  gCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCaa                                                             <  2:338297/101‑1 (MQ=255)                                                                                                                                  gCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATc                                            <  4:183946/118‑1 (MQ=255)                                                                                                                                  gCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATc                                            >  3:183946/1‑118 (MQ=255)                                                                                                                                     gCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTAc                       >  4:86840/1‑136 (MQ=255)                                                                                                                                     gCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTAc                       <  3:86840/136‑1 (MQ=255)                                                                                                                                        ggTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGCACCCGCGAc        <  5:201083/148‑1 (MQ=255)                                                                                                                                        ggTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGCACCCGCGAc        >  6:201083/1‑148 (MQ=255)                                                                                                                                             aaGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGc                 <  4:56663/134‑1 (MQ=255)                                                                                                                                             aaGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGc                 >  3:56663/1‑134 (MQ=255)                                                                                                                                             aaGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGCACCCGCGACTTGTac  >  8:56346/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                               AACCACAACGTGCTCAACGCCCTGGCGACCATCGCCATCGCCACTGACGAAGGCATTACTGACGAGGCCATCGTCCAGGGCCTGTCCGG‑CTTCCAGGGTGTCGGTCGTCGCTTCCAGGTATACGGCGAGCTGCCGGTCGAAGGCGGCAGTGTGATGCTGGTCGACGACTACGGTCATCACCCCACCGAAGTGGCGGCGGTGATCAAGGCTGTGCGCGGTGGCTGGCCAAGCCGCCGCCTGGTGATCGTCTACCAGCCGCACCGCTACAGCCGCACCCGCGACTTGTAC  >  NC_002947/1525230‑1525517

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