Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,808,547	A→G	8.3%	noncoding (413/434 nt)	PP_mr35 ←	cspA regulatory region

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,808,547	0	A	G	8.3%	166.3 / 3.9	60	noncoding (413/434 nt)	PP_mr35	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base A (26/29);  new base G (2/3);  total (28/32)
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.53e-01

GCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTACG  >  NC_002947/2808399‑2808695                                                                                                                                                     |                                                                                                                                                      gcAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACaa                                                                                                                                                       <  6:279079/149‑1 (MQ=255)   aaaaTGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAgg                                                                                                                                                     >  3:142416/1‑149 (MQ=255)          gAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtt                                                                                                                                              >  8:67943/1‑149 (MQ=255)           aaaaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGACCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGatt                                                                                                                                             >  7:84227/1‑149 (MQ=255)             aaGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGttttt                                                                                                                                           <  7:67943/149‑1 (MQ=255)                   gcgcGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGTTTTTTTatgg                                                                                                                                     <  3:343570/149‑4 (MQ=39)                    cgcgGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCTTGTGGCGGGTTTTTTTatgg                                                                                                                                     <  6:139416/148‑4 (MQ=39)                           gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTTGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                                                                                                              <  6:364666/148‑1 (MQ=255)                           gCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCg                                                                                                                              <  5:151786/148‑1 (MQ=255)                             cTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCggg                                                                                                                                                >  4:209712/1‑128 (MQ=255)                             cTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCggg                                                                                                                                                <  3:209712/128‑1 (MQ=255)                                         ggTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGa                                                                                                              <  3:294157/149‑1 (MQ=255)                                           ttCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCCCAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAAc                                                                                                            >  3:233525/1‑149 (MQ=255)                                                 cgcgAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCtgtg                                                                                                       >  4:319613/1‑148 (MQ=255)                                                       cacGGGCTTTGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTg                                                                                                >  5:199601/1‑149 (MQ=255)                                                         cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCtgt                                                                                                        <  1:228165/139‑1 (MQ=255)                                                         cGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCtgt                                                                                                        >  2:228165/1‑139 (MQ=255)                                                            gCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGc                                                                                           <  6:92944/149‑1 (MQ=255)                                                                  cccGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCt                                                                                     <  3:319613/149‑1 (MQ=255)                                                                       gCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                          >  7:267848/1‑92 (MQ=255)                                                                       gCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                          <  8:267848/92‑1 (MQ=255)                                                                        ccATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGGCACGAACCTGTGCAGCCTGCATGCCCTTCTGGcc                                                                                 <  1:224292/148‑1 (MQ=255)                                                                              acacTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTAATCCGG‑ACGAACCTTTGCAGCCTGAATGCCCTTCTGGCCGCGCTCgg                                                                         >  8:214174/4‑149 (MQ=255)                                                                               cccTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGc                                                                        <  1:11742/149‑1 (MQ=255)                                                                                    cAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacg                                                                    >  4:158561/1‑148 (MQ=255)                                                                                     aGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGcgacg                                                                    >  7:118471/1‑147 (MQ=255)                                                                                       tCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAg                                                                 >  5:49111/1‑148 (MQ=255)                                                                                         ggTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaa                                                              <  6:199601/149‑1 (MQ=255)                                                                                          gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                             <  3:253298/149‑1 (MQ=255)                                                                                          gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                             <  1:206046/149‑1 (MQ=255)                                                                                          gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaa                                                              >  2:336586/1‑148 (MQ=255)                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAATGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCa                                                                                                     <  7:150834/104‑1 (MQ=255)                                                                                               gCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCa                                                                                                     >  8:150834/1‑104 (MQ=255)                                                                                                 aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                      >  3:174028/1‑149 (MQ=255)                                                                                                 aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                      >  5:141442/1‑149 (MQ=255)                                                                                                  ggCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑AAGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTgg                                                     >  1:265499/1‑149 (MQ=255)                                                                                                   gCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCt                                                                            <  6:307181/125‑1 (MQ=255)                                                                                                   gCAATCCCGCCACTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCt                                                                            >  5:307181/1‑125 (MQ=255)                                                                                                      aTCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCctt                                                 <  4:233525/149‑1 (MQ=255)                                                                                                            ccccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                           <  8:118471/149‑1 (MQ=255)                                                                                                              ccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAgg                                         >  8:21015/1‑149 (MQ=255)                                                                                                              ccTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAgg                                         <  4:174028/149‑1 (MQ=255)                                                                                                              ccTTCAGATGATGCATCACCCAT‑AAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGc                                        >  5:38910/1‑149 (MQ=255)                                                                                                               cTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCctt                                                 >  6:130713/1‑140 (MQ=255)                                                                                                               cTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCctt                                                 <  5:130713/140‑1 (MQ=255)                                                                                                                 tCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAg                                                                                                    >  1:51128/1‑87 (MQ=255)                                                                                                                 tCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAg                                                                                                    <  2:51128/87‑1 (MQ=255)                                                                                                                   aGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCtt                                    <  6:141442/149‑1 (MQ=255)                                                                                                                     atgatgCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGa                                  <  7:214174/149‑1 (MQ=255)                                                                                                                           cATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                            <  3:158561/149‑1 (MQ=255)                                                                                                                               aCCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGt                            >  3:23193/1‑145 (MQ=255)                                                                                                                                 ccATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGAtt                      >  5:130415/1‑149 (MQ=255)                                                                                                                                     aaaaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                           <  3:74660/124‑1 (MQ=39)                                                                                                                                     aaaaaaaCCCGCCCTTGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCa                                           >  4:74660/1‑124 (MQ=39)                                                                                                                                      aaaaaaCCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGAt                  <  7:21015/148‑1 (MQ=255)                                                                                                                                      aaaaaaACCGCCAGAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑TCTAACCTGTGCAGCCTTCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGATCTTGAAGCCGTCAGATTGGAt                  <  8:84227/148‑1 (MQ=255)                                                                                                                                       aaaaaCCCGCCCTGGTGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg                >  5:212506/1‑149 (MQ=39)                                                                                                                                            cccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGACCCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCGCTTGAAGCCGTCAGATTCGATgg                <  8:318811/144‑1 (MQ=255)                                                                                                                                            cccGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATgg                >  7:318811/1‑144 (MQ=255)                                                                                                                                              cGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGaa         <  1:336586/149‑1 (MQ=255)                                                                                                                                                  acaAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAgtgt     <  6:38910/149‑1 (MQ=255)                                                                                                                                                    aaGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATg                 <  2:138829/135‑1 (MQ=255)                                                                                                                                                    aaGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATg                 >  1:138829/1‑135 (MQ=255)                                                                                                                                                     aGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACTAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTACg  <  2:265499/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GCAAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTACG  >  NC_002947/2808399‑2808695

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