Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,364,420	C→G	8.8%	G72R (GGT→CGT)	folP ←	dihydropteroate synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,364,420	0	C	G	8.8%	241.5 / 9.8	79	G72R (GGT→CGT)	folP	dihydropteroate synthase
Reads supporting (aligned to +/- strand):  ref base C (37/34);  new base G (3/4);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 7.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

CCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364277‑5364568                                                                                                                                                |                                                                                                                                                      cccGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGgc                                                                                                                                                  <  7:256267/149‑1 (MQ=255)  ccGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGgcg                                                                                                                                                 >  6:227132/1‑149 (MQ=255)     gCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGATGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                              <  1:223265/149‑1 (MQ=255)     gCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGt                                                                                                                                              <  1:228490/149‑1 (MQ=255)       ccccAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCGGGTGCGCGTgc                                                                                                                                            <  6:266567/149‑2 (MQ=255)        cccAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCGGGTGCGCGTgca                                                                                                                                           <  2:246160/149‑3 (MQ=255)        cccAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGa                                                                                                                                           >  7:299550/1‑149 (MQ=255)                 gcgGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCAcc                                                                                                                                  >  2:209736/1‑149 (MQ=255)                        cACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGatat                                                                                                                             <  5:207483/147‑1 (MQ=255)                         aCGCATGACAGCAGGCGTGGAGGTATCGAACGAAATGACCACGTCGAGGCGACTGTTGATGGACTCGACCATTGGCGCCACACGATCCAGCTCCTCGGTCACCGAGACTGCACGCGCACACGGTCGCGTGGACTAGCCACCGATATCGa                                                                                                                          >  1:248849/1‑149 (MQ=255)                           gCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCGGGTGCGCGTgcagt                                                                                                                                         <  8:161215/132‑5 (MQ=255)                           gCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCCCGGGTGCGCGTgcagt                                                                                                                                         >  7:161215/1‑128 (MQ=255)                             aTGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCaa                                                                                                                      >  8:57682/1‑149 (MQ=255)                             aTGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCaa                                                                                                                      >  7:210752/1‑149 (MQ=255)                             aTGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCaa                                                                                                                      >  5:51405/1‑149 (MQ=255)                               gACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAAc                                                                                                                     <  4:138753/148‑1 (MQ=255)                               gACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAAc                                                                                                                     >  3:138753/1‑148 (MQ=255)                                  agcagGCGTTGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTTATGGCGTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTg                                                                                                                 <  3:80174/149‑1 (MQ=255)                                    cagGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGACACACGCTCCAGCTCCTCGGTAACCGAGACTGCACGCGCTCCCGGGCGCGTGGACTCGGCAGCGATATAGATAAACGttcag                                                                                                               >  8:273531/1‑146 (MQ=255)                                     agGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTgcgc                                                                                                              >  3:260718/1‑149 (MQ=255)                                               gTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCaa                                                                                                     <  8:299550/148‑1 (MQ=255)                                                 aTCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑Ca                                                                                                 >  4:140704/1‑149 (MQ=255)                                                  tCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGtt                                                                                                                  <  8:316333/132‑1 (MQ=255)                                                  tCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGtt                                                                                                                  >  7:316333/1‑132 (MQ=255)                                                    gACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATc                                                                                               <  7:47287/148‑1 (MQ=255)                                                     aCCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGc                                                                                             >  5:357928/1‑149 (MQ=255)                                                     aCCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTCGACTCGCCACCGATAACCATCAACGTAGCGCCGGCGGCAAC‑CATCGc                                                                                             >  6:71105/1‑149 (MQ=255)                                                          attGACCACGGCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCAGCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGGTTCgg                                                                                        <  5:71105/147‑1 (MQ=255)                                                             gACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGt                                                                                     >  3:240863/1‑149 (MQ=255)                                                                cACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGgc                                                                                  >  5:207010/1‑149 (MQ=255)                                                                   gTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCa                                                                               <  6:51405/149‑1 (MQ=255)                                                                      gAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAggg                                                                            <  3:140704/149‑1 (MQ=255)                                                                              tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGGGGCAAC‑CATCGCTTCGGCGTgg                                                                                   >  2:224333/1‑134 (MQ=255)                                                                              tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTgg                                                                                   <  1:224333/134‑1 (MQ=255)                                                                              tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCtcgtcg                                                                    >  8:288135/1‑149 (MQ=255)                                                                              tGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCATCGGTCACCGAGACTGCACGCGCACCCGGGCGAGTGGACTCGCCACCCAAATCGATCAACGATGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGAAGGGCCtcgtcg                                                                    >  6:375644/1‑149 (MQ=255)                                                                                  gATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGct                                                                >  6:217754/1‑149 (MQ=255)                                                                                      gCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgac                                                             <  1:209736/148‑1 (MQ=255)                                                                                       ccTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgactg                                                           <  1:329639/149‑1 (MQ=255)                                                                                        cTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGctgactg                                                           >  3:299451/1‑148 (MQ=255)                                                                                           gACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGc                                                       >  1:245527/1‑149 (MQ=255)                                                                                             ccATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGa                                                     >  4:316098/1‑149 (MQ=255)                                                                                              cATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGAc                                                    >  6:320446/1‑149 (MQ=255)                                                                                               aTTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGAcc                                                   <  5:320446/149‑1 (MQ=255)                                                                                               aTTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGAcc                                                   >  2:265579/1‑149 (MQ=255)                                                                                                   gcgcCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCa                                               >  2:293020/1‑149 (MQ=255)                                                                                                     gcCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATc                                             <  4:260718/149‑1 (MQ=255)                                                                                                      ccACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCg                                            <  7:288135/149‑1 (MQ=255)                                                                                                        acacGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCAGCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGga                                          <  3:147678/149‑1 (MQ=255)                                                                                                          acGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGgaga                                        <  2:248849/149‑1 (MQ=255)                                                                                                           cGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGaa                                       <  7:57682/149‑1 (MQ=255)                                                                                                           cGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCAATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGaa                                       <  5:227132/149‑1 (MQ=255)                                                                                                            gCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAg                                      <  3:316098/149‑1 (MQ=255)                                                                                                              tCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTg                                                                                    >  7:348763/1‑101 (MQ=255)                                                                                                              tCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTg                                                                                    <  8:348763/101‑1 (MQ=255)                                                                                                                       tcGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGt                           >  4:29010/1‑149 (MQ=255)                                                                                                                        cGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGt                           >  2:140489/1‑148 (MQ=255)                                                                                                                        cGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGt                           >  3:260164/1‑148 (MQ=255)                                                                                                                             aCCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     >  3:146190/1‑149 (MQ=255)                                                                                                                             aCCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     <  1:293020/149‑1 (MQ=255)                                                                                                                              ccGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     >  2:104398/1‑148 (MQ=255)                                                                                                                               cGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTga                   <  4:240863/149‑1 (MQ=255)                                                                                                                                gagaCTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     >  6:202868/1‑146 (MQ=255)                                                                                                                                    gtccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     >  8:181096/4‑142 (MQ=255)                                                                                                                                    gtccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATAtt                     <  7:181096/139‑1 (MQ=255)                                                                                                                                     tccACGCGCCCGGGTGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATAc             >  4:157739/3‑149 (MQ=255)                                                                                                                                     tGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATAc             >  8:166258/1‑149 (MQ=255)                                                                                                                                      gCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAACGAAACGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCGGACTGAAGCGACCGCCAACGGAGAAGGAATAGGGGGTGATATTGAGAATAc             >  1:295413/1‑149 (MQ=255)                                                                                                                                      gCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATAcc            >  1:233272/1‑149 (MQ=255)                                                                                                                                       cACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGACGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATAccc           <  2:233272/149‑1 (MQ=255)                                                                                                                                       cACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggggg                            <  7:226926/132‑1 (MQ=255)                                                                                                                                       cACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCggggg                            >  8:226926/1‑132 (MQ=255)                                                                                                                                        aagcgcACCCGGGCGGGTGGACTCGGCAACGATATCGATGAACGTTGCGGCGTCGGCAAC‑CATCGCTTAGGCGTGGCGCAGGGCATCGTAGCGCGTACTGAAGCGAACGCCATCGGAGAAGGAATCGGGGGTTATATTGAGAATAccc           <  6:118457/146‑1 (MQ=255)                                                                                                                                        aCGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCa          >  4:83045/1‑149 (MQ=255)                                                                                                                                         cgcgCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCAt         <  3:29010/149‑1 (MQ=255)                                                                                                                                         cgcgCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCAt         <  4:146190/149‑1 (MQ=255)                                                                                                                                         cgcgCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCAt         <  6:247293/149‑1 (MQ=255)                                                                                                                                         cgcgCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCAt         <  5:202868/149‑1 (MQ=255)                                                                                                                                             cACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACa     <  4:220409/149‑1 (MQ=255)                                                                                                                                                ccGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATgg  >  7:327281/1‑149 (MQ=255)                                                                                                                                                ccGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATgg  <  4:299451/149‑1 (MQ=255)                                                                                                                                                ccGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATgg  <  6:357928/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CCCGGCCCCCAGGCGCGCGGCTTCACGCATGACAGCAGGCGTGGAGGTATCGACCGAAATGACCACGTCGAGGCGACTGTTGATGGCCTCGACCATTGGCGCCACACGCTCCAGCTCCTCGGTCACCGAGACTGCACGCGCACCCGGGCGCGTGGACTCGCCACCGATATCGATCAACGTTGCGCCGGCGGCAAC‑CATCGCTTCGGCGTGGCGCAGGGCCTCGTCGCGCTGACTGAAGCGACCGCCATCGGAGAAGGAATCGGGGGTGATATTGAGAATACCCATGACATGG  >  NC_002947/5364277‑5364568

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