Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,379,150	T→C	7.8%	L117P (CTG→CCG) ‡ A63A (GCA→GCG)	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,150	0	T	C	7.8%	180.2 / 5.0	64	L117P (CTG→CCG) ‡ A63A (GCA→GCG)	PP_5639 PP_3855	hypothetical protein hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (31/28);  new base C (3/2);  total (34/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

CCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTC  >  NC_002947/4379005‑4379296                                                                                                                                                   |                                                                                                                                                   ccATAGGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGcc                                                                                                                                                 >  7:340596/1‑149 (MQ=255) ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAGTACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGcc                                                                                                                                                 >  8:191388/1‑149 (MQ=255) ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTTTGcc                                                                                                                                                 >  6:257110/1‑149 (MQ=255) ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                                  <  5:274030/148‑1 (MQ=255) ccATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGc                                                                                                                                                  >  6:274030/1‑148 (MQ=255)    tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCg                                                                                                                                                >  3:351833/1‑147 (MQ=255)    tACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCAGAGCCTGCGCCGGa                                                                                                                                              <  4:191357/149‑1 (MQ=255)           ttCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCg                                                                                                                                       >  3:44658/1‑149 (MQ=255)            tCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGACCCTTCg                                                                                                                                       >  7:123551/1‑148 (MQ=255)                   gACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATc                                                                                                                               <  1:154392/149‑1 (MQ=255)                   gACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATc                                                                                                                               <  7:338000/149‑1 (MQ=255)                      aTGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGac                                                                                                                            <  2:83386/149‑1 (MQ=255)                        ggCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACCCAAGCTATCAGCAATAAAAGCTAGGCAATATGGTGAGACCAAGC‑TTCAGCACTAAGATGCCTCAGGATCCGGCGCAGGCTCTGACGGAGCATTCGCAAGGATCGaca                                                                                                                           >  1:216971/1‑148 (MQ=255)                         gCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACAcc                                                                                                                         >  4:55249/1‑149 (MQ=255)                                   gCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCt                                                                                                                 <  7:274019/147‑1 (MQ=255)                                   gCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑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                                           <  6:107258/149‑1 (MQ=255)                                                                                                       ccAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACgg                                           >  7:258037/1‑149 (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                >  1:43504/1‑149 (MQ=255)                                                                                                                                   gCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTggagcggag                >  7:28478/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                                                                       >  6:182633/1‑85 (MQ=255)                                                                                                                                            cAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCgg                                                                       <  5:182633/85‑1 (MQ=255)                                                                                                                                             aGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgcg                                                    <  3:177782/103‑1 (MQ=255)                                                                                                                                             aGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATcgcg                                                    >  4:177782/1‑103 (MQ=255)                                                                                                                                                 tctGCCGGAGCCTTCGCAATGATCGTCAGCTTAGCGCTGTAATCCTTCAGTACCTTGTCTGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGTCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTc  <  8:123551/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                   CCATACGAGATTCTGATCGACATGGCTGCGTCTGGCGTTGCTGACACAATTGTTCTGCCCTACACCAGCTATCAGCAATACAAGCTAGGCAATATGGTGATACCAAGC‑TTCAGCCCGCAGATGCCTCAGGCTCCGGCGCAGGCTCTGCCGGAGCCTTCGCAATGATCGACACCTTAGCGCTGTAATCCTTCAGTACCTGGTCGGTTAATACCTGCGCAGTCGGGAACTGACTCAGGATCGCGCGGGCACGGGTATCAGCCCCCTCCTGTGGAGCGGAGCGGCTGTTGTTTTC  >  NC_002947/4379005‑4379296

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