Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,622,479	A→G	6.3%	L420P (CTG→CCG)	PP_4941 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,622,479	0	A	G	6.3%	188.4 / 2.6	63	L420P (CTG→CCG)	PP_4941	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/31);  new base G (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.56e-01

CTTACAGCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCACTGTCGCCCATCGGCCAGCGGTCGT  >  NC_002947/5622333‑5622626                                                                                                                                                          |                                                                                                                                                    cTTACAGCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAgg                                                                                                                                                   <  2:215438/149‑1 (MQ=255) cTTACAGCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAg                                                                                                                                                    <  3:115417/148‑1 (MQ=255)     cAGCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGt                                                                                                                                               <  7:262151/149‑1 (MQ=255)       gCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGc                                                                                                                                             <  5:275766/149‑1 (MQ=255)           aCGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGcatc                                                                                                                                          <  7:211241/148‑1 (MQ=255)               cgccAGAATATTTCGTGTCGGCGCTCGGCCGTGCGGTCGAAATCGTTCTCCCCTTCTATGTAGCCCGTTGGTTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGcat                                                                                                                                           <  1:259980/149‑1 (MQ=255)                gccAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGcc                                                                                                                                     >  6:315446/1‑148 (MQ=255)                 ccAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGcc                                                                                                                                     >  4:87140/1‑147 (MQ=255)                  cAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAg                                                                                                                                  <  2:136743/149‑1 (MQ=255)                  cAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAg                                                                                                                                  >  2:27707/1‑149 (MQ=255)                    gAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑TCGATGGCTGTCGACCGGCAAGCCAGCGTTGCAGCAGGAGGCGCAGGCGC‑CGCTTCAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAggg                                                                                                                                >  6:220824/1‑149 (MQ=255)                            cGTGGCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGc                                                                                                                         <  7:105248/148‑1 (MQ=255)                             gtgtCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGcc                                                                                                                        <  4:29250/148‑1 (MQ=255)                             gtgtCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCt                                                                                                                       <  5:315446/149‑1 (MQ=255)                                tCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑TCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTTAACTGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTg                                                                                                                     <  5:220824/148‑1 (MQ=255)                                   gcgcACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGTGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTcgc                                                                                                                 >  6:234605/1‑149 (MQ=255)                                                    agaaTTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTcc                                                                                                <  3:87140/149‑1 (MQ=255)                                                        ttCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGc                                                                                                                                             >  8:152541/1‑100 (MQ=255)                                                        ttCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGc                                                                                                                                             <  7:152541/100‑1 (MQ=255)                                                           gTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGGTCGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCAc                                                                                           >  1:94608/1‑148 (MQ=255)                                                           gTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCAc                                                                                           >  1:264961/1‑147 (MQ=255)                                                            ttCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGc                                                                                                                               >  6:181368/1‑110 (MQ=255)                                                            ttCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGc                                                                                                                               <  5:181368/110‑1 (MQ=255)                                                                   c‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGTGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAg                                                                                 <  5:234605/149‑1 (MQ=255)                                                                        aT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAgg                                                                                <  4:139081/149‑1 (MQ=255)                                                                           aaG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGt                                                                              >  2:151999/1‑149 (MQ=255)                                                                            aG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTg                                                                             <  2:57153/149‑1 (MQ=255)                                                                             g‑gCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAgg                                                                                       >  3:196653/1‑138 (MQ=255)                                                                             g‑gCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAgg                                                                                       <  4:196653/138‑1 (MQ=255)                                                                                      cTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCCGCACCAGGCTGCAGGGTGCCCAAGCCt                                                                    >  8:141461/1‑149 (MQ=255)                                                                                      cTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAGGctt                                                                    <  1:27707/149‑3 (MQ=255)                                                                                       tGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTg                                                                   >  2:131995/1‑149 (MQ=255)                                                                                       tGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTg                                                                   >  6:86235/1‑149 (MQ=255)                                                                                           gCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTggg                                                                 >  1:486/1‑147 (MQ=255)                                                                                            cccGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCGGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTg                                                                                    <  4:285967/127‑1 (MQ=255)                                                                                            cccGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCGGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTg                                                                                    <  4:290562/127‑1 (MQ=255)                                                                                            cccGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCGGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTg                                                                                    >  3:290562/1‑127 (MQ=255)                                                                                            cccGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCGGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTg                                                                                    >  3:285967/1‑127 (MQ=255)                                                                                              cGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGAGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAgcg                                                            >  4:48824/1‑149 (MQ=255)                                                                                                 cAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCt                                                          >  7:117286/1‑148 (MQ=255)                                                                                                     ccAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAAtt                                                     >  2:307300/1‑149 (MQ=255)                                                                                                      cAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTg                                                    >  3:40291/1‑149 (MQ=255)                                                                                                      cAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTg                                                    >  8:97984/1‑149 (MQ=255)                                                                                                       aGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTgc                                                   <  1:151999/149‑1 (MQ=255)                                                                                                       aGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTgc                                                   <  2:61310/149‑1 (MQ=255)                                                                                                        gCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTACAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTgcg                                                  >  6:276478/1‑149 (MQ=255)                                                                                                           ttGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTgcggcg                                               <  3:48824/149‑1 (MQ=255)                                                                                                              cagcagGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGt                                             >  8:294938/1‑148 (MQ=255)                                                                                                               agcagGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGACCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTa                                            >  1:189126/1‑148 (MQ=255)                                                                                                                gcaggcggcgcaggcgg‑cGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTaaa                                          <  2:264961/149‑1 (MQ=255)                                                                                                                    gcggcgcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCa                                                               <  1:14018/124‑1 (MQ=255)                                                                                                                    gcggcgcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCa                                                               >  2:14018/1‑124 (MQ=255)                                                                                                                    gcggcgcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCa                                      >  6:313793/1‑149 (MQ=255)                                                                                                                     cggcgcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCAc                                     <  8:291734/149‑1 (MQ=255)                                                                                                                         gcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCa                                                                                  <  3:148647/100‑1 (MQ=255)                                                                                                                         gcaggcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCa                                                                                  >  4:148647/1‑100 (MQ=255)                                                                                                                             gcgg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGc                             <  8:27064/149‑1 (MQ=255)                                                                                                                               gg‑cgcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCa                           <  2:94608/149‑1 (MQ=255)                                                                                                                                   gcTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCACTg                        >  3:111842/1‑149 (MQ=255)                                                                                                                                        aaCGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCa                           <  7:191572/141‑1 (MQ=255)                                                                                                                                        aaCGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCa                           >  8:191572/1‑141 (MQ=255)                                                                                                                                            gTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCACTGTCGCCCATc               <  2:119175/149‑1 (MQ=255)                                                                                                                                                         cAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCACTGTCGCCCATCGGCCAGCGGTCGt  <  2:486/149‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                    CTTACAGCTCACGCCGCCAGAACAGTTCGTGTCGGCGCACGGCCTTGCGGAAGAATTCGTTCTCCCC‑‑‑‑AT‑AAG‑GCGATGGCTGTCGCCCGGCAAGCCAGCGTTGCAGCAGGCGGCGCAGGCGG‑CGCTTGAACGGTGCCGGCGGCTGCAGGTCGTGCATCATGCCCAGGGCCATCGCCTTGTCGCGCTGGATCGCGCAGTCCAGCACCAGGCTGCAGGGTGCCCAAGCCTGGGCAGCGCTCAATTGCGGCGGTAAAGCCACGCCATCGCCACTGTCGCCCATCGGCCAGCGGTCGT  >  NC_002947/5622333‑5622626

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