Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,633,942	A→G	6.9%	S1158G (AGC→GGC)	putA →	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,633,942	0	A	G	6.9%	168.0 / 2.4	59	S1158G (AGC→GGC)	putA	proline dehydrogenase/1‑pyrroline‑5‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (27/27);  new base G (2/2);  total (30/29)
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.97e-01

CCCGGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCCGCGG  >  NC_002947/5633794‑5634089                                                                                                                                                     |                                                                                                                                                     cccgGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCa                                                                                                                                                      <  1:129806/149‑1 (MQ=255)           gCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gc                                                                                                                                           <  1:141361/149‑1 (MQ=255)           gCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gc                                                                                                                                           >  8:35080/1‑149 (MQ=255)             cATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑gcca                                                                                                                                         <  6:113435/149‑1 (MQ=255)                          ttCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGG‑CGAACTGGCTGCACAGCGCGGCCAGGTCGgccagctg                                                                                                                                     <  3:326682/140‑3 (MQ=37)                          ttCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGG‑CGAACTGGCTGCACAGCGCGGCCAGGTCGgccagctg                                                                                                                                     >  4:326682/1‑138 (MQ=37)                            ccAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCCCGAACAGCTGGTAACGCCGCTGCACGGCCTGACGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGACGCCAGTTC‑GCCAGCCAGACGCAGAGCg                                                                                                                          >  3:216379/1‑149 (MQ=255)                              agcagcagGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCgg                                                                                                                         >  3:217281/1‑148 (MQ=255)                                       aCGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATcgcccgc                                                                                                               >  3:25138/1‑149 (MQ=255)                                           cGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCctgc                                                                                                           >  6:599/1‑149 (MQ=255)                                            gAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCCGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCctgct                                                                                                          <  5:269121/149‑1 (MQ=255)                                              aGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCtgctgc                                                                                                        <  6:169929/149‑1 (MQ=255)                                                     ccGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGccc                                                                                                 <  7:35080/149‑1 (MQ=255)                                                     ccGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGc                                                                                                   <  1:142959/147‑1 (MQ=255)                                                      cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAg                                                                                                                                                 <  1:50462/101‑1 (MQ=255)                                                      cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAg                                                                                                                                                 >  2:50462/1‑101 (MQ=255)                                                      cGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGccc                                                                                                 >  5:156706/1‑148 (MQ=255)                                                         aCCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGAcc                                                                                             >  5:91620/1‑149 (MQ=255)                                                            gCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGc                                                                                          <  5:599/149‑1 (MQ=255)                                                              aCCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGa                                                                                        <  4:72390/149‑1 (MQ=255)                                                                ccTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGGgc                                                                                      >  8:90049/1‑149 (MQ=255)                                                                ccTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                                      >  6:94963/1‑149 (MQ=255)                                                                 cTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCgg                                                                                           >  3:41053/1‑143 (MQ=255)                                                                 cTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCgg                                                                                           <  4:41053/143‑1 (MQ=255)                                                                      cGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                                      >  8:276190/1‑143 (MQ=255)                                                                      cGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAgc                                                                                      <  7:276190/143‑1 (MQ=255)                                                                        aaCAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGc                                                                              <  6:156706/149‑1 (MQ=255)                                                                          cAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTa                                                                            <  1:11086/149‑1 (MQ=255)                                                                             cTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTacac                                                                         <  5:94963/149‑1 (MQ=255)                                                                                       cGCTGCACGACCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCcg                                                               >  5:59093/1‑149 (MQ=255)                                                                                               ggCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCa                                                                                                                               >  1:314089/1‑77 (MQ=255)                                                                                               ggCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCa                                                                                                                               <  2:314089/77‑1 (MQ=255)                                                                                                     aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCAc                                                      <  8:219115/144‑1 (MQ=255)                                                                                                     aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCAc                                                      >  7:219115/1‑144 (MQ=255)                                                                                                     aaGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCACAACAGCTACACCATTCTGCCGCGTGAGCACGTg                                                   <  2:70183/147‑1 (MQ=255)                                                                                                        gCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCtgtg                                              >  1:333926/1‑149 (MQ=255)                                                                                                           tGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGt                                                            >  5:136279/1‑132 (MQ=255)                                                                                                           tGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGt                                                            <  6:136279/132‑1 (MQ=255)                                                                                                             ggCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTg                                          <  6:59093/148‑1 (MQ=255)                                                                                                              gCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                                        <  4:25138/149‑1 (MQ=255)                                                                                                                     acaacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAAc                                 <  6:91620/149‑1 (MQ=255)                                                                                                                      caacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAAc                                 >  7:240746/1‑148 (MQ=255)                                                                                                                       aacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                                        <  5:184756/140‑1 (MQ=255)                                                                                                                       aacCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGc                                        >  6:184756/1‑140 (MQ=255)                                                                                                                              tGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGa                         >  3:198841/1‑148 (MQ=255)                                                                                                                                gCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAc                        <  4:217281/147‑1 (MQ=255)                                                                                                                                gCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACct                      >  8:216409/1‑149 (MQ=255)                                                                                                                                 gcGAACTGGCTGCACAGCGCGGCCAGGTCGGCCAGCTAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAc                        <  7:201295/146‑1 (MQ=38)                                                                                                                                   gaacTGGCTGCACAGCGCGGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACctgct                   >  5:99084/4‑149 (MQ=255)                                                                                                                                    aCCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACctgctg                  >  5:253077/1‑149 (MQ=255)                                                                                                                                      cTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgc                >  1:350341/1‑149 (MQ=255)                                                                                                                                         gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAcc                       <  4:65816/139‑1 (MQ=255)                                                                                                                                         gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGAcc                       >  3:65816/1‑139 (MQ=255)                                                                                                                                         gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgcg               >  2:246821/1‑147 (MQ=255)                                                                                                                                         gCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGgcg               >  3:172187/1‑147 (MQ=255)                                                                                                                                          ccGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGAGGCATCGCCCTCCTGCTGCCAGGCCCGACCGGCGATCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAg            <  4:216379/149‑1 (MQ=255)                                                                                                                                           cgcgCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGAGCAGt           >  6:252194/1‑149 (MQ=255)                                                                                                                                                  tgCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTcgccg     <  1:246821/148‑1 (MQ=255)                                                                                                                                                    cAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCcgcg   <  2:333926/148‑1 (MQ=255)                                                                                                                                                    cAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCCGCgg  >  7:329100/1‑149 (MQ=255)                                                                                                                                                     aGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCcgcg   >  6:80548/1‑147 (MQ=255)                                                                                                                                                     aGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCcgcg   <  5:80548/147‑1 (MQ=255)                                                                                                                                                     aGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCCGCgg  >  4:167573/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                     CCCGGCCGACGCCATTGGCCGCCACTTCCAGCAGCAGGACGGCGAAGGCACGCCGGACCGCACCCTGCACGAACAGCTGGTAAAGCCGCTGCACGGCCTGAAGGCCTGGGCCGAGAACAACCAGCTGGCCGACCTGGCCGCGCTGTGCAGCCAGTTC‑GCCAGCCAGTCGCAGAGCGGCATCGCCCGCCTGCTGCCAGGCCCGACCGGCGAGCGCAACAGCTACACCATTCTGCCGCGTGAGCACGTGCTGTGCCTGGCCGACAACGAAACCGACCTGCTGGCGCAGTTCGCCGCGG  >  NC_002947/5633794‑5634089

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