Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_bme_pgi	1,927,949	A→C	R297R (CGT→CGG)	ptrB ←	protease II

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_bme_pgi	1,927,949	0	A	C	98.3%	178.6 / NA	60	R297R (CGT→CGG)	ptrB	protease II
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base C (32/27);  minor base T (1/0);  total (33/27)
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.86e-01

TTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAG  >  NC_000913_3_bme_pgi/1927803‑1928089                                                                                                                                                    |                                                                                                                                             ttGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCg                                                                                                                                            <  1:153356/149‑1 (MQ=255)     ggTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGtat                                                                                                                                        >  1:128703/1‑149 (MQ=255)      gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGtata                                                                                                                                       >  2:1145930/1‑149 (MQ=255)      gTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGtata                                                                                                                                       >  2:1138868/1‑149 (MQ=255)       ttAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCGCGCATAC‑GGGTCCGGTATaa                                                                                                                                      <  2:1061374/149‑1 (MQ=255)                cgcAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACCGTTGCTCATCAAGCATAC‑GGGTCCGGTATAAGACAAAGtt                                                                                                                             >  1:234394/1‑149 (MQ=255)                  cAAACTGGTTAACCCGCGCTGACGCTCTGCAACCCCCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAATTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGtttt                                                                                                                           <  2:340601/149‑1 (MQ=255)                   aaaCTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGttttt                                                                                                                          <  2:1015018/149‑1 (MQ=255)                               ccGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCCGCCTGCTGCACCCACTTGGCGCATTTTACGCTGCCCACTCTTGATCACCACGCACAC‑GGGGCCGGTTTAAGCCAAAGTTTTTGCCGTGGCggct                                                                                                              >  2:878438/1‑147 (MQ=255)                                    cTGACGCTCTTCAACCACCAGCCAGTCGGTAAACGGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCCCTGTTGCTCATCACGCATAC‑GGGTCCGGAATAACCCAAAGTCTTTGCCCTGGCGGGTGGAAc                                                                                                         >  2:370432/1‑149 (MQ=255)                                    cTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAAc                                                                                                         <  1:1079789/149‑1 (MQ=255)                                    cTGACGCTACTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAAc                                                                                                         <  1:246780/149‑1 (MQ=255)                                         gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                                                                    <  1:956011/149‑1 (MQ=255)                                         gCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGa                                                                                                    <  1:738205/149‑1 (MQ=255)                                             ttCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaa                                                                                                >  1:678819/1‑149 (MQ=255)                                               cAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaaa                                                                                              <  1:713927/149‑1 (MQ=255)                                                aaCCACGAGGCAGTCGGTAAACAGCGTAAACCCTGCCAGCATGATGTTTTCGAGTGGCGGAACTAACTCTTCCCACTGTGGCTCATCACGCATAC‑GGGTCTGGTAGAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                             <  2:1176891/149‑1 (MQ=255)                                                aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCCTCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGCCCGTTGGCACGCAGATAAacac                                                                                             >  2:819989/1‑149 (MQ=255)                                                aaCCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATaaaatc                                                                                             >  1:447871/1‑147 (MQ=255)                                                  cccccagccagTCGGTAAACAGCGTAAACCCTGCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCAGCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGa                                                                                           <  1:1100524/146‑1 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                                                          >  1:772062/1‑149 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGAt                                                                                          >  1:527050/1‑149 (MQ=255)                                                   caccaGCCAGTCGGTAAACAGCGTAAACCCTTCAAGCATGATGTTTTCGCGTGGCGGAATTACCTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGAAGATAAAAACAAt                                                                                          >  2:171052/1‑149 (MQ=255)                                                    accaGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATACGGGGTCCGGTATAAGCCAAAGTTTTTGCCGTGACGGGTGGAACGCAGAAAAAAACggt                                                                                          >  2:1061291/1‑147 (MQ=255)                                                      cagccagTAGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCt                                                                                       <  2:355774/149‑1 (MQ=255)                                                       agccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                                                      <  1:1138868/149‑1 (MQ=255)                                                       agccagTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTg                                                                                      <  1:1145930/149‑1 (MQ=255)                                                         ccagcaggTAAACAGCGCAATCCCTCTCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGt                                                                                    <  2:208186/143‑1 (MQ=255)                                                            gTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGTTTGGAACGCAGATAAAAACGATGCTGGTAGt                                                                                 >  2:473794/1‑149 (MQ=255)                                                               ggTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTGTTTGCCGTGGCGGTTGGAACTCAGATAAAAACGATGCTGGTAGTGAt                                                                              >  1:841698/1‑149 (MQ=255)                                                               ggTAAACAGCGTAAACCCCTCCAGCATGATGTTTTCGCGTGGCGCAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGAt                                                                              <  1:36183/149‑1 (MQ=255)                                                                    aCAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGGTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGAGCACgg                                                                         >  2:743162/1‑149 (MQ=255)                                                                       gCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCCCGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTg                                                                      <  2:32791/149‑1 (MQ=255)                                                                             aCCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCg                                                                >  1:670976/1‑149 (MQ=255)                                                                              cccTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGt                                                               <  2:678819/149‑1 (MQ=255)                                                                               ccgtCCAGCAAGATGTTTTCGCGTGGCGGAACTAATGCTTCCCACTGATGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGCTTTTGCCGGGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTg                                                              <  2:37275/146‑1 (MQ=255)                                                                                cTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGACCAAGGCTGTATTCGTGa                                                             >  1:1287625/1‑149 (MQ=255)                                                                                       aatggtgATTTCGGGTGGCGGAATTAGCTTGTCCCACTGTTGCTCATCACTCATAC‑GGGTCCGGTAGAAGCCAAGGGTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGTTCCCTGAGGTACTGGTGATCTTTgc                                                      <  2:823130/148‑1 (MQ=255)                                                                                           atgTTTTCGCGTGGCGGAATTAACTCTTCCCACTGCTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCGGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcg                                                  >  2:425845/1‑149 (MQ=255)                                                                                              ttttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTgcggcgcgg                                               >  1:36605/1‑149 (MQ=255)                                                                                                ttCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCAGTCTACGCCAAAGTTTTTGCCGTGGCGGCTGGAACGCAGATAAAAACGATACAGGTGGTGATCAGCGCTCAATTCGTGATCTTTccggccgcgca                                             >  2:697755/1‑143 (MQ=255)                                                                                                       ggCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAaaac                                      <  2:31686/149‑1 (MQ=255)                                                                                                         cGGAATTAACTCTTCGCGATGTTGCTCATGACGCATAC‑CGGTTCAGTATCAGCCAAGGTGATGGCCGTGGCAGTTGGAACGCAGATAAAAACGAGGCTGGTAGTGATGAAGGCGATATTCGTGAGCTTTGGAGCGCGGGAGAAaagcaa                                    >  2:282134/1‑149 (MQ=255)                                                                                                                 aCTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGTACGCAGATAAAAACGATGCGGGTAGTGATCAAGGCTGTATCCGTGATCTTTGCGGCGCGGCAGAAAAACAACCGGCTCg                            >  2:1090595/1‑149 (MQ=255)                                                                                                                 aCTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAAAAAATGGCTCg                            >  2:105419/1‑149 (MQ=255)                                                                                                                  ctctTCCCACTGTTGCTCATCACGTATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTTCGCCCCGGCAGAAAAACCAACgccccag                           >  2:93788/1‑143 (MQ=255)                                                                                                                    ctTCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGAATAAGCCAAAGTTTTTGCCGTGGCGGTTGGACCGCAGATAAAAACGATGCTGGTAGTGATCAAGGTTCTATTTGTGATCTTTGCTGCGCGGGAGAAAAGCAAACGACTgggaa                         >  2:830211/1‑147 (MQ=255)                                                                                                                     ttCCCACTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGACAAACAAACGGCTCGGCAt                        >  1:134205/1‑149 (MQ=255)                                                                                                                          aCTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                             >  2:974202/1‑75 (MQ=255)                                                                                                                          aCTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAAc                                                                                             <  1:974202/75‑1 (MQ=255)                                                                                                                          aCTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaa                                         <  2:391356/127‑1 (MQ=255)                                                                                                                          aCTGTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGaaa                                         >  1:391356/1‑127 (MQ=255)                                                                                                                             gTTGCTCATCACGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCAtt                <  2:670976/149‑1 (MQ=255)                                                                                                                                  tcatcaCGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCcgc           >  1:1053211/1‑149 (MQ=255)                                                                                                                                     tcaCGCATAC‑GGGTCCGTGATAGGCCAAAGTTTTGGCCGTGGCGGTTGTAACGCAGATAAAAACGATGCTGGTAGTGTTCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTc        <  2:92239/149‑1 (MQ=255)                                                                                                                                     tcaCGCATAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTc        >  1:812678/1‑149 (MQ=255)                                                                                                                                          cctAC‑GGGCCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagca   <  1:361041/147‑1 (MQ=255)                                                                                                                                          cATAC‑GGGCCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTCATCAGATCGTTATGCGGGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagca   <  2:878108/149‑1 (MQ=255)                                                                                                                                           aTAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCcagcag  >  1:826535/1‑149 (MQ=255)                                                                                                                                           aTAC‑GGGTCCGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCCTCcagcac  >  2:194803/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                             TTGCGGTTAATTTGGCGCAAACTGGTTAACCCGCGCTGACGCTCTTCAACCACCAGCCAGTCGGTAAACAGCGTAAACCCTTCCAGCATGATGTTTTCGCGTGGCGGAATTAACTCTTCCCACTGTTGCTCATCACGCATAC‑GGGTACGGTATAAGCCAAAGTTTTTGCCGTGGCGGTTGGAACGCAGATAAAAACGATGCTGGTAGTGATCAAGGCTGTATTCGTGATCTTTGCGGCGCGGCAGAAAAACAAACGGCTCGGCATCGGCCATTTCCGCGTCCAGCAG  >  NC_000913_3_bme_pgi/1927803‑1928089

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 8 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑