Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,534,893	C→T	6.3%	intergenic (+84/‑99)	secA → / → argJ	protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,534,893	0	C	T	6.3%	286.3 / 5.5	95	intergenic (+84/‑99)	secA/argJ	protein translocase subunit SecA/glutamate N‑acetyltransferase/amino acid acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (42/47);  new base T (4/2);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

GCAACGAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCC  >  NC_002947/1534745‑1535029                                                                                                                                                     |                                                                                                                                           gCAACGAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCCCAGGTACCGCGCTCGTc                                                                                                                                            >  4:109868/1‑149 (MQ=255)   aaCGAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCTGACCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTcc                                                                                                                                           >  6:216831/1‑148 (MQ=255)   aaCGAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCa                                                                                                                                          >  5:5490/1‑149 (MQ=255)       aGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGa‑‑a                                                                                                                                    >  4:94147/1‑149 (MQ=255)         ccGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACg                                                                                                                                  >  6:64071/1‑149 (MQ=255)          cGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACgg                                                                                                                                 <  5:373306/149‑1 (MQ=255)           gTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACggg                                                                                                                                <  5:64071/149‑1 (MQ=255)                      gTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtg                                                                                                                     >  3:118176/1‑149 (MQ=255)                          gggCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTgg                                                                                                                  <  3:93942/148‑1 (MQ=255)                           ggCAAGAAGTTCAAGCATTGCCTCGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGGCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTgggg                                                                                                                <  8:355661/149‑1 (MQ=255)                            gCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTggggg                                                                                                               >  7:13133/1‑149 (MQ=255)                              aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcg                                                                                                             >  4:249753/1‑149 (MQ=255)                              aagaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcg                                                                                                             >  5:56288/1‑149 (MQ=255)                               agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcgg                                                                                                            >  5:345770/1‑149 (MQ=255)                               agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGggcgg                                                                                                            >  5:284679/1‑149 (MQ=255)                               agaagTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagcgg                                                                                                            <  4:285183/149‑1 (MQ=21)                                  acttCAAGCATTGCCACGGTCAGATCGAGTGATCTTAACCTTACGTTGTACAAGCCTGGGGTCTCTTTGCGGCCCTTTCGCAACACACGGCCGCTCCCACAGGGACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCt                                                                                                         <  5:216831/147‑1 (MQ=255)                                   gTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCtt                                                                                                        <  8:211371/149‑1 (MQ=255)                                    ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtg                                                                                                       >  2:308843/1‑149 (MQ=255)                                    ttCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtg                                                                                                       <  8:13133/149‑1 (MQ=255)                                      cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtgtg                                                                                                     >  5:287629/1‑149 (MQ=255)                                      cAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTtgt                                                                                                      <  6:345770/148‑1 (MQ=255)                                       aaGCATTGCCACGGTCGGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTgtgt                                                                                                    <  4:261588/149‑1 (MQ=255)                                        aGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTc                                                                                                   >  7:30468/1‑149 (MQ=255)                                         gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                  <  6:287629/149‑1 (MQ=255)                                         gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                  >  5:263703/1‑149 (MQ=255)                                         gCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcg                                                                                                  >  1:57928/1‑149 (MQ=255)                                          cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                 >  4:270980/1‑149 (MQ=255)                                          cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                 <  4:276394/149‑1 (MQ=255)                                          cATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTcgc                                                                                                 >  4:275894/1‑149 (MQ=255)                                              gCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATg                                                                                             >  5:28961/1‑149 (MQ=255)                                               ccACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgg                                                                                            <  6:56288/149‑1 (MQ=255)                                                cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                           >  4:7119/1‑149 (MQ=255)                                                cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                           <  3:7119/149‑1 (MQ=255)                                                cACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                           <  3:322427/149‑1 (MQ=255)                                                 aCGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATgggg                                                                                          >  1:237045/1‑149 (MQ=255)                                                   ggTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTg                                                                                        >  4:87562/1‑149 (MQ=255)                                                         aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGc                                                                                  <  1:2591/149‑1 (MQ=255)                                                         aTCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCAATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGc                                                                                  >  5:152057/1‑149 (MQ=255)                                                            gAGTGATCTGAACCTTACGTTGTCCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCCACACCAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGC‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCgg                                                                                <  8:260623/148‑1 (MQ=255)                                                             aGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagc                                                                                                              >  5:94782/1‑117 (MQ=14)                                                             aGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGagc                                                                                                              <  6:94782/117‑1 (MQ=14)                                                                     gAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCCCGATGTGGTGCGAAGCGGCCCCTTTTTc                                                                      >  5:219089/1‑149 (MQ=255)                                                                      aaCCTTACGTTGTACATGCCTGGGGTCGCTTTGCGGCCCTTTCGCATCACACGGCCGGTCCCACAGGTAGCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCttttt                                                                       <  3:94147/147‑1 (MQ=255)                                                                      aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                      <  7:295360/148‑1 (MQ=255)                                                                      aaCCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTc                                                                      <  3:368668/148‑1 (MQ=255)                                                                            aCGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAcatc                                                                <  6:28961/148‑1 (MQ=255)                                                                                   aCAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttt                                                        >  4:250641/1‑149 (MQ=255)                                                                                     aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                      <  3:275894/149‑1 (MQ=255)                                                                                     aaGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGtttttt                                                      <  3:270980/149‑1 (MQ=255)                                                                                       gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTc                                                     >  4:6543/1‑148 (MQ=255)                                                                                       gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTc                                                     >  6:166860/1‑148 (MQ=255)                                                                                       gCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCt                                                    <  2:237045/149‑1 (MQ=255)                                                                                         cTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAg                                                  <  3:249753/149‑1 (MQ=255)                                                                                          tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtgtg                                                                                                                   <  8:250951/83‑1 (MQ=255)                                                                                          tGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTtgtg                                                                                                                   >  7:250951/1‑83 (MQ=255)                                                                                             ggTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAgc                                              >  1:5645/1‑149 (MQ=255)                                                                                              gTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCATTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAgcg                                             <  5:143941/149‑1 (MQ=255)                                                                                                 gCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGctc                                          <  6:5490/149‑1 (MQ=255)                                                                                                      gCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAAt                                      <  3:250641/148‑1 (MQ=255)                                                                                                             tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATTA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                <  3:160090/147‑1 (MQ=255)                                                                                                             tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                              >  1:155973/1‑149 (MQ=255)                                                                                                             tttCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGt                                >  4:160090/1‑147 (MQ=255)                                                                                                               tCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                              <  5:166860/147‑1 (MQ=255)                                                                                                               tCGCAACACAAGGCCGCTCCCACACGTACCGCGCCCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTATCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttg                              <  4:343232/147‑1 (MQ=255)                                                                                                                 gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                           <  5:185496/148‑1 (MQ=255)                                                                                                                 gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                           <  3:87562/148‑1 (MQ=255)                                                                                                                 gCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                           <  2:5645/148‑1 (MQ=255)                                                                                                                  aaaCAAACGGCCGCTCCCACCGGTACCGCGCTCGTCATTTACCACGCCATCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtc                           <  1:316238/148‑1 (MQ=255)                                                                                                                   aaCACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                        >  6:301811/1‑149 (MQ=255)                                                                                                                     cacaAGGCCGCTCCAACAGGTTCCGCGCTCGTCGATGT‑‑ACGGGCGCGATCTTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttg                        <  6:219089/147‑1 (MQ=255)                                                                                                                      acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                           >  8:310024/1‑79 (MQ=255)                                                                                                                      acaAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATggg                                                                                           <  7:310024/79‑1 (MQ=255)                                                                                                                       caAGGCCGCCCCCACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTTGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTcc                    >  3:153127/1‑149 (MQ=18)                                                                                                                        aaGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                   >  5:363619/1‑149 (MQ=255)                                                                                                                         aGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCtt                  <  3:6543/149‑1 (MQ=255)                                                                                                                          ggCCGCTCCCACAGGTACCGCGCTCGGCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCttt                 <  8:296424/149‑1 (MQ=255)                                                                                                                           gCCGCTCCCACAGGTACCGCGCGCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTg                <  8:30468/149‑1 (MQ=255)                                                                                                                            ccGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGc               >  2:369941/1‑149 (MQ=255)                                                                                                                               cccccACAAGGATCGCGCCCGTTCATGG‑‑ACGAGCGCGGTACCTGTGGGAGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCa            >  1:355360/3‑149 (MQ=21)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAc                                                                   >  6:42453/1‑91 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGAc                                                                   <  5:42453/91‑1 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAccc         >  3:162065/1‑149 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAcc          <  6:263703/148‑1 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCAc           <  6:236407/147‑1 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCC‑TTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttggtc                     <  5:51857/136‑1 (MQ=255)                                                                                                                                  ccACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCC‑TTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtcttggtc                     >  6:51857/1‑136 (MQ=255)                                                                                                                                    acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                          <  7:258767/130‑1 (MQ=255)                                                                                                                                    acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGttggtct                          >  8:258767/1‑130 (MQ=255)                                                                                                                                    acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg       >  2:225290/1‑149 (MQ=255)                                                                                                                                    acaGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCC‑TTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTg       >  4:198298/1‑148 (MQ=255)                                                                                                                                      aGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa     >  5:200450/1‑149 (MQ=255)                                                                                                                                      aGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCa     >  1:66487/1‑149 (MQ=255)                                                                                                                                        gTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                   >  7:64071/1‑133 (MQ=255)                                                                                                                                        gTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCt                   <  8:64071/133‑1 (MQ=255)                                                                                                                                         gatcgcgCCCGTTCATGG‑‑ACGAGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc  >  7:243391/4‑149 (MQ=39)                                                                                                                                          aCCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCTATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCAccc  <  4:118176/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           GCAACGAGCCGTGCTGGTGTGGTTCGGGCAAGAAGTTCAAGCATTGCCACGGTCAGATCGAGTGATCTGAACCTTACGTTGTACAAGCCTGGGGTCGCTTTGCGGCCCTTTCGCAACACAAGGCCGCTCCCACAGGTACCGCGCTCGTCCATGA‑‑ACGGGCGCGATCCTTGTGGGGGCGGCCTTGTGTCGCGATGGGGTGCGAAGCGGCCCCTTTTTCGACATCATCGTTTTTTCTAGGAGCGCTCTAATGGCTGTTGGTCTTGGTCCTTTGCCCACCCTGCACCC  >  NC_002947/1534745‑1535029

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