Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	40,952	T→C	6.4%	W162R (TGG→CGG)	oprP →	porin P

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	40,952	0	T	C	6.4%	181.2 / 2.2	62	W162R (TGG→CGG)	oprP	porin P
Reads supporting (aligned to +/- strand):  ref base T (24/34);  new base C (2/2);  total (26/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.01e-01

CTTCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGT  >  NC_002947/40804‑41091                                                                                                                                                       |                                                                                                                                            cTTCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACt                                                                                                                                             >  5:118175/1‑149 (MQ=255) cTTCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACt                                                                                                                                             <  8:171794/149‑1 (MQ=255) cTTCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACt                                                                                                                                             <  5:203995/149‑1 (MQ=255) cTTCGACGAAGCCAGCCTGACCTACACGGGTATCAAGCCGCTGC‑AACTGCGCTTTGGACGTTTCTACACCGACTACGGTCTGGAGTAAGCAACCAGTTCGAAGTTGA‑CAACCGACAAGGCGCGGAACCTGTCTTACGACCTGGCGGCCt                                                                                                                                             >  2:287783/1‑149 (MQ=255)  ttCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGCAAACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACt                                                                                                                                             >  1:282956/1‑149 (MQ=255)    cgacgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGa                                                                                                                                          >  2:304479/1‑149 (MQ=255)     gacgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCTGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGAt                                                                                                                                         <  3:35513/149‑1 (MQ=255)     gacgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGa                                                                                                                                          >  6:91935/1‑148 (MQ=255)     gacgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGa                                                                                                                                          <  5:252832/148‑1 (MQ=255)     gacgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGa                                                                                                                                          <  8:10683/148‑1 (MQ=255)      acgaAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGAt                                                                                                                                         >  7:33696/1‑148 (MQ=255)         acgccagccTGACCTACACGGGTGTCAAGCCGGTGC‑AACTGCGCGTGGGCCGTGTCTACCCCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATTGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAAc                                                                                                                                     <  6:380307/147‑1 (MQ=255)         aaGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAAc                                                                                                                                     <  8:33696/149‑1 (MQ=255)          agccagccTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGACCCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAAc                                                                                                                                     <  1:304479/149‑1 (MQ=255)          agccagccTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACg                                                                                                                                    <  8:325644/149‑1 (MQ=255)           gccagccTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGa                                                                                                                                   >  1:74902/1‑149 (MQ=255)           gccagccTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGa                                                                                                                                   <  7:50245/149‑1 (MQ=255)                      cTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCa                                                                                                                        >  8:157225/1‑149 (MQ=255)                       tACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCAt                                                                                                                       >  1:318990/1‑149 (MQ=255)                         cacGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATgg                                                                                                                     <  3:294435/149‑1 (MQ=255)                         cacGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATgg                                                                                                                     <  6:118175/149‑1 (MQ=255)                              gTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGtt                                                                                                                >  3:338182/1‑149 (MQ=255)                                       cGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGccacc                                                                                                       <  8:331134/149‑1 (MQ=255)                                         gTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGAGTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAg                                                                                                     >  5:352150/1‑149 (MQ=255)                                          tGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGc                                                                                                    <  6:352150/149‑1 (MQ=255)                                               aCTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCa                                                                                                >  3:163204/1‑149 (MQ=255)                                                cTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCa                                                                                                >  2:97684/1‑148 (MQ=255)                                                     cTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTg                                                                                          <  2:282956/149‑1 (MQ=255)                                                        tGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGa                                                                                         <  2:318990/147‑1 (MQ=255)                                                           ccGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATgg                                                                                    <  5:91935/149‑1 (MQ=255)                                                            cGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGc                                                                                   >  2:302626/1‑149 (MQ=255)                                                              tttCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCt                                                                                  >  8:128959/1‑148 (MQ=255)                                                                tCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTAt                                                                               >  2:397373/1‑149 (MQ=255)                                                                  tACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGt                                                                             >  2:266928/1‑149 (MQ=255)                                                                    cacaGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGttt                                                                           >  7:353146/1‑149 (MQ=255)                                                                    cacCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGtt                                                                            >  8:112435/1‑148 (MQ=255)                                                                      ccGACTTCGGTCTTGAGATAGCCACAAGTTCGAAGTTGA‑CCACCGCCATGGAGCGTAACCTCTAGTACGACCTGGCCGACTGGATCAACGACAACGCCGGCATGTGTGTTCAAGCCAACAGAACCATCGCTGACATGTCTTATGTTTcc                                                                         <  8:296220/149‑1 (MQ=255)                                                                      ccGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTcc                                                                         <  1:65109/149‑1 (MQ=255)                                                                      ccGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTcc                                                                         <  4:105617/149‑1 (MQ=255)                                                                          cTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCgg                                                                       >  2:26833/1‑147 (MQ=255)                                                                          cTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTa                                                                     <  2:74902/149‑1 (MQ=255)                                                                                    gagaAAGCCACCAGGTCGAAGTGGA‑CCCCCGGCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCa                                                            <  4:338182/148‑1 (MQ=255)                                                                                     agaAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGc                                                          >  5:225901/1‑149 (MQ=255)                                                                                         aGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGaaa                                                      <  1:97684/149‑1 (MQ=255)                                                                                          gCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAAc                                                     <  8:353146/149‑1 (MQ=255)                                                                                           ccaccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACCGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGt                                                                      >  8:309051/1‑131 (MQ=255)                                                                                           ccaccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACCGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGt                                                                      <  7:309051/131‑1 (MQ=255)                                                                                           ccaccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACCGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGt                                                                      <  7:314773/131‑1 (MQ=255)                                                                                           ccaccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACCGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGt                                                                      >  8:314773/1‑131 (MQ=255)                                                                                            caccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAAcc                                                    >  2:342487/1‑148 (MQ=255)                                                                                             accaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATTGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCaa                                                  <  1:266928/149‑1 (MQ=255)                                                                                             accaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCgg                                                                       <  6:343968/128‑1 (MQ=255)                                                                                             accaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCgg                                                                       >  5:343968/1‑128 (MQ=255)                                                                                              ccaGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAAt                                                 <  7:157225/149‑1 (MQ=255)                                                                                               caGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATa                                                <  7:128959/149‑1 (MQ=255)                                                                                                        gTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCg                                        <  4:307442/148‑1 (MQ=255)                                                                                                        gTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCg                                        >  3:307442/1‑148 (MQ=255)                                                                                                                   gCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAg                             <  6:225901/149‑1 (MQ=255)                                                                                                                    ccATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGc                            >  6:224332/1‑149 (MQ=255)                                                                                                                     cATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCa                           <  7:112435/149‑1 (MQ=255)                                                                                                                              gAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATaa                                               >  5:226933/1‑120 (MQ=255)                                                                                                                              gAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATaa                                               <  6:226933/120‑1 (MQ=255)                                                                                                                                 ccTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACa               <  5:224332/149‑1 (MQ=255)                                                                                                                                  ctctCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACaa              <  4:163204/149‑1 (MQ=255)                                                                                                                                       gTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAacgta          <  1:26833/148‑1 (MQ=255)                                                                                                                                              cTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCgt  >  4:193353/1‑149 (MQ=255)                                                                                                                                              cTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCgt  <  1:342487/149‑1 (MQ=255)                                                                                                                                                       tGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAg                                                                    >  7:300081/1‑74 (MQ=255)                                                                                                                                                       tGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAg                                                                    <  8:300081/74‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                            CTTCGACGAAGCCAGCCTGACCTACACGGGTTTCAAGCCGGTGC‑AACTGCGCTTTGGCCGTTTCTACACCGACTTCGGTCTTGAGAAAGCCACCAGTTCGAAGTGGA‑CCACCGCCATGGAGCGGAACCTCTCGTACGACCTGGCGGACTGGATCAACGACAACGCCGGCATGGGTGTTCAAGCCACCAGCACCATCGCTGACATGGCTTATGTTTCCGGTAGCGTTTTCAGCGAAACCAATAACAACTCGGATGGCGACAGCACCAAGCGCTACAACGTACGTGGCGT  >  NC_002947/40804‑41091

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