Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,331,580	G→T	6.7%	G94C (GGC→TGC)	PP_3803 →	cation ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,331,580	0	G	T	6.7%	186.4 / 3.1	60	G94C (GGC→TGC)	PP_3803	cation ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base T (2/2);  total (32/28)
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.69e-01

TGGGCACCTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCT  >  NC_002947/4331434‑4331728                                                                                                                                                    |                                                                                                                                                     tGGGCACCTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCgg                                                                                                                                                     >  3:157426/1‑148 (MQ=255)      aCCTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCt                                                                                                                                               <  1:58920/149‑1 (MQ=255)        cTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGc                                                                                                                                             >  4:44589/1‑149 (MQ=255)        cTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGc                                                                                                                                             >  8:126053/1‑149 (MQ=255)                  gtCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTcc                                                                                                                                    <  7:126053/148‑1 (MQ=255)                   tCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTg                                                                                                                                  >  3:156265/1‑149 (MQ=255)                              gCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTa                                                                                                                        <  6:129381/148‑1 (MQ=255)                                     tctcTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTATCCCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCt                                                                                                                 >  1:81373/1‑149 (MQ=255)                                       tctGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTc                                                                                                              <  6:300503/149‑1 (MQ=255)                                          gAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGcc                                                                                                           >  2:219055/1‑149 (MQ=255)                                                 tCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTCCGCAACCGTACCGGGCTGTCGCCGGATaca                                                                                                    >  7:272447/1‑149 (MQ=255)                                                   ggCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATacac                                                                                                   <  5:18972/148‑1 (MQ=255)                                                       aTGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGAtt                                                                                              <  2:81373/149‑1 (MQ=255)                                                            gcgcTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCgt                                                                                         >  8:231774/1‑149 (MQ=255)                                                               cTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGtt                                                                                      >  8:201735/1‑149 (MQ=255)                                                                 gACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTg                                                                                                 >  6:337763/1‑136 (MQ=255)                                                                 gACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTg                                                                                                 <  5:337763/136‑1 (MQ=255)                                                                     ggCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCTGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGc                                                                                           >  4:72224/1‑138 (MQ=255)                                                                     ggCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCTGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGc                                                                                           <  3:72224/138‑1 (MQ=255)                                                                      gCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCg                                                                               >  2:261042/1‑149 (MQ=255)                                                                      gCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCg                                                                               <  4:157426/149‑1 (MQ=255)                                                                       cGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCgg                                                                              <  1:219055/149‑1 (MQ=255)                                                                             cATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTc                                                                                <  5:162533/141‑1 (MQ=255)                                                                             cATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTc                                                                                >  6:162533/1‑141 (MQ=255)                                                                                  gTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGgcg                                                                   >  3:17985/1‑149 (MQ=255)                                                                                  gTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGgcg                                                                   >  6:15446/1‑149 (MQ=255)                                                                                     t‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCtt                                                                <  4:58161/149‑1 (MQ=255)                                                                                         tgctgGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCTGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGgc                                                             >  1:308450/1‑149 (MQ=255)                                                                                              ggTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCTGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAg                                                        <  2:308450/149‑1 (MQ=255)                                                                                              ggTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCa                                                         <  7:231774/148‑1 (MQ=255)                                                                                              ggTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAg                                                        >  4:26246/1‑149 (MQ=255)                                                                                               gTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCCAAGCGCCAGc                                                       <  4:17985/149‑1 (MQ=255)                                                                                                     ccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                 >  1:255469/1‑149 (MQ=255)                                                                                                     ccTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCctgctg                                                 >  7:322930/1‑149 (MQ=255)                                                                                                          aCCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAg                                                        >  6:132952/1‑137 (MQ=255)                                                                                                          aCCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAg                                                        <  5:132952/137‑1 (MQ=255)                                                                                                                ccTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTgg                                       >  3:22051/1‑148 (MQ=255)                                                                                                                 cTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggca                                     <  7:201735/149‑1 (MQ=255)                                                                                                                 cTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggca                                     <  2:255469/149‑1 (MQ=255)                                                                                                                  ttACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCct                                                     >  4:282310/1‑132 (MQ=255)                                                                                                                  ttACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCct                                                     <  3:282310/132‑1 (MQ=255)                                                                                                                  ttACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGAGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcag                                    >  1:282765/1‑149 (MQ=255)                                                                                                                    aCGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCggg                                                                                                                   <  4:271461/68‑1 (MQ=255)                                                                                                                    aCGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCggg                                                                                                                   >  3:271461/1‑68 (MQ=255)                                                                                                                      ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGc                                                       >  4:3192/1‑126 (MQ=255)                                                                                                                      ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGc                                                       <  3:3192/126‑1 (MQ=255)                                                                                                                      ggCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTggcaggc                                  >  1:27049/1‑147 (MQ=255)                                                                                                                       gCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAg                               <  3:44589/149‑1 (MQ=255)                                                                                                                        cAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                              <  4:22051/149‑1 (MQ=255)                                                                                                                               ttCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGt                       <  1:261042/149‑1 (MQ=255)                                                                                                                                   gCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGcac                   <  2:282765/149‑1 (MQ=255)                                                                                                                                     tACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACAt                 >  2:217241/1‑149 (MQ=255)                                                                                                                                       ctgcctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                              <  3:9679/134‑1 (MQ=255)                                                                                                                                       ctgcctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGa                              >  4:9679/1‑134 (MQ=255)                                                                                                                                          cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCg            <  2:81155/149‑1 (MQ=255)                                                                                                                                          cctgcTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGCGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTTGCAGGCAAGATCAACGTGCACATCCTc             >  4:122049/1‑148 (MQ=255)                                                                                                                                            tgctgTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTc             >  7:170018/1‑146 (MQ=255)                                                                                                                                             gctgTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACTCTGATTGGCCTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGAACGATACACGTGCGCATCCACggg          >  5:278567/1‑146 (MQ=255)                                                                                                                                              ctgTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGaa        >  8:168911/1‑149 (MQ=255)                                                                                                                                                gTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACg      >  8:297628/1‑149 (MQ=255)                                                                                                                                                   cGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTg    <  5:15446/148‑1 (MQ=255)                                                                                                                                                    ggCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCt  >  6:341693/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TGGGCACCTTGGTGGTGGTCAAGCGCTTTGCGTTCTTCTCTGAAGCGGTCGGCCATGCGGCGCTGACCGGCGTGGCCATCGGTAT‑CCTGCTGGGTGAGCCCTATACCGGCCCTTACGGCAGCCTGTTCGGCTACTGCCTGCTGTTCGGCATCCTGCTCAACTTCCTGCGCAACCGTACCGGGCTGTCGCCGGATACACTGATTGGCGTGTTCCTGTCGGTGTCGCTGGCGCTTGGCGCCAGCCTGCTGCTGATGCTGGCAGGCAAGATCAACGTGCACATCCTCGAGAACGTGCT  >  NC_002947/4331434‑4331728

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