Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	994,947	T→A	11.2%	intergenic (+26/+99)	PP_0859 → / ← PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	994,947	0	T	A	11.2%	165.4 / 8.2	62	intergenic (+26/+99)	PP_0859/PP_0860	ketoglutaramate omega‑amidase/sulfite reductase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (28/27);  new base A (4/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCT  >  NC_002947/994801‑995089                                                                                                                                                   |                                                                                                                                                 ccTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCAGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTg                                                                                                                                                >  6:111807/1‑149 (MQ=255)    gCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGccc                                                                                                                                             >  4:283526/1‑149 (MQ=255)          tGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                               <  5:15666/141‑1 (MQ=255)          tGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGc                                                                                                                                               >  6:15666/1‑141 (MQ=255)            ccGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCAAGCCCGCGaaga                                                                                                                                      <  3:283526/148‑3 (MQ=255)                   ggCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGggcc‑‑ggcc                                                                                                                            >  6:65589/1‑149 (MQ=255)                     cGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCt                                                                                                                          <  6:193787/149‑1 (MQ=255)                              tgtTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACGGCAGGGCCGGCCCTTTCGCGGGCAAGCCCGCGAAGAAGCC‑‑GGGCCTgctgtcagt                                                                                                                 <  4:127087/149‑8 (MQ=21)                                   aCTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTc                                                                                                             <  8:307686/148‑1 (MQ=255)                                       cAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCg                                                                                                        >  5:25252/1‑149 (MQ=255)                                            cTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGaaa                                                                                                                                       <  8:81269/115‑1 (MQ=255)                                             tGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCg                                                                                                        >  6:95546/1‑143 (MQ=255)                                             tGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCg                                                                                                        <  5:95546/143‑1 (MQ=255)                                                     gcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTg                                                                                          >  6:264927/1‑149 (MQ=255)                                                     gcggAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTg                                                                                          >  4:267992/1‑149 (MQ=255)                                                            tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGcc                                                                                                     <  1:64100/131‑1 (MQ=255)                                                            tGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGcc                                                                                                     >  2:64100/1‑131 (MQ=255)                                                              gcggcgTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCCAAGGGGAGGCCAGCCAATTTCCACCGGCGGTTAACCGTTgcggccgc                                                                                   >  1:116226/1‑149 (MQ=255)                                                                         ggCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgcg                                                                      >  4:174872/1‑149 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                     <  5:65589/149‑1 (MQ=255)                                                                          gCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGCCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCgg                                                                     <  1:110111/149‑1 (MQ=255)                                                                             aaGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGt                                                                   >  7:269958/1‑148 (MQ=255)                                                                               gTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTaaa                                                                >  6:231505/1‑149 (MQ=255)                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGt                                                                   >  3:306595/1‑142 (MQ=255)                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGt                                                                   <  4:306595/142‑1 (MQ=255)                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCg                                                            >  7:168747/1‑149 (MQ=255)                                                                                   ccGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCg                                                            <  5:231505/149‑1 (MQ=255)                                                                                    cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCAc                                                                                                          <  3:78334/102‑1 (MQ=255)                                                                                    cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCAc                                                                                                          >  4:78334/1‑102 (MQ=255)                                                                                    cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCg                                                                              <  3:285958/130‑1 (MQ=255)                                                                                    cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCg                                                                              >  4:285958/1‑130 (MQ=255)                                                                                    cGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGc                                                           <  3:267992/149‑1 (MQ=255)                                                                                      gCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAACGGCC‑‑GGCCCTGCCAGCCAATTTCAACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCt                                                          >  3:207922/1‑148 (MQ=255)                                                                                         aaTCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTa                                                      <  8:168747/149‑1 (MQ=255)                                                                                                   ccGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATc                                            <  3:174872/149‑1 (MQ=255)                                                                                                    cGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCa                                           <  5:111807/149‑1 (MQ=255)                                                                                                      aTACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAAt                                         <  6:25252/149‑1 (MQ=255)                                                                                                          cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATacac                                     <  5:264927/149‑1 (MQ=255)                                                                                                          cTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATaca                                      >  7:67727/1‑148 (MQ=255)                                                                                                           ttCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCGCCTACACACCCATCAATacaa                                     >  5:2668/1‑147 (MQ=255)                                                                                                                 cTGCACTGACAGCAGGCCCGGCTTCTTCACGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                              >  2:169452/1‑85 (MQ=255)                                                                                                                 cTGCACTGACAGCAGGCCCGGCTTCTTCACGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCAcc                                                                                              <  1:169452/85‑1 (MQ=255)                                                                                                                    cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAtctc                                >  2:189012/1‑144 (MQ=255)                                                                                                                    cACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACAtctc                                <  1:189012/144‑1 (MQ=255)                                                                                                                      cTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTcgc                                                                       <  8:139960/103‑1 (MQ=255)                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGa                           <  8:248264/146‑1 (MQ=255)                                                                                                                       tGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGa                           >  7:248264/1‑146 (MQ=255)                                                                                                                            gcagGGCCGGCCCTTTCGCGGGTAAGCCCGCGAAGGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGcc                    >  1:260233/1‑148 (MQ=21)                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGTCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCg                      >  1:309601/1‑142 (MQ=255)                                                                                                                                gCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGTCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCg                      <  2:309601/142‑1 (MQ=255)                                                                                                                                  ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGTCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc             >  7:74504/1‑149 (MQ=255)                                                                                                                                  ccGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTc             >  5:362837/1‑149 (MQ=255)                                                                                                                                   cGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATCTCGcaccacag              >  4:24042/1‑140 (MQ=255)                                                                                                                                   cGGCTTCTTCGCGGGCTTGCACGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCTGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCt              <  3:24042/147‑1 (MQ=255)                                                                                                                                    ggccctttcGCGGGCAAGCCCGCGAAGAAGCC‑‑GGGCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGa           >  7:297590/7‑149 (MQ=21)                                                                                                                                      cttcttCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                     >  7:58663/1‑137 (MQ=255)                                                                                                                                      cttcttCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGc                     <  8:58663/137‑1 (MQ=255)                                                                                                                                           tcGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATCTCGCCCATCCTCGATCAActg    >  3:199334/1‑149 (MQ=255)                                                                                                                                           tcGCGGGCAAGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                             >  2:120585/1‑108 (MQ=255)                                                                                                                                           tcGCGGGCAAGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGAt                                             <  1:120585/108‑1 (MQ=255)                                                                                                                                            cgcgcgCTTGCCCGCGATAGGGCC‑‑GGCCCTTCCAGCCAATTTCCATCGTCAGTTGAACGTTGTCGCCGCCACCGTATTCGCGGGTAATCCCGATGCTAGACACCGATCAATAGACATCTCGCAGATATCGCCCATCCTCGATCAActt    <  4:207922/148‑2 (MQ=25)                                                                                                                                             gcgGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAActgct  <  8:67727/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 CCTGCTCAGTGCCGGGGAGGCGGACGGGGTGTTCACTGCAGTACTGAGCGCGGCGGAACTGGCGGCGTACCGGGCCAAGTTCCCGGCCAATCTGGATGCCGATACCTTCGAGCTGCACTGACAGCAGGCCCGGCTTCTTCGCGGGCTTGCCCGCGAAAGGGCC‑‑GGCCCTGCCAGCCAATTTCCACCGGCCGTTCACCGTTGCCGCCGCCACCGTATTCGCGGGTAAACCCGCTCCTACACACCGATCAATACACATCTCGCCGATATCGCCCATCCTCGATCAACTGCT  >  NC_002947/994801‑995089

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