Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,808,548	G→T	9.2%	noncoding (412/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,548	0	G	T	9.2%	156.1 / 4.7	62	noncoding (412/434 nt)	PP_mr35	cspA regulatory region
Reads supporting (aligned to +/- strand):  ref base G (26/30);  new base T (2/4);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 6.81e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTACG  >  NC_002947/2808401‑2808695                                                                                                                                                    |                                                                                                                                                     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                                                                                                                                          <  8:267848/92‑1 (MQ=255)                                                                     gCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGtttttt                                                                                                                                          >  7:267848/1‑92 (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                                                             <  1:206046/149‑1 (MQ=255)                                                                                        gTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGaaa                                                             <  3:253298/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                                                      >  5:141442/1‑149 (MQ=255)                                                                                               aGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTg                                                      >  3:174028/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                                                                                                    <  2:51128/87‑1 (MQ=255)                                                                                                               tCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAg                                                                                                    >  1:51128/1‑87 (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)                                                                                                                                                    |                                                                                                                                                     AAAATGTGAAAAGCCTGCGCGGTCGCCTTGTAGGAGCGGGTTCACCCGCGAACACGGGCAAAGCCCGTGCCATTCACCCTACAGTCGGTTTCGCAGGCAATCCCGCCCCTTCAGATGATGCATCACCCATAAAAAAACCCGCCACAAGGGCGGGTTTTTTTACAGCTATCCGAATTACTCCGG‑ACGAACCTGTGCAGCCTGCATGCCCTTCTGGCCGCGCTCGGCGACGAAGGAAACAGTCTGGCCTTCTTTCAGGCTCTTGAAGCCGTCAGATTCGATGGCTTTGAAGTGTACG  >  NC_002947/2808401‑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: ‑