Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,426,462	C→G	6.6%	Y32* (TAC→TAG)	PP_3920 →	putative Phage integrase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,426,462	0	C	G	6.6%	251.0 / 6.6	76	Y32* (TAC→TAG)	PP_3920	putative Phage integrase
Reads supporting (aligned to +/- strand):  ref base C (31/40);  new base G (3/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.66e-01

CCCACCACCAACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCT  >  NC_002947/4426317‑4426599                                                                                                                                                  |                                                                                                                                           cccACCACCAACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGtacta                                                                                                                                          >  3:167885/1‑148 (MQ=255)  ccaccaccaACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGtacta                                                                                                                                          >  1:350291/1‑147 (MQ=255)  ccaccaccaACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGtacta                                                                                                                                          >  3:220714/1‑147 (MQ=255)    accaccaccGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAgg                                                                                                                                      <  8:137408/149‑1 (MQ=255)       accaACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTAGTACAgggcc                                                                                                                                   >  3:30266/1‑147 (MQ=255)       accaACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGAc                                                                                                                                   >  7:162670/1‑149 (MQ=255)        ccaACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTAGTACAgggcct                                                                                                                                  <  4:30266/149‑4 (MQ=255)         caACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACa                                                                                                                                        <  2:46060/142‑1 (MQ=255)         caACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACa                                                                                                                                        >  1:46060/1‑142 (MQ=255)           aCGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑t                                                                                                                              >  2:291294/1‑149 (MQ=255)              cccgccTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCgg                                                                                                                           >  7:333580/1‑149 (MQ=255)                  ccTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATGGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCTAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCaa                                                                                                                       <  2:166065/149‑1 (MQ=255)                   cTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAg                                                                                                                      <  6:295854/149‑1 (MQ=255)                         cccTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAgg                                                                                                                                      >  8:144593/1‑128 (MQ=255)                         cccTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAgg                                                                                                                                      <  7:144593/128‑1 (MQ=255)                         cccTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATc                                                                                                                <  4:167885/149‑1 (MQ=255)                           cTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCtt                                                                                                              >  5:286477/1‑149 (MQ=255)                            tGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTc                                                                                                             >  1:289109/1‑149 (MQ=255)                              aTTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAACAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTc                                                                                                             <  1:235899/147‑1 (MQ=255)                                 tttCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCg                                                                                                         <  5:286254/148‑1 (MQ=255)                                  ttCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATATATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCgg                                                                                                        >  3:55140/1‑148 (MQ=255)                                     ccGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAAt                                                                                                    <  4:220714/149‑1 (MQ=255)                                         gtgaagtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGAcc                                                                                                >  7:37177/1‑149 (MQ=255)                                          tgaagtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCggg                                                                                                                          >  2:155259/1‑122 (MQ=255)                                          tgaagtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCggg                                                                                                                          <  1:155259/122‑1 (MQ=255)                                            aagtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGGCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTgg                                                                                             >  1:303293/1‑149 (MQ=255)                                            aagtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTg                                                                                              <  6:185707/148‑1 (MQ=255)                                              gtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGtacta                                                                                                                                          <  7:196591/103‑1 (MQ=255)                                              gtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGtacta                                                                                                                                          >  8:196591/1‑103 (MQ=255)                                              gtgaaATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATAACTTCCTCGGCAATGACCTGGt                                                                                            >  3:156856/1‑148 (MQ=255)                                                 aaaTGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAAACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGa                                                                                          >  1:106819/1‑147 (MQ=255)                                                   aTGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGc                                                                                      >  8:2276/1‑149 (MQ=255)                                                      cgcCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTc                                                                                                          >  7:344988/1‑126 (MQ=255)                                                      cgcCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTc                                                                                                          <  8:344988/126‑1 (MQ=255)                                                      cgcCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCaa                                                                                   <  8:333580/149‑1 (MQ=255)                                                      cgcCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCa                                                                                    >  4:56140/1‑148 (MQ=255)                                                         ccAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACt                                                                                >  2:22356/1‑149 (MQ=255)                                                            aGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGaa                                                                             <  1:22356/149‑1 (MQ=255)                                                                   ccGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGaa                                                                             <  6:121728/142‑1 (MQ=255)                                                                   ccGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGaa                                                                             >  5:121728/1‑142 (MQ=255)                                                                     gagaATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGa                                                                                                                    >  8:68991/1‑101 (MQ=255)                                                                     gagaATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGa                                                                                                                    <  7:68991/101‑1 (MQ=255)                                                                          tCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTgg                                                               >  3:83699/1‑149 (MQ=255)                                                                          tCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTg                                                                >  3:337223/1‑148 (MQ=255)                                                                            gcgACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGa                                                              <  4:337223/148‑1 (MQ=255)                                                                            gcgACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAg                                                             <  1:324594/149‑1 (MQ=255)                                                                            gcgACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAg                                                             <  2:289109/149‑1 (MQ=255)                                                                            gcgACCTGCCACCAGGCATGTATCGCCGCAAGAGAAAAAAAGCCAATGGCAAGGGGTGGGAGGCCCTGTACTACAGGGACAGG‑TCTGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAg                                                             <  7:269610/149‑1 (MQ=255)                                                                              gACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAg                                                                            >  5:287447/1‑132 (MQ=255)                                                                              gACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAg                                                                            <  6:287447/132‑1 (MQ=255)                                                                              gACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGc                                                           <  7:2276/149‑1 (MQ=255)                                                                               aCCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGcc                                                          <  3:56140/149‑1 (MQ=255)                                                                                 cTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGcc                                                          >  5:71602/1‑147 (MQ=255)                                                                                      accaGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCgg                                                   >  7:321/1‑149 (MQ=255)                                                                                       ccaGGCATATATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGt                                                  <  4:55140/149‑1 (MQ=255)                                                                                         aGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGc                                                <  6:71602/149‑1 (MQ=255)                                                                                             aTGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGaa                                            <  2:350291/149‑1 (MQ=255)                                                                                                         aaGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGGTTCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCa                                 >  1:101235/1‑149 (MQ=255)                                                                                                            agaCAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGAt                             <  8:162670/149‑1 (MQ=255)                                                                                                               cAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGaa                          <  6:286477/149‑1 (MQ=255)                                                                                                               cAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGaa                          <  8:37177/149‑1 (MQ=255)                                                                                                                aaaaaGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCTAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAg                         <  1:291294/149‑1 (MQ=255)                                                                                                                 aaaaGCCAATGGCAAGGTGTGGGAGGCCCTTTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAgg                        <  3:133765/149‑1 (MQ=255)                                                                                                                 aaaaGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAgg                        <  4:156856/149‑1 (MQ=255)                                                                                                                  aaaGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAgg                        >  6:222076/1‑148 (MQ=255)                                                                                                                   aaGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAggg                       >  3:317762/1‑148 (MQ=255)                                                                                                                      ccAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATc                   <  2:106819/149‑1 (MQ=255)                                                                                                                         aTGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTc                <  8:321/149‑1 (MQ=255)                                                                                                                                  tgtgGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTAc       <  2:101235/149‑1 (MQ=255)                                                                                                                                    tgGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACct     >  5:328202/1‑149 (MQ=255)                                                                                                                                      ggAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAgg                        <  2:229355/128‑1 (MQ=255)                                                                                                                                      ggAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACctgc   <  5:222076/149‑1 (MQ=255)                                                                                                                                       ctgtcccTGTAGTACAGGGACAGG‑TTGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACctgct  >  6:379796/5‑149 (MQ=255)                                                                                                                                        tgtcccTGTAGTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAg                                                             <  2:160335/86‑1 (MQ=255)                                                                                                                                        tgtcccTGTAGTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAg                                                             >  1:160335/4‑89 (MQ=255)                                                                                                                                             cTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGaa                                                                             >  7:151835/1‑68 (MQ=255)                                                                                                                                             cTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGaa                                                                             <  8:151835/68‑1 (MQ=255)                                                                                                                                              tGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTAc       >  3:155663/1‑137 (MQ=255)                                                                                                                                              tGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTAc       <  4:155663/137‑1 (MQ=255)                                                                                                                                               gTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCtt                 <  3:34895/126‑1 (MQ=255)                                                                                                                                               gTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCtt                 >  4:34895/1‑126 (MQ=255)                                                                                                                                                  |                                                                                                                                           CCCACCACCAACGCCCGCCTGGACCCCTGATTTTTCCCGAGTGAAGTGAAATGCGCCCAAGGAAGACCGAGAATCGCGACCTGCCACCAGGCATGTATCGCCGCAAGAGACAAAAAGCCAATGGCAAGGTGTGGGAGGCCCTGTACTACAGGGACAGG‑TCGGGCAAGGATATCTTCCTCGGCAATGACCTGGTGAAAGCCAAACTGAAGTGGGCGGAACTGGAGGCCAAGTCGGTGCCGAAGGAACTGACCACGATGAAGGGAATCTTCGACGAGTACCTGCT  >  NC_002947/4426317‑4426599

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑