Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,284,555	A→G	7.7%	S165G (AGC→GGC)	PP_1121 →	OmpA family protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,284,555	0	A	G	7.7%	182.4 / 4.7	65	S165G (AGC→GGC)	PP_1121	OmpA family protein
Reads supporting (aligned to +/- strand):  ref base A (27/33);  new base G (2/3);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

GAAGTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAAC  >  NC_002947/1284411‑1284698                                                                                                                                                    |                                                                                                                                                  gaagTCATTGTCATCCGCGATGTGCTCTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGt                                                                                                                                               <  4:189711/149‑1 (MQ=255) gaagTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGt                                                                                                                                               >  1:268072/1‑149 (MQ=255)  aagTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTc                                                                                                                                              <  3:40237/149‑1 (MQ=255)   agTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCTCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTc                                                                                                                                              <  3:59611/148‑1 (MQ=255)      cATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGt                                                                                                                                           >  8:231923/1‑148 (MQ=255)       aTTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTc                                                                                                                                         <  1:11618/149‑1 (MQ=255)       aTTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTc                                                                                                                                         >  7:143565/1‑149 (MQ=255)       aTTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTc                                                                                                                                         <  2:23909/149‑1 (MQ=255)        ttGTCATCCGCGATGTGCACTTCTAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTcc                                                                                                                                        <  4:183398/149‑1 (MQ=255)          gTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CGCTGTCGGtatgg                                                                                                                                       <  3:174137/148‑5 (MQ=255)           tCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑a                                                                                                                                    <  4:296645/148‑1 (MQ=255)             aTCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CGCTGTCGGtatggccg                                                                                                                                    <  2:267725/148‑8 (MQ=255)                cgcgATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTa                                                                                                                              >  3:229861/1‑149 (MQ=255)                    aTGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAAc                                                                                                                          <  2:268072/149‑1 (MQ=255)                    aTGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAAc                                                                                                                          <  4:270524/149‑1 (MQ=255)                       tgCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑Ac                                                                                                                                   >  2:222660/1‑137 (MQ=255)                       tgCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑Ac                                                                                                                                   <  1:222660/137‑1 (MQ=255)                          aCTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGc                                                                                                                                                 >  6:43203/1‑122 (MQ=255)                          aCTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGc                                                                                                                                                 <  5:43203/122‑1 (MQ=255)                          aCTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGaaa                                                                                                                    <  5:54269/149‑1 (MQ=255)                          aCTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGaaa                                                                                                                    <  7:231923/149‑1 (MQ=255)                              cGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACt                                                                                                                  <  5:251005/147‑1 (MQ=255)                                aGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CGCTGTCGGtatgg                                                                                                                                       <  5:111257/126‑5 (MQ=255)                                aGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CGCTGTCGGtatgg                                                                                                                                       >  6:111257/1‑122 (MQ=255)                                  ttCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGa                                                                                                            >  8:87227/1‑149 (MQ=255)                                    cGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAgc                                                                                                          >  5:119618/1‑149 (MQ=255)                                     gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑GCCGCGACAAGCAGAAACTGTCCGAgcg                                                                                                         >  6:217987/1‑149 (MQ=255)                                     gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAgcg                                                                                                         >  1:19073/1‑149 (MQ=255)                                     gATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAgcg                                                                                                         >  3:180345/1‑149 (MQ=255)                                      aTTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAgcgc                                                                                                        >  4:15181/1‑149 (MQ=255)                                                 cTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTc                                                                                             >  5:286066/1‑149 (MQ=255)                                                   gACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCgg                                                                                           >  8:161227/1‑149 (MQ=255)                                                   gACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCgg                                                                                           >  4:240284/1‑149 (MQ=255)                                                    aCGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCgg                                                                                           <  3:225432/148‑1 (MQ=255)                                                      ggCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGa                                                                                        <  4:229861/149‑1 (MQ=255)                                                       gCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGAc                                                                                       <  6:119618/149‑1 (MQ=255)                                                        ccAATGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGAcc                                                                                      >  3:139330/1‑149 (MQ=255)                                                                aaaGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCAGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCt                                                                              >  4:183928/1‑149 (MQ=255)                                                                 aaGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTg                                                                             <  5:217987/149‑1 (MQ=255)                                                                   gagCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGt                                                                           <  8:35629/149‑1 (MQ=255)                                                                        ccTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCgaga                                                                      <  7:87227/149‑1 (MQ=255)                                                                         cTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGAAGGGGGCGTCGGTTCCGCTGCAGCTACAACCAGAAACTGTCCGAGCGCCGTGCACATTCGGTGACCGATTACCTGGTc                                                                          >  1:265133/1‑149 (MQ=255)                                                                             aTACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGgt                                                                 <  8:143565/149‑1 (MQ=255)                                                                                 cATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTAc                                                             >  6:194806/1‑149 (MQ=255)                                                                                   ttGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACcg                                                           >  2:146348/1‑149 (MQ=255)                                                                                      ccACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACcgcg                                                         <  4:180345/148‑1 (MQ=255)                                                                                         cgcgCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCgcagc                                                     <  3:240284/149‑1 (MQ=255)                                                                                         cgcgCCTGAAGCAGGAAGCGCCCTCTGCCCGCCCTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTTTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCgcagc                                                     <  1:318898/149‑1 (MQ=255)                                                                                               tGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTc                                               >  7:212514/1‑149 (MQ=255)                                                                                                 aaGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTc                                               <  6:286066/147‑1 (MQ=255)                                                                                                   gCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGttt                                           >  6:197795/1‑149 (MQ=255)                                                                                                     aGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCg                                         >  3:94091/1‑149 (MQ=255)                                                                                                     aGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCg                                         <  3:183928/149‑1 (MQ=255)                                                                                                        aaGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCggtg                                      <  2:19073/149‑1 (MQ=255)                                                                                                                      cGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGa                        <  3:15181/149‑1 (MQ=255)                                                                                                                        ccTTAGCGGCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGaaa                      <  2:119615/149‑1 (MQ=255)                                                                                                                                  aGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGt                          >  4:177313/1‑135 (MQ=255)                                                                                                                                  aGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGt                          <  3:177313/135‑1 (MQ=255)                                                                                                                                        aaTACCGA‑‑‑CAGCGTCGTTTCCG‑‑ACATCTACACCCAGAACCTGTCCGACCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTAACGCGCAGCAGCTGCGTTTCGGTGGGCGGAGCGGGTGAAACCCAGACGGTGGCAGa      <  3:183399/148‑1 (MQ=255)                                                                                                                                         aTACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGcgcg                             >  2:285377/1‑125 (MQ=255)                                                                                                                                         aTACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGcgcg                             <  1:285377/125‑1 (MQ=255)                                                                                                                                         aTACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGa      >  4:52279/1‑148 (MQ=255)                                                                                                                                            ccGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAAc  >  6:150604/1‑149 (MQ=255)                                                                                                                                            ccGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAAc  <  5:150604/149‑1 (MQ=255)                                                                                                                                             cGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCg                                  >  8:144258/1‑116 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GAAGTCATTGTCATCCGCGATGTGCACTTCGAGTTCGATTCTGCGCGCCTGACGGCCAGTGACAAAGAGCGCCTCAATACCATTGCCACGCGCCTGAAGCAGGAAGCGCCCTCTGCCCGCCTTAGCGTCAGCGGCCATACCGA‑‑‑CAGCGTCGGTTCCG‑‑ACAGCTACAACCAGAAACTGTCCGAGCGCCGTGCCCATTCGGTGACCGATTACCTGGTCGAGAGCGGTGTACCGCGCAGCAGCTTCGTTTCGGTGGTCGGCGCGGGTGAAACCCAGCCGGTGGCAGACAAC  >  NC_002947/1284411‑1284698

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