Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,834	C→T	6.3%	G10S (GGT→AGT)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,834	0	C	T	6.3%	283.7 / 5.2	95	G10S (GGT→AGT)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (49/39);  new base T (4/2);  total (53/42)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.06e-01

CAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017693‑1017982                                                                                                                                              |                                                                                                                                                      cAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTaaaa                                                                                                                                                  <  1:5472/147‑1 (MQ=255)  aGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTg                                                                                                                                               <  7:218520/149‑1 (MQ=255)   gTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTg                                                                                                                                               >  7:311182/1‑148 (MQ=255)      agcaagcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGTaca                                                                                                                                           >  4:269407/1‑149 (MQ=255)       gcaagcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGacac                                                                                                                                          <  3:252263/149‑1 (MQ=255)         aagcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCa                                                                                                                                        >  6:80500/1‑149 (MQ=255)          agcaGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACTTTTAGGTGGACACCAt                                                                                                                                       <  3:343740/149‑1 (MQ=255)                 aTCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGc                                                                                                                                >  2:137392/1‑149 (MQ=255)                   cAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGCAGGCCATGGTGTCCACCTAAAAGTGGACACCAGTTCTGGc‑c                                                                                                                              >  8:326636/1‑148 (MQ=255)                    aaaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGc‑cc                                                                                                                             <  5:22928/148‑1 (MQ=255)                    aaaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCa                                                                                                                            <  6:357900/149‑1 (MQ=255)                    aaaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCa                                                                                                                            >  1:217840/1‑149 (MQ=255)                     aaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTGAAAGTGGACACCATTTCTGGC‑CCAt                                                                                                                           >  8:100515/1‑149 (MQ=255)                     aaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGACCGGCCAACACATCATCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCAt                                                                                                                           >  5:356663/1‑149 (MQ=255)                        cATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTAACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACACTTTCTGGC‑CCATc                                                                                                                          >  5:97137/1‑147 (MQ=255)                             aaaaTGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTGAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGa                                                                                                                   <  7:100515/149‑1 (MQ=255)                                 tGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGAtttt                                                                                                               >  6:155594/1‑149 (MQ=255)                                    cTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                            <  5:155594/149‑1 (MQ=255)                                    cTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                            <  8:311182/149‑1 (MQ=255)                                    cTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                            <  6:331781/149‑1 (MQ=255)                                       gTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAgg                                                                                                          >  5:154713/1‑148 (MQ=255)                                       gTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAg                                                                                                           <  2:429293/147‑1 (MQ=255)                                                acaAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAgc                                                                                                >  7:346245/1‑149 (MQ=255)                                                acaAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAgc                                                                                                >  7:357965/1‑149 (MQ=255)                                                  aaaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAgcgc                                                                                              >  2:198957/1‑149 (MQ=255)                                                   aaaaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGcc                                                                                             >  8:208972/1‑149 (MQ=255)                                                     aaaaCTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGt                                                                                           <  7:209247/149‑1 (MQ=255)                                                          tGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATTGTGTCCCCCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      >  4:69653/1‑149 (MQ=255)                                                           gACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAg                                                                                                                       <  8:230467/115‑1 (MQ=255)                                                           gACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAg                                                                                                                       >  7:230467/1‑115 (MQ=255)                                                           gACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      >  6:367314/1‑148 (MQ=255)                                                            acacTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      <  7:3774/147‑1 (MQ=255)                                                              acTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACaa                                                                                   >  5:435266/1‑148 (MQ=255)                                                                tCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCATCAGCGCCGTACCTACAATa                                                                                 >  5:224955/1‑148 (MQ=255)                                                                  cccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATaaa                                                                               <  8:346245/148‑1 (MQ=255)                                                                  cccTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAAc                                                                              <  6:154713/149‑1 (MQ=255)                                                                    cTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCt                                                                            <  1:198957/149‑1 (MQ=255)                                                                    cTGATAGTCGAATACCGTTATCAGTCCGGCCAACACATCTTCCTGAAGAAACCGGCGTAGGCAATTGTGTCCAGCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCt                                                                            <  5:141734/149‑1 (MQ=255)                                                                     tGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCtt                                                                           >  2:324841/1‑149 (MQ=255)                                                                      gATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCtt                                                                           >  4:308609/1‑148 (MQ=255)                                                                      gATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCtt                                                                           >  8:392622/1‑148 (MQ=255)                                                                       aTAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCttt                                                                          <  7:358848/148‑1 (MQ=255)                                                                        tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                         <  5:3180/148‑1 (MQ=255)                                                                        tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCa                                                                        <  3:308609/149‑1 (MQ=255)                                                                        tAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACAATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTc                                                                         >  7:400646/1‑148 (MQ=255)                                                                           tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAg                                                                     >  7:141203/1‑149 (MQ=255)                                                                           tCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGAGGACACCATTTCTGGA‑CCATGAAGCGGATTTTCCAGGAATCCTCAGCGCCGTACCTACAATAAACCTTTAAAAg                                                                     >  3:320211/1‑149 (MQ=255)                                                                             gAATCCCGGTATCCGTCCGGCCTACACATATTACTGACGAAACCAGCGAAGGCCATGGTATCAACCAAAAAGAGGACACAATTTCTGGC‑CAATCAAGCGGATTTTACAGGCAACCTCAGCGCCGTAACTACAATAAACATTTCAAAGcc                                                                   >  7:271933/1‑149 (MQ=255)                                                                              aaTCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGcc                                                                   >  7:153663/1‑148 (MQ=255)                                                                              aaTCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGca                                                                   >  1:259369/1‑147 (MQ=255)                                                                               aTCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCa                                                                 <  2:259369/149‑1 (MQ=255)                                                                                  ccGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                            >  4:16094/1‑103 (MQ=255)                                                                                  ccGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                            <  3:16094/103‑1 (MQ=255)                                                                                   cGGTATCCGTCCGGCCAACACGTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGt                                                             >  6:313651/1‑149 (MQ=255)                                                                                   cGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCa                                                                                                  >  1:230818/1‑112 (MQ=255)                                                                                   cGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCa                                                                                                  <  2:230818/112‑1 (MQ=255)                                                                                   cGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGt                                                             <  2:217840/149‑1 (MQ=255)                                                                                      tATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGcg                                                          <  7:392622/149‑1 (MQ=255)                                                                                         ccgtccgGCCAAAACATCTTCCTGACGAAACCGGCGAAGGCCATTGTGTCCAACTAAAAGTGGACACCATTTCTGGC‑CAATCAAGCGGATTTTCCAGGAATCATCAGCGCCGTACCTACAATAAACCGTTCAAAGCCAAGAGTGcgac                                                        >  8:139771/1‑148 (MQ=255)                                                                                           gtccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      >  7:358163/1‑116 (MQ=255)                                                                                           gtccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTa                                                                                      <  8:358163/116‑1 (MQ=255)                                                                                           gtccgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAg                                                     >  7:140795/1‑149 (MQ=255)                                                                                              cgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGccc                                                  <  5:367314/149‑1 (MQ=255)                                                                                              cgGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGcc                                                   <  8:153663/148‑1 (MQ=255)                                                                                                gCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCgg                                                >  4:267186/1‑149 (MQ=255)                                                                                                       caTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTcc                                         <  1:324841/149‑1 (MQ=255)                                                                                                        aTCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCa                                        >  7:356839/1‑149 (MQ=255)                                                                                                          cTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATc                                      >  2:277190/1‑149 (MQ=255)                                                                                                             ccTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATc                                      <  5:313651/146‑1 (MQ=255)                                                                                                               tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATaa                                                                                >  1:133578/1‑102 (MQ=255)                                                                                                               tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATaa                                                                                <  2:133578/102‑1 (MQ=255)                                                                                                               tGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCaaa                                 >  3:117643/1‑149 (MQ=255)                                                                                                                  cGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACAt                              >  5:157530/1‑149 (MQ=255)                                                                                                                   gAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACGATTTCTGGC‑CCATCAAGCGGATGTTCCAGGCATCCTCAGCGGCGTAGCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATc                             <  6:356663/149‑1 (MQ=255)                                                                                                                   gAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATc                             >  3:33966/1‑149 (MQ=255)                                                                                                                   gAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATc                             <  3:267186/149‑1 (MQ=255)                                                                                                                        ctgCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CGATCAAGCGGAATTTCCAGGCATCCTCAGCGCCGTACCTTCAATATACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGTAAACATCGCCtt                        <  8:271933/147‑1 (MQ=255)                                                                                                                            gAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                    <  3:269407/149‑1 (MQ=255)                                                                                                                             ccagaaaTGGTGTCCACTTTTAGGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                    >  5:149096/7‑148 (MQ=255)                                                                                                                             ccagaaaTGGTGTCCACTTTTAGGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGc                    >  3:426930/7‑148 (MQ=255)                                                                                                                              aGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGAGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCca                  >  1:319545/1‑149 (MQ=255)                                                                                                                                 ccATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGTGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGc                           >  2:360540/1‑137 (MQ=255)                                                                                                                                 ccATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGTGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGc                           <  1:360540/137‑1 (MQ=255)                                                                                                                                  cATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGc              <  1:277190/149‑1 (MQ=255)                                                                                                                                   aTGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGcc             >  1:394289/1‑149 (MQ=255)                                                                                                                                   aTGGTGTCCACATAAAAGGGGACACCATTTCTGGC‑ACATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGcc             <  8:400646/149‑1 (MQ=255)                                                                                                                                     ggTGTCCACTTTTAGGTGGACACCATTTCTGGCTCCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCGTGGGCCACAGCCt            >  1:273142/1‑149 (MQ=255)                                                                                                                                     ggTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAAAATCGCATTAGGCAAAAGcca            >  6:303416/1‑147 (MQ=255)                                                                                                                                      gtgtCCACTTTTAGGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTc           >  1:249253/1‑148 (MQ=255)                                                                                                                                      gtgtCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCTCAGCCTc           >  4:411971/1‑148 (MQ=255)                                                                                                                                       tgtCCACTTTTAGGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCaa         <  6:149096/149‑1 (MQ=255)                                                                                                                                        gtCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt        <  1:53667/149‑1 (MQ=255)                                                                                                                                         tCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAAt        <  6:224955/148‑1 (MQ=255)                                                                                                                                         tCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATg       >  3:274050/1‑149 (MQ=255)                                                                                                                                           cACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGACGCAATgg      >  6:303627/1‑148 (MQ=255)                                                                                                                                              cTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAc  >  1:172715/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                      CAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAAC  >  NC_002947/1017693‑1017982

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