Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,432,715	T→C	7.8%	S17P (TCG→CCG)	PP_3049 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,432,715	0	T	C	7.8%	177.3 / 4.7	64	S17P (TCG→CCG)	PP_3049	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (32/27);  new base C (2/3);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

GTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAGGTAGGCCGC  >  NC_002947/3432569‑3432858                                                                                                                                                   |                                                                                                                                                 gTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCg                                                                                                                                                <  5:156321/149‑1 (MQ=255)    gATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTg                                                                                                                                             >  6:133163/1‑149 (MQ=255)          ccGCCTGATGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTgg                                                                                                                                         <  6:158955/147‑1 (MQ=255)            gCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGa                                                                                                                                     >  1:142839/1‑149 (MQ=255)               agagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGAGATCAACCCGCTGc                                                                                                                                            >  6:85477/1‑139 (MQ=255)               agagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGAGATCAACCCGCTGc                                                                                                                                            <  5:85477/139‑1 (MQ=255)                gagGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCg                                                                                                                                 >  6:323289/1‑149 (MQ=255)                   gACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCt                                                                                                                                              <  6:1149/133‑1 (MQ=255)                   gACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCt                                                                                                                                              >  5:1149/1‑133 (MQ=255)                         ggcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAg                                                                                                                                          <  4:168884/131‑1 (MQ=255)                         ggcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAg                                                                                                                                          >  3:168884/1‑131 (MQ=255)                          gcggcACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGccccc                                                                                                                        >  6:243875/1‑148 (MQ=255)                                ccTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑tt                                                                                                                <  4:49009/149‑1 (MQ=255)                                  tACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑tttt                                                                                                              >  7:272612/1‑149 (MQ=255)                                     acCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAg                                                                                                           <  2:142839/149‑1 (MQ=255)                                       cATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGa                                                                                                         <  5:133163/149‑1 (MQ=255)                                        aTCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGAc                                                                                                        >  2:93338/1‑149 (MQ=255)                                              ggCAGTGAGGCTGGTGGGGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCa                                                                                                                               <  4:3692/121‑1 (MQ=255)                                              ggCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCa                                                                                                                               >  3:3692/1‑121 (MQ=255)                                                 aGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTggcc                                                                                                                                       <  6:175410/110‑3 (MQ=38)                                                 aGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCTGCAGCGAGATCAACCCGCTGCTggcc                                                                                                                                       >  5:175410/1‑108 (MQ=38)                                                           gAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTccc                                                                                       <  5:243875/147‑1 (MQ=255)                                                                     ccTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCa                                                                                                    >  1:347014/1‑124 (MQ=255)                                                                     ccTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCa                                                                                                    <  2:347014/124‑1 (MQ=255)                                                                     ccTTGTAACGCTGATCTTGAGGAAGGCGGGATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCa                                                                           <  8:199287/149‑1 (MQ=255)                                                                       ttGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATg                                                                         >  3:74897/1‑149 (MQ=255)                                                                         gTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAg                                                                                                           >  3:94645/1‑113 (MQ=255)                                                                         gTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAg                                                                                                           <  4:94645/113‑1 (MQ=255)                                                                          tAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                                       <  8:272612/148‑1 (MQ=255)                                                                          tAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGACGGCCAGCAGCGCGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTACAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcg                                                                      >  5:282988/1‑149 (MQ=255)                                                                              gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                                    >  6:33403/1‑147 (MQ=255)                                                                              gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                                    >  7:140130/1‑147 (MQ=255)                                                                              gCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCcgcg                                                                    >  8:241861/1‑147 (MQ=255)                                                                                 gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTTGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGc                                                                >  4:143661/1‑148 (MQ=255)                                                                                 gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                                       <  2:154395/141‑1 (MQ=255)                                                                                 gATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                                       >  1:154395/1‑141 (MQ=255)                                                                                     ttGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACgg                                                           >  7:42948/1‑149 (MQ=255)                                                                                       gAGGAAGGCGGTATGAGTGTTGGTAGCTCTCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGAc                                                             <  7:241861/145‑1 (MQ=255)                                                                                       gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                                       >  6:6080/1‑135 (MQ=255)                                                                                       gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACggg                                                          <  5:33403/148‑1 (MQ=255)                                                                                       gAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑GTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGcc                                                                       <  5:6080/135‑1 (MQ=255)                                                                                                     gAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGAGCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                                      <  1:112940/138‑1 (MQ=255)                                                                                                     gAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                                      >  2:112940/1‑138 (MQ=255)                                                                                                       gtgtTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctg                                         >  3:269601/1‑149 (MQ=255)                                                                                                            ggTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                                    >  3:85123/1‑149 (MQ=255)                                                                                                             gTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTgg                                   >  2:133572/1‑149 (MQ=255)                                                                                                             gTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTgg                                   >  3:191486/1‑149 (MQ=255)                                                                                                                gCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCtt                                <  4:85123/149‑1 (MQ=255)                                                                                                                 cTATCACAAGCAGTTGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTc                               <  3:143661/149‑1 (MQ=255)                                                                                                                  tATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGGATTTTCAGGACTTCACGCGGCTGGTCCCCAAAgc                                                                                 <  1:46777/99‑1 (MQ=255)                                                                                                                  tATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAgc                                                                                 >  2:46777/1‑98 (MQ=255)                                                                                                                    tCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                                      <  5:160852/123‑1 (MQ=255)                                                                                                                    tCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGt                                                      >  6:160852/1‑123 (MQ=255)                                                                                                                      acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                                    >  1:322374/1‑139 (MQ=255)                                                                                                                      acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGAcgctgcgctg                                    <  2:322374/139‑1 (MQ=255)                                                                                                                      acaAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAAc                           <  4:269601/148‑1 (MQ=255)                                                                                                                         aGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCca                        >  2:8634/1‑148 (MQ=255)                                                                                                                                ggCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGcc                 <  7:249394/148‑1 (MQ=255)                                                                                                                                ggCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGcc                 >  8:249394/1‑148 (MQ=255)                                                                                                                                    agcagcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGa            <  1:93338/149‑1 (MQ=255)                                                                                                                                    agcagcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGa            <  4:191486/149‑1 (MQ=255)                                                                                                                                        gcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCCAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCGTGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAg           >  3:253962/1‑146 (MQ=255)                                                                                                                                        gcGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCGTGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAg           <  4:253962/146‑1 (MQ=255)                                                                                                                                               gATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAGGTAGgccgc  >  4:32308/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 GTCGATCTTCCGCCAGAGGACGGTGGCGGCACCTACACCATCGCCGGCAGTGAGCCTGGAGGCGATGGCCTTGTAACGCTGATCTTGAGGAAGGCGGTATGAGTGTTGGTAGCTATCACAAGCAGTCGGCCAGCAGCGGGTTGATCTCGCTGCAGGTGGACCCGCAGGCCCCCAAAGG‑TTTTCAGGACTTCACGCGGCTGGTCCCCAAAGCGGCGCATGCCGCGCAGCGACGGGCAGTCAACAAGACGCTGCGCTGGCTTCGAACCCACATCGCCCGGGAGGTAGGCCGC  >  NC_002947/3432569‑3432858

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