Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,297,132	T→A	5.9%	E77D (GAA→GAT)	TpbB ←	diguanylate cyclase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,297,132	0	T	A	5.9%	244.3 / 3.0	78	E77D (GAA→GAT)	TpbB	diguanylate cyclase
Reads supporting (aligned to +/- strand):  ref base T (34/39);  new base A (3/2);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAG  >  NC_002947/5296984‑5297272                                                                                                                                                         |                                                                                                                                             gCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGAt                                                                                                                                              >  1:54192/1‑149 (MQ=255)     aTCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGt                                                                                                                                          >  4:120313/1‑149 (MQ=255)      tCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGtt                                                                                                                                         <  5:15116/149‑1 (MQ=255)      tCCAGCAAAGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGtt                                                                                                                                         >  6:19739/1‑149 (MQ=255)       ccAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGAt                                                                                                                                              <  1:179185/143‑1 (MQ=255)       ccAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGAt                                                                                                                                              >  2:179185/1‑143 (MQ=255)        cAGCAAAGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTgg                                                                                                                                       <  5:19739/149‑1 (MQ=255)         aGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTggc                                                                                                                                      <  3:120313/149‑1 (MQ=255)         aGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTggc                                                                                                                                      >  2:267469/1‑149 (MQ=255)          gCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTggcg                                                                                                                                     >  5:13549/1‑149 (MQ=255)             aGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATACCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCATTggcggcg                                                                                                                                  >  7:60978/1‑149 (MQ=255)                     gCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGt                                                                                                                           <  2:54192/148‑1 (MQ=255)                         ccTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCg                                                                                                                      <  6:13549/149‑1 (MQ=255)                            gCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATc                                                                                                                       >  1:84972/1‑145 (MQ=255)                            gCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATc                                                                                                                       <  2:84972/145‑1 (MQ=255)                               ccTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAACGAATCGCTggcgt                                                                                                                                    >  6:257860/1‑128 (MQ=255)                               ccTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAACGAATCGCTggcgt                                                                                                                                    <  5:257860/129‑2 (MQ=255)                                cTCCAGCTGCGCGAGTATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACGTCTTCGGTACTGGCAATAA‑‑‑ACGCCAGCGAATc                                                                                                                                            >  6:295993/1‑120 (MQ=255)                                        gcgcGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTggcggcg                                                                                                                                  >  1:326116/1‑122 (MQ=255)                                        gcgcGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTggcggcg                                                                                                                                  <  2:326116/122‑1 (MQ=255)                                                     ggTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGt                                                                                          >  8:330714/1‑149 (MQ=255)                                                      gTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGtt                                                                                         <  8:42290/149‑1 (MQ=255)                                                        aTTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGa                                                                                       >  2:28117/1‑149 (MQ=255)                                                         ttGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGa                                                                                       <  7:330714/148‑1 (MQ=255)                                                         ttGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGa                                                                                       <  5:182627/148‑1 (MQ=255)                                                            tCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCACTTACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGa                                                                                       <  1:278085/148‑1 (MQ=255)                                                                gTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCCCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGCAGTTGATCgcccgc                                                                               >  1:277932/1‑144 (MQ=255)                                                                   gCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGc                                                                             >  1:299822/1‑148 (MQ=255)                                                                   gCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGc                                                                             >  6:92318/1‑148 (MQ=255)                                                                    ccAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGAACGAGCGGGCGa                                                                           >  6:289200/1‑149 (MQ=255)                                                                      aGTTCGCCAGTTTCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                         <  7:167985/149‑1 (MQ=255)                                                                          cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAAc                                                                      <  1:267469/148‑1 (MQ=255)                                                                          cgcCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCATTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGCCCGCTCGATCAAc                                                                      <  8:60978/148‑1 (MQ=255)                                                                              aGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGAACGAGCGGGCGATCAACTGCAg                                                                 <  5:289200/149‑1 (MQ=255)                                                                               gTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAgg                                                                <  1:28117/149‑1 (MQ=255)                                                                                ttGCTCGCCCTCGTTGTTGAAGACCTTGGCCCCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGCAGTTGATCGCCCGCTCGATCAACTACaccg                                                               <  2:277932/149‑4 (MQ=255)                                                                                   cTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTa                                                             <  5:92318/148‑1 (MQ=255)                                                                                    tCGCCCTCGTTGTTGAAGACCTTGGCCTCGGACCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTAt                                                            >  1:36592/1‑149 (MQ=255)                                                                                      gCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATgg                                                          <  2:299822/148‑1 (MQ=255)                                                                                       cccTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGtt                                                        >  1:203665/1‑149 (MQ=255)                                                                                         cTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                      >  3:196522/1‑149 (MQ=255)                                                                                         cTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                      >  8:216076/1‑149 (MQ=255)                                                                                         cTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTgg                                                      <  7:216076/149‑1 (MQ=255)                                                                                         cTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGGTATGGTTgg                                                      >  8:249966/1‑149 (MQ=255)                                                                                            gttgttGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGc                                                     <  2:36592/147‑1 (MQ=255)                                                                                               gttgAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTTGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                 <  1:22218/148‑1 (MQ=255)                                                                                                   aaGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTTTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                 <  6:255537/128‑1 (MQ=255)                                                                                                   aaGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAg                                                                 >  5:255537/1‑128 (MQ=255)                                                                                                            gCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAgg                                        <  5:272126/144‑1 (MQ=255)                                                                                                            gCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAgg                                        >  6:272126/1‑144 (MQ=255)                                                                                                               tCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCCGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCAACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                >  6:286467/1‑149 (MQ=255)                                                                                                               tCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                    >  5:236628/1‑145 (MQ=255)                                                                                                               tCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGa                                    <  6:236628/145‑1 (MQ=255)                                                                                                                 gg‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCAc                                <  2:203665/147‑1 (MQ=255)                                                                                                                    ccACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGcc                             >  4:188271/1‑148 (MQ=255)                                                                                                                     cACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCca                           >  4:174669/1‑149 (MQ=255)                                                                                                                      aCTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                         >  4:266296/1‑102 (MQ=255)                                                                                                                      aCTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATc                                                                         <  3:266296/102‑1 (MQ=255)                                                                                                                      aCTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACg                                            >  5:196360/1‑131 (MQ=255)                                                                                                                      aCTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACg                                            <  6:196360/131‑1 (MQ=255)                                                                                                                      aCTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCcag                          <  3:188271/149‑1 (MQ=255)                                                                                                                        ttcttcGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCg                                                                                                                      <  5:233229/55‑1 (MQ=255)                                                                                                                        ttcttcGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCg                                                                                                                      >  6:233229/1‑55 (MQ=255)                                                                                                                           ttcGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATAGAGCGGGCGATCAACTGCAGGTTATGGTTGTCATAGACAAGCAGTGCGAGCACGCCCAGCAggg                     >  4:253964/1‑149 (MQ=255)                                                                                                                              ggTACTGTCAATCA‑‑‑ACGCCAGCGATTCGTTTGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGTCGAGCACGCCCAGCAGGGTGa                  <  3:253964/149‑1 (MQ=255)                                                                                                                               gTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTg                                                       <  8:236533/111‑1 (MQ=255)                                                                                                                                   tGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                <  1:118220/146‑1 (MQ=255)                                                                                                                                   tGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAgg                >  2:118220/1‑146 (MQ=255)                                                                                                                                   tGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                 >  1:185930/1‑145 (MQ=255)                                                                                                                                   tGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAg                 <  2:185930/145‑1 (MQ=255)                                                                                                                                    ggCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                      <  1:211711/123‑1 (MQ=255)                                                                                                                                    ggCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGc                                      >  2:211711/1‑123 (MQ=255)                                                                                                                                     gCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                      >  2:292393/1‑138 (MQ=255)                                                                                                                                     gCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAgg                      <  1:292393/138‑1 (MQ=255)                                                                                                                                     gCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGgaga            <  7:249966/148‑1 (MQ=255)                                                                                                                                               aCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTa                                                                                            <  5:175172/61‑1 (MQ=255)                                                                                                                                               aCGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTa                                                                                            >  6:175172/1‑61 (MQ=255)                                                                                                                                                cGCCAACGAATCGCTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGc     >  8:230584/1‑147 (MQ=255)                                                                                                                                                 gCCAACGAATCGCTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAg  <  7:158909/149‑1 (MQ=255)                                                                                                                                                  ccAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                 >  2:174197/1‑101 (MQ=255)                                                                                                                                                  ccAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAg                                                 <  1:174197/101‑1 (MQ=255)                                                                                                                                                         |                                                                                                                                             GCTCATCCAGCAAGGCGCTGGCAACCTGCACCTCCAGCTGCGCGAGCATGCCGGTATTGTCGCGTTGCCAGTGCGCCAGTTGCTCGCCCTCGTTGTTGAAGACCTTGGCCTCGG‑CCACTTCTTCGGTACTGGCAATCA‑‑‑ACGCCAGCGATTCGTTGGCGGCGGCGCTGTCATCGAACACCACAGCCGCTTCTACGGTGTAGTTGATCGAGCGGGCGATCAACTGCAGGTTATGGTTGGCATAGACACGCAGGGCGAGCACGCCCAGCAGGGTGAGGGAGATACCCGCCAG  >  NC_002947/5296984‑5297272

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