Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,784,032	C→G	9.4%	Q15H (CAG→CAC) ‡	glnD ←	uridylyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,784,032	0	C	G	9.4%	173.4 / 5.5	61	Q15H (CAG→CAC) ‡	glnD	uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (27/28);  new base G (3/3);  total (30/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTA  >  NC_002947/1783885‑1784171                                                                                                                                                       |                                                                                                                                            cAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCt                                                                                                                                            >  3:197929/1‑149 (MQ=255)   ggATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTgg                                                                                                                                          <  1:2108/149‑1 (MQ=255)      tatTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAAc                                                                                                                                       >  1:325613/1‑149 (MQ=255)          gTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGGGGGCCAGTTCCGCCTGGAACTGGc                                                                                                                                   >  7:203867/1‑149 (MQ=255)           tccactaGCCAGGCGCGGGCCTCGATCAGCGGTCGTATGTCGCTGCCGGCGCTAAACCGCTTGACCAGCA‑CCC‑C‑CCCGGTAAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTGGc                                                                                                                                   <  1:265448/142‑1 (MQ=255)                 aGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGGGGGCCAGTTCCGCCTGGAACTGGccctcaag                                                                                                                            <  8:203867/149‑7 (MQ=255)                    cAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                    <  7:141632/138‑1 (MQ=255)                    cAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCAGGCGGAACTgg                                                                                                                                    >  8:141632/1‑138 (MQ=255)                       gcgcGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAg                                                                                                                      <  2:301135/149‑1 (MQ=255)                              ccTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCgggg                                                                                                               <  3:189401/149‑1 (MQ=255)                                 cGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTc                                                                                                             >  2:313824/1‑148 (MQ=255)                                    tgaGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCACCCTACAGTCGGCCAGGCGGATGGCTATCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTcc                                                                                                            <  1:105463/147‑1 (MQ=255)                                      aGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTg                                                                                                       <  4:197929/149‑1 (MQ=255)                                               aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTgg                                                                                                                                          >  7:92464/1‑105 (MQ=255)                                               aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTgg                                                                                                                                          <  8:92464/105‑1 (MQ=255)                                               aTGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATcgc                                                                                              >  7:300344/1‑149 (MQ=255)                                                     cTGCCGGCGCTAAAAGGTTTGTCGAGCA‑CCT‑C‑GCCGGCCTGGCGTATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGTGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACg                                                                                        <  2:325613/149‑1 (MQ=255)                                                        ccGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTc                                                                                     >  5:211948/1‑149 (MQ=255)                                                        ccGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTc                                                                                     <  3:30535/149‑1 (MQ=255)                                                         cGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTcc                                                                                    >  4:302426/1‑149 (MQ=255)                                                               gAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                              >  1:55657/1‑149 (MQ=255)                                                               gAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                              >  4:291008/1‑149 (MQ=255)                                                               gAAAACGTTTGGCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGCTGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGt                                                                              <  7:202575/149‑1 (MQ=255)                                                                aaaaCGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTg                                                                             >  6:126514/1‑149 (MQ=255)                                                                  aaCGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTACTTGAAGGCGGCGAGGGGGCTCGCCTTCAGGGCCAGTTCCGCGTGGAACTGGCCACGGACGAACAGCTCGGGGACAACCTGGGGCATCGCGTCACGTTCGTGTCGTggg                                                                           >  2:96546/1‑149 (MQ=255)                                                                   aCGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTgggg                                                                          <  7:20689/149‑1 (MQ=255)                                                                      tttGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTgg                                                                                                                                    >  2:105007/1‑88 (MQ=255)                                                                        tGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGc                                                                      <  5:126514/148‑1 (MQ=255)                                                                            gAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTc                                                                 >  2:49685/1‑149 (MQ=255)                                                                               cA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCggg                                                              <  3:302426/149‑1 (MQ=255)                                                                                   cT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                              <  4:143512/98‑1 (MQ=255)                                                                                   cT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGt                                                                                                              >  3:143512/1‑98 (MQ=255)                                                                                    t‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGc                                                          >  8:60071/1‑149 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCgg                                                        <  1:96546/149‑1 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCg                                                         <  2:55657/148‑1 (MQ=255)                                                                                        gCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCg                                                         <  8:300344/148‑1 (MQ=255)                                                                                           ggCCAGGCGTATGGCTTTCTTGAAGGCGGCGATGGGGCTGGCGTTGAGGGCCAGTTCCGCGTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCTTCGCGTCTCGTTCCTGTCTTGGGGTGGCGTGTCGGGGTGCGGTCa                                                     <  8:28518/149‑1 (MQ=255)                                                                                             ccAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAgg                                                   >  1:147462/1‑149 (MQ=255)                                                                                                         cTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTgcgcc                                       >  5:159033/1‑148 (MQ=255)                                                                                                          tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGcgcg                                       >  5:111583/1‑148 (MQ=255)                                                                                                          tttCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCgg                                      >  6:94229/1‑149 (MQ=255)                                                                                                            tCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGcgcg                                       <  6:111583/146‑1 (MQ=255)                                                                                                                 aaGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgtg                               <  7:60071/149‑1 (MQ=255)                                                                                                                 aaGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGgt                                <  1:49685/148‑1 (MQ=255)                                                                                                                     cggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgt                           >  6:148923/1‑149 (MQ=255)                                                                                                                      ggcgATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTgtcgtc                          <  1:313824/149‑1 (MQ=255)                                                                                                                                 tCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAggg               >  8:3852/1‑149 (MQ=255)                                                                                                                                     cTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTg             >  8:202118/1‑147 (MQ=255)                                                                                                                                     cTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGa            <  5:219505/148‑1 (MQ=255)                                                                                                                                     cTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGa            >  6:219505/1‑148 (MQ=255)                                                                                                                                         gtggCCAGTTCCAGGCGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCGGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGGGGGATGGTGTCGTGGTTGCGCCGGGTgatgatg       >  4:54493/3‑148 (MQ=255)                                                                                                                                         aGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtc       >  1:95552/1‑149 (MQ=255)                                                                                                                                          gggCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtct      <  7:3852/149‑1 (MQ=255)                                                                                                                                            gCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc     <  5:31654/148‑1 (MQ=255)                                                                                                                                            gCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc     >  6:31654/1‑148 (MQ=255)                                                                                                                                             ccAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc     >  7:3934/1‑147 (MQ=255)                                                                                                                                             ccAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGAtctc     <  8:3934/147‑1 (MQ=255)                                                                                                                                              cAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTa  <  5:148923/149‑1 (MQ=255)                                                                                                                                              cAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTa  <  4:94533/149‑1 (MQ=255)                                                                                                                                                     gCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTTAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTg                      <  1:39974/122‑1 (MQ=255)                                                                                                                                                     gCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTg                      >  2:39974/1‑122 (MQ=255)                                                                                                                                                       |                                                                                                                                            CAGGATATTGTCGACGAGCCAGGCGCGGGCCTCGATCAGCGGTCGGATGTCGCTGCCGGCGCGAAAACGTTTGTCGAGCA‑CCT‑C‑GCCGGCCAGGCGGATGGCTTTCTTGAAGGCGGCGATGGGGCTCGCCTTGAGGGCCAGTTCCGCCTGGAACTGGCCACGGTCGAACAGCTCGGGGTCCACCTGGGGCATCGCGTCACGTTCCTGTCGTGGGGTGGCCTGTCGGGGTGCGGTCAGGCCGAAGTGCGCGGGATGGTGTCGTCCTTGCGCAGGGTGAAGATCTCGTA  >  NC_002947/1783885‑1784171

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