Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,685,639	T→A	7.1%	Q709L (CAG→CTG)	PP_2354 ←	two‑component system sensor histidine kinase/response regulator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,685,639	0	T	A	7.1%	207.0 / 3.4	67	Q709L (CAG→CTG)	PP_2354	two‑component system sensor histidine kinase/response regulator
Reads supporting (aligned to +/- strand):  ref base T (27/35);  new base A (3/2);  total (30/37)
Fisher's exact test for biased strand distribution p-value = 6.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTAAGAAACACCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGGTTCGTTCACC  >  NC_002947/2685491‑2685782                                                                                                                                                     |                                                                                                                                                ggTAAGAAACACCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTg                                                                                                                                                >  1:224138/1‑150 (MQ=255)       aaaCACCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCAc                                                                                                                                          <  1:110321/150‑1 (MQ=255)        aaCACCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCAGCCGCAca                                                                                                                                         <  2:1309294/150‑2 (MQ=255)         acacCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCt                                                                                                                                        <  1:475593/150‑1 (MQ=255)            ccaccGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGt                                                                                                                                     <  2:224138/150‑1 (MQ=255)                 ggCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCg                                                                                                                                >  1:561929/1‑150 (MQ=255)                      gCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTg                                                                                                                           <  1:606512/150‑1 (MQ=255)                        ggggATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGc                                                                                                                         <  1:349661/150‑1 (MQ=255)                            aTACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCa                                                                                                                                           >  2:27941/1‑128 (MQ=255)                            aTACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCa                                                                                                                                           <  1:27941/128‑1 (MQ=255)                            aTACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTaa                                                                                                                     <  1:153486/150‑1 (MQ=255)                            aTACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTaa                                                                                                                     <  2:561929/150‑1 (MQ=255)                            aTACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTaa                                                                                                                     <  1:873627/150‑1 (MQ=255)                                cGTTCACGTGCAGCCGGCGCCAGTTGTGCACCATTGAACTCGGGCGTGACAAAATCCGTCTGCAGCAGATAAATTTAGTCAGCCACCAAGGCCAAGTCCCACTCGTAACTCTGTACCAGCCACACCTGGTAACCGTATTGGAGTAacaga                                                                                                                 >  2:6746/1‑148 (MQ=255)                                 gTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACAt                                                                                                                >  1:559635/1‑150 (MQ=255)                                  ttCGCGGGCAGCCAGCGCCAGTTGTGCCCCATTGCAATCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATc                                                                                                               <  2:944724/150‑1 (MQ=255)                                             ccagcGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACGTCGGAGCCGTATTGGCGTAACACATCGCCTAGCATGt                                                                                                    >  2:1404374/1‑150 (MQ=255)                                                       ttGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAg                                                                                          <  2:801708/150‑1 (MQ=255)                                                          tgCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAggtg                                                                                       <  2:559635/150‑1 (MQ=255)                                                            cACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGc                                                                                                                         <  1:559794/114‑1 (MQ=255)                                                            cACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGc                                                                                                                         >  2:559794/1‑114 (MQ=255)                                                             aCCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGgtc                                                                                    <  2:675426/150‑1 (MQ=255)                                                               cATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGgtcgt                                                                                  >  2:425165/1‑150 (MQ=255)                                                                         gggCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACcag                                                                        <  2:785964/150‑1 (MQ=255)                                                                         gggCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACcag                                                                        >  2:492725/1‑150 (MQ=255)                                                                           gCATGGCGATATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACcagca                                                                      >  1:903989/1‑150 (MQ=255)                                                                                 cGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCg                                                                >  1:1264882/1‑150 (MQ=255)                                                                                     aTCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTTGCGTAACACATCTCCTAGCATGTCGCGCACCAGGTGGTCGTCTTCCACCAGCAGTACCGtccc                                                            >  2:654839/1‑147 (MQ=255)                                                                                     aTCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACTAGGTGGTCGTCGTCCACCAGCAGTACCGTGCg                                                            <  1:1404374/150‑1 (MQ=255)                                                                                     aTCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCg                                                            <  2:1264882/150‑1 (MQ=255)                                                                                       cagtcagCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGt                                                          <  2:214259/150‑1 (MQ=255)                                                                                          tcagCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTAATGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTAt                                                       >  1:854075/1‑150 (MQ=255)                                                                                            agcagcagGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTAATGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATgg                                                     <  2:854075/150‑1 (MQ=255)                                                                                             gcagcagGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGgtcgtc                                                                                 <  2:597415/121‑1 (MQ=255)                                                                                             gcagcagGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGgtcgtc                                                                                 >  1:597415/1‑121 (MQ=255)                                                                                                gcagGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGc                                                 >  2:403099/1‑150 (MQ=255)                                                                                                    gTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTc                                             <  1:403099/150‑1 (MQ=255)                                                                                                       aaTGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCt                                          <  1:1375167/150‑1 (MQ=255)                                                                                                        aTGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTg                                         >  2:601526/1‑150 (MQ=255)                                                                                                            cgtcATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTaa                                                                                                                     <  1:836310/70‑1 (MQ=255)                                                                                                            cgtcATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTaa                                                                                                                     >  2:836310/1‑70 (MQ=255)                                                                                                             gtcATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCga                                    <  1:130462/150‑1 (MQ=255)                                                                                                              tcATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATg                                                                                                     <  1:915860/84‑1 (MQ=255)                                                                                                              tcATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATg                                                                                                     >  2:915860/1‑84 (MQ=255)                                                                                                                aTCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAggtg                                                                                       <  2:689448/96‑1 (MQ=255)                                                                                                                aTCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAggtg                                                                                       >  1:689448/1‑96 (MQ=255)                                                                                                                 tCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCGGCTTTCGCTGATCgactga                                >  2:936743/1‑150 (MQ=255)                                                                                                                     gCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACcag                                                                        <  1:1093463/106‑1 (MQ=255)                                                                                                                     gCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACcag                                                                        >  2:1093463/1‑106 (MQ=255)                                                                                                                        aGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCa                         >  1:136029/1‑150 (MQ=255)                                                                                                                             aaTGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCCGAATTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCActgct                    <  1:36789/150‑1 (MQ=255)                                                                                                                             aaTGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCActgct                    <  1:601526/150‑1 (MQ=255)                                                                                                                              aTGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCActgctg                   <  1:741049/150‑1 (MQ=255)                                                                                                                                   tGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCActgct                    >  2:926082/1‑144 (MQ=255)                                                                                                                                   tGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCActgct                    <  1:926082/144‑1 (MQ=255)                                                                                                                                    gCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTCGCCGTCGTCCACCAGCAGTACCATCCAAGTAAGCCGACTTTCGCTGATCGACTGACTGAGCACTGCTGGCACgg             <  1:248265/150‑1 (MQ=255)                                                                                                                                     cTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACggg            >  2:248265/1‑150 (MQ=255)                                                                                                                                        gccgcTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCgactgact                              >  2:1007388/1‑129 (MQ=255)                                                                                                                                        gccgcTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCgactgact                              <  1:1007388/129‑1 (MQ=255)                                                                                                                                        gccgcTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGgttc         >  2:443377/1‑150 (MQ=255)                                                                                                                                        gccgcTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGgttc         <  2:136029/150‑1 (MQ=255)                                                                                                                                          ccaggtgCGCCAGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGcc                                                   <  1:1083669/102‑1 (MQ=255)                                                                                                                                          ccaggtgCGCCAGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGcc                                                   >  2:1083669/5‑106 (MQ=255)                                                                                                                                            cTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGgttcgttc     <  1:443377/150‑1 (MQ=255)                                                                                                                                            aggtgCGCCAGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGgttcgttc     >  1:936377/3‑150 (MQ=255)                                                                                                                                             tgtgCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGGTTCGTTCa    >  2:1331614/1‑150 (MQ=255)                                                                                                                                               tgCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGGACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGGTTCGTTCAcc  >  1:400603/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                GGTAAGAAACACCACCGGCAGGCGGGGATACCGTTCGCGGGCAGCCAGCGCCAGTTGTGCACCATTGAACTCGGGCATGGCGAAATCAGTCAGCAGCAGGTCAATGTCGTCATCCAGCAAGGCCAATGCCTGCTCGCCGCTGTGCGCCTGGCGCACCTGGTAGCCGTATTGGCGTAACACATCGCCTAGCATGTCGCGCACCAGGTGGTCGTCGTCCACCAGCAGTACCGTGCGGTTATGGCCGCTTTCGCTGATCGACTGACTGAGCACTGCTGGCACGGGTTCGTTCACC  >  NC_002947/2685491‑2685782

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