Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	663,973	T→C	5.9%	intergenic (‑10/+129)	dacA ← / ← rlpA	D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)/septal ring protein, suppressor of prc, minor lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	663,973	0	T	C	5.9%	208.3 / 3.1	68	intergenic (‑10/+129)	dacA/rlpA	D‑alanyl‑D‑alanine carboxypeptidase (penicillin‑binding protein 5)/septal ring protein, suppressor of prc, minor lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (32/32);  new base C (2/2);  total (34/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.89e-01

CGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCG  >  NC_000913/663826‑664111                                                                                                                                                    |                                                                                                                                           cGCTTCGATCTGCTGTACACGCGGGATCATAGTTTTTATTTTCAGGTCATCGGCATGTGCAGCAGAGATAATGGCTGTGGAAAGAGCCGTGGTGAGCGCGAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACt                                                                                                                                            <  2:86229/148‑1 (MQ=255)    aTCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACg                                                                                                                                         <  5:128058/148‑1 (MQ=255)    aTCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACg                                                                                                                                         <  8:106174/148‑1 (MQ=255)         tCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATc                                                                                                                                    >  5:288826/1‑148 (MQ=255)                   acCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGaa                                                                                                                          >  1:371904/1‑148 (MQ=255)                      cGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAAtt                                                                                                                        >  7:11217/1‑147 (MQ=255)                      cGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAAtt                                                                                                                        >  7:41316/1‑147 (MQ=255)                       gggATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTa                                                                                                                       >  4:229238/1‑147 (MQ=255)                         gATCATATTTTTGATATTGAGGTCATCGGCATGTTCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTaaaa                                                                                                                    <  6:288826/148‑1 (MQ=255)                         gATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAATAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTaaaa                                                                                                                    <  4:250391/148‑1 (MQ=255)                         gATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTaaaa                                                                                                                    <  5:4077/148‑1 (MQ=255)                          aTCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTaaaaa                                                                                                                   <  3:229238/148‑1 (MQ=255)                           tCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTaaaa                                                                                                                    >  1:358777/1‑146 (MQ=255)                               aGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGc                                                                                                              >  3:367737/1‑148 (MQ=255)                                   ttGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCct                                                                                                           <  5:338028/147‑1 (MQ=255)                                     gATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCctac                                                                                                         <  8:337999/147‑1 (MQ=255)                                      atatTCAGGGCATCGGCATGTGCAGCAGTTATAAAGGCTGGTCATAGAGCCGTGGTGAGCGCCAGGCGCTTCATTATCCGAGCGTAAACTATGGGATTCATCGTCTGAAGTACGACATCCGTGATGGAATTAAAAAAGTGCCctact                                                                                                        <  3:26106/147‑1 (MQ=255)                                            aGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCaa                                                                                                 <  5:305279/148‑1 (MQ=255)                                                       aTGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACtat                                                                                                      <  6:164659/132‑1 (MQ=255)                                                       aTGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACtat                                                                                                      >  5:164659/1‑132 (MQ=255)                                                         gtgCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCGCTACCgg                                                                                     >  8:337366/1‑147 (MQ=255)                                                         gtgCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCgg                                                                                     >  3:182767/1‑147 (MQ=255)                                                               cagAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACCACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCAt                                                                              >  5:112517/1‑148 (MQ=255)                                                                agagATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATc                                                                             >  6:165296/1‑148 (MQ=255)                                                                  agaTAAAGGCTGTGCAAAGAGCCGTTGTGAGCGCCAGGCGCTTCATGATACGAGTGGAAAAAATGGGATTCATGGTCTGAACTACGACATCCGTGGTGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTg                                                                           <  5:102383/148‑1 (MQ=255)                                                                   gaTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAg                                                                                                                 >  1:159221/1‑109 (MQ=255)                                                                   gaTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAg                                                                                                                 <  2:159221/109‑1 (MQ=255)                                                                    aTAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGAc                                                                         <  8:408317/148‑1 (MQ=255)                                                                     tAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCaa                                                                                                 >  8:381058/1‑123 (MQ=255)                                                                     tAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCaa                                                                                                 >  8:273519/1‑123 (MQ=255)                                                                     tAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCaa                                                                                                 <  7:381058/123‑1 (MQ=255)                                                                     tAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCaa                                                                                                 <  7:273519/123‑1 (MQ=255)                                                                     tAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACt                                                                        >  4:362135/1‑148 (MQ=255)                                                                      aaaGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACtt                                                                       >  8:50343/1‑148 (MQ=255)                                                                           cTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCcgc                                                                  >  8:364947/1‑148 (MQ=255)                                                                             gtgCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAATAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGGATGTGACTTTCCGCGt                                                                >  4:61528/1‑148 (MQ=255)                                                                               gCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGa                                                              >  2:269533/1‑148 (MQ=255)                                                                                 aaaGTGCCGTGGTGAGCGCCAGGCGATTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTg                                                                                                               <  1:34366/97‑1 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGaa                                                                                                                          <  5:283598/85‑1 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGaa                                                                                                                          >  6:283598/1‑85 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACtt                                                           <  4:367737/148‑1 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACtt                                                           <  8:11217/148‑1 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACtt                                                           <  2:172353/148‑1 (MQ=255)                                                                                  aaGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGGCTGAACTACGACATACGGGATGGAATTAAAAAAGTGc                                                                                                              >  2:34366/1‑97 (MQ=255)                                                                                     agCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGt                                                        <  4:182767/148‑1 (MQ=255)                                                                                          tggtgAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCAGCGTGACTTTGTTAACg                                                   >  6:389610/1‑148 (MQ=255)                                                                                           ggtgAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACCACGACATCCGTGATGGAATTAAAAAAGTGcc                                                                                                             <  8:195850/89‑1 (MQ=255)                                                                                           ggtgAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACCACGACATCCGTGATGGAATTAAAAAAGTGcc                                                                                                             >  7:195850/1‑89 (MQ=255)                                                                                               aGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACCACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCAt                                               <  6:112517/147‑1 (MQ=255)                                                                                                  gcCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTAc                                           <  5:165296/148‑1 (MQ=255)                                                                                                     aGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAg                                         >  6:90813/1‑147 (MQ=255)                                                                                                      ggCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCAc                                                                                           >  7:376797/1‑96 (MQ=255)                                                                                                      ggCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCAc                                                                                           <  8:376797/96‑1 (MQ=255)                                                                                                                 aTACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATc                            <  2:369070/148‑1 (MQ=255)                                                                                                                    cGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCGCTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGa                         <  7:337366/148‑1 (MQ=255)                                                                                                                    cGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCTGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          <  2:371904/147‑1 (MQ=255)                                                                                                                    cGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          <  1:269533/147‑1 (MQ=255)                                                                                                                    cGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          <  7:50343/147‑1 (MQ=255)                                                                                                                                aTGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTAAAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTAACAGAAATTGCCACATCAGATGCCTGCTTTAc             >  1:290070/1‑148 (MQ=255)                                                                                                                                     aTTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTAc        >  4:221926/1‑148 (MQ=255)                                                                                                                                       tCATGGTCTGAACTACGACATCCGTGATGGAATTTAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGCCTTTGTGAACGACATTTACAGAAATTGACACATCAGATGCCTGCTTTACGcaccgg      >  2:293568/1‑143 (MQ=255)                                                                                                                                       tCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          <  5:29816/128‑1 (MQ=255)                                                                                                                                       tCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          >  6:29816/1‑128 (MQ=255)                                                                                                                                       tCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTg      >  3:122718/1‑148 (MQ=255)                                                                                                                                         aTGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTgcg    >  5:296713/1‑148 (MQ=255)                                                                                                                                           ggTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGcgcg  >  5:148668/1‑148 (MQ=255)                                                                                                                                            gTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGcgcg  <  2:358777/147‑1 (MQ=255)                                                                                                                                                gAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          >  3:31471/1‑119 (MQ=255)                                                                                                                                                gAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAg                          <  4:31471/119‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                           CGCATCGATCTGCGGTACACCCGGGATCATAGTTTTGATATTCAGGTCATCGGCATGTGCAGCAGAGATAAAGGCTGTGCAAAGAGCCGTGGTGAGCGCCAGGCGCTTCATGATACGAGCGGAAAAAATGGTATTCATGGTCTGAACTACGACATCCGTGATGGAATTAAAAAAGTGCCCTACTATAGCAAATGCACTACCGGCAGGCATCTGACTTTCCGCGTGACTTTGTTAACGTCATTTACAGAAATTGACACATCAGATGCCTGCTTTACGCTACTGCGCG  >  NC_000913/663826‑664111

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