Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,379,146	G→C	7.6%	A116P (GCT→CCT)  P65A (CCT→GCT)	PP_5639 → PP_3855 ←	hypothetical protein hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,379,146	0	G	C	7.6%	192.5 / 6.7	66	A116P (GCT→CCT)  P65A (CCT→GCT)	PP_5639 PP_3855	hypothetical protein hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/29);  new base C (3/2);  total (35/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.75e-01

CAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCT  >  NC_002947/4379000‑4379287                                                                                                                                                     |                                                                                                                                              cAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGc                                                                                                                                              <  3:126618/148‑1 (MQ=255)  aGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGGACACAATTTTTATGCCCTCCACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGATCCGGCGCAGGct                                                                                                                                             >  1:83386/1‑149 (MQ=255)      ccATAGGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGcc                                                                                                                                        >  7:340596/1‑149 (MQ=255)      ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAGTACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGcc                                                                                                                                        >  8:191388/1‑149 (MQ=255)      ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTTTGcc                                                                                                                                        >  6:257110/1‑149 (MQ=255)      ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                         >  6:274030/1‑148 (MQ=255)      ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                         <  5:274030/148‑1 (MQ=255)         tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCg                                                                                                                                       >  3:351833/1‑147 (MQ=255)         tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGa                                                                                                                                     <  4:191357/149‑1 (MQ=255)                ttCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCg                                                                                                                              >  3:44658/1‑149 (MQ=255)                 tCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGACCCTTCg                                                                                                                              >  7:123551/1‑148 (MQ=255)                        gACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATc                                                                                                                      <  1:154392/149‑1 (MQ=255)                        gACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATc                                                                                                                      <  7:338000/149‑1 (MQ=255)                           aTGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGac                                                                                                                   <  2:83386/149‑1 (MQ=255)                             ggCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACCCAAGCTATCAGCAATAAAAGCTAGGCAATATGGTGAGACCAAGC‑TTCAGCACTAAGATGCCTCAGGATCCGGCGCAGGCTCTGACGGAGCATTCGCAAGGATCGaca                                                                                                                  >  1:216971/1‑148 (MQ=255)                              gCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACAcc                                                                                                                >  4:55249/1‑149 (MQ=255)                                        gCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCt                                                                                                        <  7:274019/147‑1 (MQ=255)                                        gCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTg                                                                                                       <  4:351833/148‑1 (MQ=255)                                                  actcaATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                             <  8:339530/146‑1 (MQ=255)                                                  acacAATTTTTCTGCCCTACACCCGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                             <  2:216971/149‑1 (MQ=255)                                                  acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTc                                                                                              <  8:340596/148‑1 (MQ=255)                                                  acacAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCa                                                                                             <  5:64822/149‑1 (MQ=255)                                                      aaTTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑ATCAGCCAGCAGATCCCTCAGCCTCCGGCGCAGGCTCTAATCAACCCAACGAAATGATCGACACCTTACAGATGAAATCCTTCAGTAc                                                                                         >  1:316780/1‑149 (MQ=255)                                                           ttCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGt                                                                                    >  4:156439/1‑149 (MQ=255)                                                             cTGCCCTACACCAGCTATCAGCACTACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTc                                                                                   <  7:191388/148‑1 (MQ=255)                                                                cccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCgg                                                                                 <  5:197943/147‑1 (MQ=255)                                                                cccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                <  3:156439/148‑1 (MQ=255)                                                                cccTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGt                                                                                >  4:366298/1‑148 (MQ=255)                                                                     cacCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTTTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                          <  5:257110/149‑1 (MQ=255)                                                                        cAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGCTTTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCt                                                                        >  1:268828/1‑149 (MQ=255)                                                                                aGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTACGCAATGATCGACACCTAAGCGCTGTAATCCTTCAGTACATGGTCGGTTAATACCTGCGCAGTCg                                                               >  1:362918/1‑149 (MQ=255)                                                                                     tACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCATCAGTACCTGGTCGGTTAATACCTGCGCAGAAGGGaa                                                           >  8:243991/1‑148 (MQ=255)                                                                                              ggCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAgg                                                 >  5:107258/1‑149 (MQ=255)                                                                                                 aaTATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                          >  8:183990/1‑121 (MQ=255)                                                                                                 aaTATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATAc                                                                          <  7:183990/121‑1 (MQ=255)                                                                                                            aCCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACg                                   >  6:311505/1‑149 (MQ=255)                                                                                                             gcAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                  <  3:366298/148‑1 (MQ=255)                                                                                                             ccAAGC‑TTCAGCCCGCATATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATAGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGTCACg                                   <  7:243991/148‑1 (MQ=255)                                                                                                             ccAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                  >  7:258037/1‑149 (MQ=255)                                                                                                             ccAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                  <  6:107258/149‑1 (MQ=255)                                                                                                              cAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTg                                                                                      >  8:223328/1‑96 (MQ=255)                                                                                                              cAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTg                                                                                      <  7:223328/96‑1 (MQ=255)                                                                                                              cAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACggg                                 >  4:23718/1‑149 (MQ=255)                                                                                                                  c‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTAt                              >  2:373491/1‑148 (MQ=255)                                                                                                                  c‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATc                             <  1:373491/149‑1 (MQ=255)                                                                                                                       aGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGcc                         <  2:268828/149‑1 (MQ=255)                                                                                                                        gCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGccc                        >  3:203662/1‑149 (MQ=255)                                                                                                                          ccGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATCAACGACACCTTAGCGCTCTAATCCTTCACTACCTCGTCGGTTAATACCAGCGCAGTCGGGAAATGACTCAGGAGCGCGCGGGCACGGGTATCAGcccc                       >  4:991/1‑148 (MQ=255)                                                                                                                           cGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCcct                     >  7:167649/1‑149 (MQ=255)                                                                                                                                    cgaaGGCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggag            >  7:213463/4‑149 (MQ=255)                                                                                                                                    cTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggag            >  2:361723/1‑149 (MQ=255)                                                                                                                                     gaaGGCTCCGGCAGAGCCTGCGCCGGAGCGGTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGGTCGCGCGGGCACGGGTATCAGAACACACCTGAggagc           >  4:257759/3‑149 (MQ=255)                                                                                                                                       aGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcg          <  4:44658/148‑1 (MQ=255)                                                                                                                                        ggCTCCGGCGCAGGCTCTGCCGGAGGGTTCGCTATGATCGACACCTTAGCGCTGTTATCCTTCAGTACCTGGTCGGTTTATACCTGCGCTGTCGGGTACTGACTCAGGATCGCGGGGGGACGGGTATCAGCCCCCTGCTGTggagcgga        <  3:991/149‑1 (MQ=255)                                                                                                                                        ggCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGGCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcgga        <  5:311505/149‑1 (MQ=255)                                                                                                                                         gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggag       >  7:28478/1‑149 (MQ=255)                                                                                                                                         gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggag       >  1:43504/1‑149 (MQ=255)                                                                                                                                         gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                >  1:135574/1‑124 (MQ=255)                                                                                                                                         gCTCCGGCAGAGCCTGCGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGt                                <  2:135574/124‑1 (MQ=255)                                                                                                                                          cTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggagc      <  4:203662/149‑1 (MQ=255)                                                                                                                                             cGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGc   >  4:90086/1‑149 (MQ=255)                                                                                                                                              ggCGCAGGCTCTGCCGTAGCCTTAGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCt  <  2:362918/149‑1 (MQ=255)                                                                                                                                                  cAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                              <  5:182633/85‑1 (MQ=255)                                                                                                                                                  cAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                              >  6:182633/1‑85 (MQ=255)                                                                                                                                                   aGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgcg                                           >  4:177782/1‑103 (MQ=255)                                                                                                                                                   aGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgcg                                           <  3:177782/103‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                              CAGACCCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCT‑GACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCT  >  NC_002947/4379000‑4379287

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