Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,383,070	T→G	5.6%	E77D (GAA→GAC)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,383,070	0	T	G	5.6%	257.3 / 4.3	84	E77D (GAA→GAC)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (35/43);  new base G (3/2);  total (38/46)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01

TACTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACAT  >  NC_002947/1382922‑1383216                                                                                                                                                     |                                                                                                                                                    tACTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAt                                                                                                                                                     >  4:137860/1‑149 (MQ=255)  aCTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAtt                                                                                                                                                    >  7:183086/1‑149 (MQ=255)  aCTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAt                                                                                                                                                     <  4:20638/148‑1 (MQ=255)   cTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTc                                                                                                                                                   <  8:183086/149‑1 (MQ=255)   cTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTc                                                                                                                                                   <  3:190967/149‑1 (MQ=255)   cTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTc                                                                                                                                                   >  2:300557/1‑149 (MQ=255)   cTTTGCGGCAGGGTACCCAACGCGAAGACGTTCACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTc                                                                                                                                                   <  8:220358/149‑1 (MQ=255)    tttGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCg                                                                                                                                                  <  8:200726/149‑1 (MQ=255)    tttGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGGAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCACAACGACGGGTCACTTCGAATTCg                                                                                                                                                  >  6:352433/1‑149 (MQ=255)       gCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGt                                                                                                                                                <  5:407882/148‑1 (MQ=255)        cGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGa                                                                                                                                              <  8:48739/149‑1 (MQ=255)         ggCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATTTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGAc                                                                                                                                             >  4:160465/1‑149 (MQ=255)          gCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAGTcaccgaat                                                                                                                                            <  6:259353/149‑8 (MQ=255)           cAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACtt                                                                                                                                           <  8:225047/149‑1 (MQ=255)            aGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTc                                                                                                                                          <  4:17925/149‑1 (MQ=255)                  cccAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAAGTCACCGAATTCGAAgt                                                                                                                                     <  5:397651/148‑2 (MQ=255)                   ccAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAgcgc                                                                                                                                   <  5:444636/149‑1 (MQ=255)                     aaGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCTCGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCtt                                                                                                                                 >  5:309327/1‑149 (MQ=255)                     aaGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCtt                                                                                                                                 <  6:371372/149‑1 (MQ=255)                      aGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATATCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCttt                                                                                                                                <  5:320408/149‑1 (MQ=255)                       ggCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCtttt                                                                                                                               <  2:238561/149‑1 (MQ=255)                          gAAGACGTTGACCACATTGGCGCCGACACCCGCCACCTCGGCCACTGCACCCAGCTGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCt                                                                                                                            <  2:410221/149‑1 (MQ=255)                          gAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAgcgc                                                                                                                                   >  7:391176/1‑142 (MQ=255)                          gAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAgcgc                                                                                                                                   <  8:391176/142‑1 (MQ=255)                           aaGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTc                                                                                                                           <  1:64398/149‑1 (MQ=255)                           aaGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTc                                                                                                                           <  1:303157/149‑1 (MQ=255)                            aGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCa                                                                                                                          >  8:252931/1‑149 (MQ=255)                            aGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGACATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCa                                                                                                                          >  1:145707/1‑149 (MQ=255)                                gTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCa                                                                                                                      >  8:127142/1‑149 (MQ=255)                                gTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCa                                                                                                                      >  4:102869/1‑149 (MQ=255)                                 ttGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCAc                                                                                                                     >  5:45710/1‑149 (MQ=255)                                    aCCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCAcgc                                                                                                                   <  1:221148/148‑1 (MQ=255)                                     ccACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTc                                                                                                                           >  8:15690/1‑139 (MQ=255)                                     ccACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTc                                                                                                                           <  7:15690/139‑1 (MQ=255)                                              cgcCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGc                                                                                                        >  3:6942/1‑149 (MQ=255)                                               gcCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCt                                                                                                       >  1:106029/1‑149 (MQ=255)                                                 cGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGc                                                                                                     <  3:413634/149‑1 (MQ=255)                                                   acacCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCgg                                                                                                   <  7:127142/149‑1 (MQ=255)                                                         gCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTc                                                                                             <  2:106029/149‑1 (MQ=255)                                                          ccaccaCGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGACATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTcc                                                                                            <  2:145707/149‑1 (MQ=255)                                                            accacGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCAt                                                                                          <  7:252931/149‑1 (MQ=255)                                                               acGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCAc                                                                                       >  6:104305/1‑149 (MQ=255)                                                                  gCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCtt                                                                                    <  6:45710/149‑1 (MQ=255)                                                                   ccACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTc                                                                                   <  3:102869/149‑1 (MQ=255)                                                                     aCTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCgg                                                                                 <  3:137860/149‑1 (MQ=255)                                                                       tGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCt                                                                               >  5:386099/1‑149 (MQ=255)                                                                         cACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCg                                                                             >  2:327467/1‑149 (MQ=255)                                                                          aCCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGc                                                                            >  4:241549/1‑149 (MQ=255)                                                                            ccAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGcc                                                                           >  7:245098/1‑148 (MQ=255)                                                                            ccAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCg                                                                          <  7:27815/149‑1 (MQ=255)                                                                             cAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGa                                                                         >  8:263824/1‑149 (MQ=255)                                                                                   gATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGa                                                                  <  2:275493/149‑1 (MQ=255)                                                                                       tcGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTg                                                              >  3:279881/1‑149 (MQ=255)                                                                                       tcGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACATCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTg                                                              >  3:278012/1‑149 (MQ=255)                                                                                        cGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTg                                                              >  2:16264/1‑148 (MQ=255)                                                                                                tcCCGCCAGGGTTGGTAAGGCGTGTAGAGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATg                                                     >  1:62705/1‑149 (MQ=255)                                                                                                  ccGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCt                                                   >  6:81726/1‑149 (MQ=255)                                                                                                   cGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTc                                                  <  8:146896/149‑1 (MQ=255)                                                                                                     ccAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGa                                                <  1:18399/149‑1 (MQ=255)                                                                                                         ggTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCa                                                                                                                      <  6:74655/76‑1 (MQ=255)                                                                                                                  ggCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTc                                   >  4:229153/1‑149 (MQ=255)                                                                                                                   gCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCg                                  >  5:201759/1‑149 (MQ=255)                                                                                                                     gtgtAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGAc                                <  3:241549/149‑1 (MQ=255)                                                                                                                           cGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGctc                          <  7:263824/149‑1 (MQ=255)                                                                                                                                aTAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCAt                     <  7:313165/149‑1 (MQ=255)                                                                                                                                aTAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCAt                     >  1:109517/1‑149 (MQ=255)                                                                                                                                  aaCGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATgg                   >  8:347654/1‑149 (MQ=255)                                                                                                                                     gacgGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCACTCCTCATGGTcg                <  8:245098/149‑1 (MQ=255)                                                                                                                                      aagggTCACTTCGAAATCGGGGACTTCGTAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACATCCGATCTGCAGGGGATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcg                <  1:113968/147‑1 (MQ=255)                                                                                                                                       cgGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgct              <  1:327467/149‑1 (MQ=255)                                                                                                                                        aaagcgcTTCGAAGTCACCGAATTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctg             >  6:88662/7‑149 (MQ=255)                                                                                                                                          gTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgcg           >  2:249061/1‑149 (MQ=255)                                                                                                                                           tCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgcgc          >  4:146846/1‑149 (MQ=255)                                                                                                                                             aCTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgcgctg        <  7:223023/149‑1 (MQ=255)                                                                                                                                             aCTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgcgct         <  2:62705/148‑1 (MQ=255)                                                                                                                                              cTTTGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       >  3:426845/1‑149 (MQ=255)                                                                                                                                              cTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGATCCATCACCTTCGGGTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTcgctgcgcttg       <  4:278012/149‑3 (MQ=255)                                                                                                                                              cTTCGAAGTCACCGAATTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       >  1:107228/1‑149 (MQ=255)                                                                                                                                              cTTCGAAGTCACCGAATTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       >  6:124884/1‑149 (MQ=255)                                                                                                                                               ttCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       >  3:33329/1‑148 (MQ=255)                                                                                                                                                 cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTggg      <  5:104305/147‑1 (MQ=255)                                                                                                                                                 cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       >  2:350442/1‑146 (MQ=255)                                                                                                                                                 cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTgg       <  1:350442/146‑1 (MQ=255)                                                                                                                                                 cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGa     >  6:125370/1‑148 (MQ=255)                                                                                                                                                 cGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGa     <  5:81726/148‑1 (MQ=255)                                                                                                                                                  gAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCAGCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACa   >  4:80075/1‑149 (MQ=255)                                                                                                                                                   aaTTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACAt  <  1:300557/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    TACTTTGCGGCAGGTTACCCAAGGCGAAGACGTTGACCACATTGGCGCCGACACCCGCCACCACGGCCACTGCACCCAGCGGGATCTCGTAGATCTCCCGCCAGGGTTGGTAAGGCGTGTAGCGGTCATAACGACGGGTCACTTCGAATTCGGTGACTTCGAAGCGCTTTTGCTCATGCACGCGCACACGTCGCTGCGGAAGCTCCATCACCTTCGGCTCGCCGACAT‑CGATCTGCAGGGTATGCTCGAGCAGTTTGCGCTCGACGCGCTCCTCATGGTCGCTGCGCTGGGACAT  >  NC_002947/1382922‑1383216

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