Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,017,839	A→G	6.3%	L8P (CTT→CCT)	PP_5461 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,017,839	0	A	G	6.3%	274.8 / 6.2	95	L8P (CTT→CCT)	PP_5461	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (50/39);  new base G (4/2);  total (54/41)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.15e-01

CAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCTC  >  NC_002947/1017693‑1017985                                                                                                                                                   |                                                                                                                                                    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                                                                                                                               >  1:217840/1‑149 (MQ=255)                    aaaGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCa                                                                                                                               <  6:357900/149‑1 (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                                                                                                               <  6:331781/149‑1 (MQ=255)                                    cTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCa                                                                                                               <  8:311182/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                                                                                                                          >  7:230467/1‑115 (MQ=255)                                                           gACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAg                                                                                                                          <  8:230467/115‑1 (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)                                                                                                                                                 aaaGTGGACACCATTTCTGGC‑GCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCTc  >  7:13304/1‑149 (MQ=255)                                                                                                                                                 aaaGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCTc  >  7:309894/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGTCAGCAAGCAGGTATCAAAGCATTGAAAATGCCTTGTCCCGTACACAAAAAAACTGACACTCCCCTGATAGTCGAATCCCGGTATCCGTCCGGCCAACACATCTTCCTGACGAAACCGGCGAAGGCCATGGTGTCCACCTAAAAGTGGACACCATTTCTGGC‑CCATCAAGCGGATTTTCCAGGCATCCTCAGCGCCGTACCTACAATAAACCTTTCAAAGCCCAGAGTGCGACGAGCCCGGCGCCTCCATCGCAAACATCGCCTTGGGCCACAGCCTCAATGGCAACCTC  >  NC_002947/1017693‑1017985

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