Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,030,801	A→T	8.6%	T287S (ACT→TCT)	PP_0887 →	sensor histidine kinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,030,801	0	A	T	8.6%	192.5 / 7.0	69	T287S (ACT→TCT)	PP_0887	sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base A (32/31);  new base T (2/4);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCT  >  NC_002947/1030656‑1030940                                                                                                                                                   |                                                                                                                                            ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACtt                                                                                                                                          >  8:212873/1‑149 (MQ=255) ccAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCACTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACtt                                                                                                                                          >  5:231561/1‑149 (MQ=255)               ccgcaCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTTGCAAAGTCTCTTTGC‑ACAGCtttac                                                                                                                            <  5:63207/149‑5 (MQ=255)                cgccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTGTTGCAAAGTCTCTTTGC‑ACAGCtttacc                                                                                                                           <  5:358229/149‑6 (MQ=255)                 gccCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcg                                                                                                                          >  2:95232/1‑149 (MQ=255)                  cccATGAACTGGGTACGCCACTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgc                                                                                                                         <  6:231561/149‑1 (MQ=255)                   ccATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGTTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgcg                                                                                                                        >  8:217597/1‑149 (MQ=255)                    cATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTgcgcg                                                                                                                        >  7:22728/1‑148 (MQ=255)                     aTGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCa                                                                                                                      <  1:95232/149‑1 (MQ=255)                        aaCTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   >  6:231088/1‑149 (MQ=255)                              gTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGc                                                                                                              >  8:267295/1‑148 (MQ=255)                                         gCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCt                                                                                                  <  5:231088/149‑1 (MQ=255)                                         gCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCt                                                                                                  >  8:163940/1‑149 (MQ=255)                                          cGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTg                                                                                                 <  3:332940/149‑1 (MQ=255)                                           gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                >  7:106160/1‑149 (MQ=255)                                           gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                <  7:212873/149‑1 (MQ=255)                                           gACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTgg                                                                                                >  2:284885/1‑149 (MQ=255)                                            aCCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGc                                                                                               >  4:300498/1‑149 (MQ=255)                                              cATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCACCAGCTGGTGCGCGCAGCCGAAGCCCACCGGCGCCTGGc                                                                                               >  8:87182/1‑147 (MQ=255)                                                tGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATc                                                                                           >  8:91723/1‑149 (MQ=255)                                                   gCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTg                                                                                        <  4:56692/149‑1 (MQ=255)                                                      tGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAg                                                                                     <  7:267295/149‑1 (MQ=255)                                                            tCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGCAGCAGGAt                                                                               <  8:22728/149‑1 (MQ=255)                                                                     tGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      >  6:332296/1‑149 (MQ=255)                                                                           aGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                >  4:340348/1‑149 (MQ=255)                                                                           aGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                <  7:217597/149‑1 (MQ=255)                                                                              aCCACACCGCCCCGCTGTTGCTGGAGGACGTGCAAATCCTCCAGGACCTGGTATCGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCATCCGCCGCCTTGCGATCGTGGAGGAGGATGGGACCGCCTGGCTGGAc                                                             <  4:27313/149‑1 (MQ=255)                                                                               ccACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                             >  6:165269/1‑148 (MQ=255)                                                                                cacaCCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGAc                                                             >  6:33500/1‑147 (MQ=255)                                                                                  cacCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAgg                                                         <  7:91723/149‑1 (MQ=255)                                                                                     cGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGgcgc                                                      <  1:284885/149‑1 (MQ=255)                                                                                       aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      >  6:280413/1‑131 (MQ=255)                                                                                       aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGcc                                                                      <  5:280413/131‑1 (MQ=255)                                                                                       aCCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCt                                                     >  5:249194/1‑148 (MQ=255)                                                                                           gCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAAcc                                                >  4:292192/1‑149 (MQ=255)                                                                                               ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   >  1:169998/1‑78 (MQ=255)                                                                                               ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGcc                                                                                                                   <  2:169998/78‑1 (MQ=255)                                                                                               ttGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             <  4:265719/148‑1 (MQ=255)                                                                                                 gCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGc                                          >  7:176249/1‑149 (MQ=255)                                                                                                  cAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCa                                         >  3:120926/1‑149 (MQ=255)                                                                                                       ggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGa                                    <  5:33500/149‑1 (MQ=255)                                                                                                       ggaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGa                                    >  5:165851/1‑149 (MQ=255)                                                                                                        gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                >  1:321840/1‑120 (MQ=255)                                                                                                        gaCCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCctggctg                                                                <  2:321840/120‑1 (MQ=255)                                                                                                          ccTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgc                                 >  3:192963/1‑149 (MQ=255)                                                                                                             gCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGcc                              >  7:66404/1‑149 (MQ=255)                                                                                                              cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGt                                                                                         >  1:128790/1‑89 (MQ=255)                                                                                                              cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGt                                                                                         <  2:128790/89‑1 (MQ=255)                                                                                                              cAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATgcgc                               <  5:165269/147‑1 (MQ=255)                                                                                                                aaTCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCgg                           <  7:190797/149‑1 (MQ=255)                                                                                                                  tcctccAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGaa                         <  6:165851/149‑1 (MQ=255)                                                                                                                       cAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCa                     <  6:249194/148‑1 (MQ=255)                                                                                                                          gACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa                 <  8:66404/149‑1 (MQ=255)                                                                                                                          gACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTa                 >  4:129939/1‑149 (MQ=255)                                                                                                                             cAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCg              <  8:106160/149‑1 (MQ=255)                                                                                                                              aGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGa                                                           <  4:123713/103‑1 (MQ=255)                                                                                                                              aGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGa                                                           >  3:123713/1‑103 (MQ=255)                                                                                                                              aGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             >  8:64974/1‑117 (MQ=255)                                                                                                                              aGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCt                                             <  7:64974/117‑1 (MQ=255)                                                                                                                                 tAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTc          <  7:163940/149‑1 (MQ=255)                                                                                                                                   caGCTGTTGCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                   <  4:104829/121‑1 (MQ=255)                                                                                                                                   caGCTGTTGCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGAGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGAt                                   >  3:104829/2‑122 (MQ=255)                                                                                                                                   aaGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCa        <  4:192963/149‑1 (MQ=255)                                                                                                                                     gctgttgCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCt            <  6:167885/138‑1 (MQ=255)                                                                                                                                     gctgttgCAAAGTCTCTTTGC‑ACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCt            >  5:167885/6‑143 (MQ=255)                                                                                                                                     gCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgc      <  3:300498/149‑1 (MQ=255)                                                                                                                                      cTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAgcg     >  5:39962/1‑149 (MQ=255)                                                                                                                                          tgCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCt  <  2:164280/149‑1 (MQ=255)                                                                                                                                          tgCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCt  <  6:39962/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                            CCAGGCCGCCGGTGCCGCCCATGAACTGGGTACGCCGCTGGCGACCATGAGCGTGCTGATCAACGAAATGCGCCAGGACCACACCGACCCGCTGTTGCAGGAGGACCTGCAAATCCTCCAGGACCAGGTAAAGCTG‑TGCAAAGAGACTTTGCAACAGCTGGTGCGCGCAGCCGAAGCCAACCGGCGCCTGGCCATCGTGGAGCAGGATGTGACCGCCTGGCTGGACGAGGCGCTCAACCGCTGGCACCTGATGCGCCCGGAAGCCAGCTACCGCTTCCAGCGCCT  >  NC_002947/1030656‑1030940

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