Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	936,058	G→A	6.1%	G3188G (GGC→GGT)	PP_0806 ←	surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	936,058	0	G	A	6.1%	198.3 / 2.2	65	G3188G (GGC→GGT)	PP_0806	surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base A (2/2);  total (34/31)
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.30e-01

GGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTC  >  NC_002947/935911‑936205                                                                                                                                                    |                                                                                                                                                    ggTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCGGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCCCAGTCCCGTCGGCGGGGCTGACCTGCCCGGTGGTGc                                                                                                                                                    >  8:319764/1‑149 (MQ=255)       gCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGtt                                                                                                                                              >  2:45007/1‑149 (MQ=255)       gCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCCGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGt                                                                                                                                               <  3:62329/148‑1 (MQ=255)        ccAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTc                                                                                                                                             <  8:233358/149‑1 (MQ=255)        ccAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTc                                                                                                                                             >  4:262178/1‑149 (MQ=255)                 tCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGAcc                                                                                                                                    >  5:73359/1‑149 (MQ=255)                    cgcAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGt                                                                                                                                 <  3:262178/149‑1 (MQ=255)                    cgcAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCAGTTCGCCGCGACCGGt                                                                                                                                 >  6:170270/1‑149 (MQ=255)                         gACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCa                                                                                                                            <  6:295803/149‑1 (MQ=255)                            cACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGc                                                                                                                         <  4:318452/149‑1 (MQ=255)                            cACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGc                                                                                                                         >  4:236185/1‑149 (MQ=255)                             aCGTCAAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCt                                                                                                                        >  7:156380/1‑149 (MQ=255)                              cGTCCAATTACTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTg                                                                                                                       >  1:93473/1‑149 (MQ=255)                               gTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTACCCGGTTCGCCGCGACCGGTGAGCAGGCTGc                                                                                                                      >  6:78757/1‑149 (MQ=255)                                   aaTTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGc                                                                                                                             >  4:53606/1‑138 (MQ=255)                                   aaTTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGc                                                                                                                             <  3:53606/138‑1 (MQ=255)                                   aaTTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATc                                                                                                                  >  6:219996/1‑149 (MQ=255)                                     ttCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                <  1:88427/149‑1 (MQ=255)                                     ttCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                <  5:343623/149‑1 (MQ=255)                                     ttCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCAGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                <  5:170270/149‑1 (MQ=255)                                     ttCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGGTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                <  3:32347/149‑1 (MQ=255)                                         tGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACt                                                                                                            >  7:258616/1‑149 (MQ=255)                                           gCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCa                                                                                                          <  4:315399/149‑1 (MQ=255)                                              gTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTACCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGtt                                                                                                       <  5:78757/149‑1 (MQ=255)                                               tCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTg                                                                                                      >  2:20959/1‑149 (MQ=255)                                               tCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTg                                                                                                      <  5:219996/149‑1 (MQ=255)                                               tCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTg                                                                                                      <  3:236185/149‑1 (MQ=255)                                                       cgcTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTACCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGc                                                                                              >  3:33246/1‑149 (MQ=255)                                                        gcTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCt                                                                                             <  1:20959/149‑1 (MQ=255)                                                              ttCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                >  1:192264/1‑124 (MQ=255)                                                              ttCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCgg                                                                                                                <  2:192264/124‑1 (MQ=255)                                                                 aGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCAt                                                                                     <  8:258616/148‑1 (MQ=255)                                                                   gACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTgg                                                                                  <  4:320570/149‑1 (MQ=255)                                                                   gACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTgg                                                                                  >  3:253446/1‑149 (MQ=255)                                                                      cACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGAc                                                                                                                                     >  7:54836/1‑95 (MQ=255)                                                                         gCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCgggg                                                                            <  4:253446/149‑1 (MQ=255)                                                                            gATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGc                                                                                                     <  6:238679/121‑1 (MQ=255)                                                                            gATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGc                                                                                                     >  5:238679/1‑121 (MQ=255)                                                                               cacaCCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTggg                                                                                 <  5:231847/138‑1 (MQ=255)                                                                               cacaCCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTggg                                                                                 >  6:231847/1‑138 (MQ=255)                                                                                acacCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGGGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGt                                                                           >  5:213418/1‑143 (MQ=255)                                                                                 cacCGACGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCgg                                                                                                                                                <  5:333228/73‑1 (MQ=255)                                                                                 cacCGACGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCgg                                                                                                                                                >  6:333228/1‑73 (MQ=255)                                                                                  acCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCgg                                                                   >  2:201378/1‑149 (MQ=255)                                                                                   ccGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCggg                                                                  <  2:93473/149‑1 (MQ=255)                                                                                        ggCACCGACCGGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCgg                                                             <  2:342061/149‑1 (MQ=255)                                                                                               aCCTGGGCGTTGCCCAGCACAGTCCCGTGGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCTGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTg                                                      >  1:192969/1‑149 (MQ=255)                                                                                                ccTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATGTgg                                                     >  1:173125/1‑149 (MQ=255)                                                                                                ccTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGCGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTgg                                                     >  4:60027/1‑149 (MQ=255)                                                                                                   gggCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTggg                                                    >  7:55964/1‑147 (MQ=255)                                                                                                   gggCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTggg                                                    >  7:60602/1‑147 (MQ=255)                                                                                                       gTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCa                                              <  6:73359/149‑1 (MQ=255)                                                                                                           cccAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCg                                           >  7:17877/1‑148 (MQ=255)                                                                                                               gCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTg                                      >  7:83870/1‑149 (MQ=255)                                                                                                                 acaGTCCCTGCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGGTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGtt                                    <  5:188464/149‑1 (MQ=255)                                                                                                                                cTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAAc                     >  3:129446/1‑149 (MQ=255)                                                                                                                                  gatgaCCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCg                   >  7:19833/1‑149 (MQ=255)                                                                                                                                    tgaCCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGGACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGcc                 >  1:36241/1‑149 (MQ=255)                                                                                                                                      aCCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCa                <  6:307861/148‑1 (MQ=255)                                                                                                                                      aCCTGCACGGTGGTACCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCa                <  4:33246/148‑1 (MQ=255)                                                                                                                                        cTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGa             >  3:77032/1‑149 (MQ=255)                                                                                                                                        cTGCACGGTGGCGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGa             <  3:60027/149‑1 (MQ=255)                                                                                                                                          gCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAgc                              >  7:157888/1‑130 (MQ=255)                                                                                                                                          gCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAgc                              <  8:157888/130‑1 (MQ=255)                                                                                                                                           cACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGAc            <  4:77032/147‑1 (MQ=255)                                                                                                                                              tgtggtgCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGAAGTTCGCTGGATTCAGTTGGCGCGGTTATATCGGGGGATGTGATCTGTGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCCCGACATc         >  5:252579/2‑147 (MQ=255)                                                                                                                                                   tgCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGtc  >  7:63575/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                    GGTGTTGCCAGCGGGGTCGCGCAGGACCACGTCCAATTCCTGGCCGTCGGTCTGCGCTGGGTTCAGGACCACGCTGATCACACCGTCGGCACCGACCTGGGCGTTGCCCAGCACAGTCCCGTCGGCGCTGATGACCTGCACGGTGGTGCCCGGTTCGCCGCGACCGGTGAGCAGGCTGCCATCGGCACTCAGTTGCAGTTCGCTGGCTTCATTGGGCGGGGTGATATCGGGGGCGGTGATCTGGGTAGGCAGCGAACTGTTGCCAGCGCCGTCAACCGCCACGACATCCAGGGTC  >  NC_002947/935911‑936205

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