Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,447,753	A→G	10.6%	A173A (GCA→GCG)	PP_5556 →	bacteriophage lysis protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,447,753	0	A	G	10.6%	173.9 / 9.9	66	A173A (GCA→GCG)	PP_5556	bacteriophage lysis protein
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base G (3/4);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 6.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

CTTCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGAATATCTT  >  NC_002947/3447606‑3447901                                                                                                                                                    |                                                                                                                                                      cTTCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCAc                                                                                                                                                      >  7:290560/1‑149 (MQ=255)  ttCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGGCCGAGGACTGATCGCACt                                                                                                                                                     >  6:366923/1‑149 (MQ=255)  ttCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACt                                                                                                                                                     >  7:159088/1‑149 (MQ=255)  ttCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACt                                                                                                                                                     >  4:112237/1‑149 (MQ=255)   tCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGACGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTAAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACt                                                                                                                                                     >  7:127064/1‑148 (MQ=255)     gccgccTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCa                                                                                                                                                  <  2:292189/149‑1 (MQ=255)      ccgccTCGGGTTGTGACAGTGGGGAGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAg                                                                                                                                                 >  4:193256/1‑149 (MQ=255)       cgccTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAgg                                                                                                                                                <  5:364921/149‑1 (MQ=255)               gTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGcc                                                                                                                                          >  6:245341/1‑147 (MQ=255)                  gtgACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCAGGCCTGCCAAGc                                                                                                                                      <  6:347189/148‑1 (MQ=255)                       aGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCt                                                                                                                                <  6:144670/149‑1 (MQ=255)                       aGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCt                                                                                                                                >  1:346348/1‑149 (MQ=255)                                cGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCAGGCCTGCCAAGCCTATCTGCGCGAAAt                                                                                                                       >  6:256646/1‑149 (MQ=255)                                 gagaAGCCACCGGCACCGGTAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAAtc                                                                                                                      <  5:28640/149‑1 (MQ=255)                                 gagaAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAAt                                                                                                                       >  6:67782/1‑148 (MQ=255)                                    aaGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATctct                                                                                                                   <  5:245341/149‑1 (MQ=255)                                        caccggcaccggAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctgact                                                                                                               >  6:192728/1‑149 (MQ=255)                                          ccggcaccggAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctgactg                                                                                                              >  3:301849/1‑148 (MQ=255)                                          ccggcaccggAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATAACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctgactg                                                                                                              >  3:207676/1‑148 (MQ=255)                                                      gCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGc                                                                                                 >  8:1549/1‑148 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGTCCTGCCAAGCGTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCaa                                                                                              <  3:193256/149‑1 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctgac                                                                                                                >  5:13832/1‑132 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctgac                                                                                                                <  6:13832/132‑1 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctga                                                                                                                 >  3:294527/1‑131 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTctga                                                                                                                 <  4:294527/131‑1 (MQ=255)                                                        gTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCaa                                                                                              <  3:112237/149‑1 (MQ=255)                                                          ggTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGc                                                                                            >  3:355021/1‑149 (MQ=255)                                                               aTGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTcgc                                                                                        <  8:290560/148‑1 (MQ=255)                                                               aTGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATAACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTcgcg                                                                                       <  4:207676/149‑1 (MQ=255)                                                                               ccAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAg                                                                       <  2:346348/149‑1 (MQ=255)                                                                                 aaCTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAgtt                                                                     <  7:1549/149‑1 (MQ=255)                                                                                   cTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAgttgt                                                                   >  5:165020/1‑149 (MQ=255)                                                                                         cGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATCTCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAgg                                                             <  5:256646/149‑1 (MQ=255)                                                                                                aTGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAAtt                                                      >  8:282132/1‑149 (MQ=255)                                                                                                aTGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAAtt                                                      >  1:213881/1‑149 (MQ=255)                                                                                                   cTCAACGAATTGTAGGCATCAGCGAGTAAGGTGAACGAGGACTGATCGCACTGGAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTcc                                                    <  8:78993/148‑1 (MQ=255)                                                                                                    tCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCgt                                                  >  5:282049/1‑149 (MQ=255)                                                                                                    tCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCgt                                                  >  1:202043/1‑149 (MQ=255)                                                                                                    tCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGATGTTCGAGGGGGAATTCCgt                                                  >  7:222685/1‑149 (MQ=255)                                                                                                      aaCGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTgg                                                                                <  4:198604/117‑1 (MQ=255)                                                                                                      aaCGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCgtgt                                                >  5:20071/1‑149 (MQ=255)                                                                                                      aaCGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCgtgt                                                >  3:90565/1‑149 (MQ=255)                                                                                                      aaCGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCgtg                                                 <  5:192728/148‑1 (MQ=255)                                                                                                         gAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTg                                             <  7:282132/149‑1 (MQ=255)                                                                                                               tCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCGCTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCt                                       >  3:217985/1‑149 (MQ=255)                                                                                                                 gCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCtt                                      <  6:165020/148‑1 (MQ=255)                                                                                                                 gCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGAATGAAA‑ATATCGAGCCAAGCTCTCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCtt                                      <  8:127064/148‑1 (MQ=255)                                                                                                                  ccATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCACGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAg                                    >  2:294480/1‑149 (MQ=255)                                                                                                                  ccATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAg                                    <  8:159088/149‑1 (MQ=255)                                                                                                                  ccATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAg                                    <  5:67782/149‑1 (MQ=255)                                                                                                                   cATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGTGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAgg                                   >  1:310629/1‑149 (MQ=255)                                                                                                                   cATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAgg                                   <  2:202043/149‑1 (MQ=255)                                                                                                                    aTCACCTTCGAAGGTGACCGATGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGt                                  <  8:222685/149‑1 (MQ=255)                                                                                                                    aTCACCGACGAAGGAGGCCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGt                                  <  5:366923/149‑1 (MQ=255)                                                                                                                       aCCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGa                                >  2:328066/1‑148 (MQ=255)                                                                                                                        ccGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACa                              <  4:355021/149‑1 (MQ=255)                                                                                                                         cgacgaAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCaa                                                                                              <  6:101404/84‑1 (MQ=255)                                                                                                                               aGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTg                       <  2:213881/149‑1 (MQ=255)                                                                                                                                      cGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAgggg                                                           >  3:248382/1‑106 (MQ=255)                                                                                                                                      cGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAgggg                                                           <  4:248382/106‑1 (MQ=255)                                                                                                                                      cGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGa                                                         >  2:329988/1‑108 (MQ=255)                                                                                                                                      cGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGa                                                         <  1:329988/108‑1 (MQ=255)                                                                                                                                          gACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATgg            >  6:262048/1‑149 (MQ=255)                                                                                                                                                tCGCGCTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGAata      >  2:113827/1‑149 (MQ=255)                                                                                                                                                 cGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAACATATCGACCCAAGCTCGCGTAATTGGCAAGCAAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGataa      >  1:266769/1‑146 (MQ=255)                                                                                                                                                  gCACTGCAGGCCTGCCAAGCCTATCTGAGCGAAATCTCTGACTGAAA‑ATATCCAGCCAAGCTCGCGGAATTCGCACTCAAAGTTGTTCGAGGGGGAATTCAGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGAatata    >  7:11422/1‑148 (MQ=255)                                                                                                                                                    aCTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGCGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGAATAtctt  <  4:301849/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CTTCGCCGCCTCGGGTTGTGACAGTGGGGTGCGAGAAGCCACCGGCACCGGAAGCGTGGTTCATGGCGCCGTACGCGCCCAACTTGACCGAGCGCATGCTCAACGAATTGTCGCCATCACCGACGAAGGTGACCGAGGACTGATCGCACTGCAGGCCTGCCAAGCCTATCTGCGCGAAATCTCTGACTGAAA‑ATATCGAGCCAAGCTCGCGGAATTGGCACTCCAAGTTGTTCGAGGGGGAATTCCGTGTTTGGGGTCTTAGGTGACAGTCACTGGCGGCCCATGGGGAATATCTT  >  NC_002947/3447606‑3447901

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