Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,647,897	C→A	6.4%	D430Y (GAC→TAC) ‡	PP_4112 ←	sulfate transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,647,897	0	C	A	6.4%	248.7 / 4.8	78	D430Y (GAC→TAC) ‡	PP_4112	sulfate transporter
Reads supporting (aligned to +/- strand):  ref base C (35/38);  new base A (2/3);  total (37/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

AGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGA  >  NC_002947/4647752‑4648037                                                                                                                                                  |                                                                                                                                             aGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATTCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCg                                                                                                                                          >  2:366275/1‑150 (MQ=255)       ttGCACCGCTCCGACAAATGCACCACGCGCACATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGGACCCGGcccag                                                                                                                                    <  2:112992/150‑4 (MQ=255)       ttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgtg                                                                                                                                    >  1:989170/1‑150 (MQ=255)       ttGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGGACCCGGcccag                                                                                                                                    <  1:844405/150‑4 (MQ=255)        tacACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGAAGTGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCGGTGTGACCTGGGCCGGGGACCGGGcccaag                                                                                                                                   <  2:317541/148‑5 (MQ=255)         gCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGGACCCGGcccaggt                                                                                                                                  >  2:1023221/1‑145 (MQ=255)           aCCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGa                                                                                                                                <  2:50407/150‑1 (MQ=255)                tCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                           <  1:1248101/150‑1 (MQ=255)                tCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                           <  1:538113/150‑1 (MQ=255)                   gACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                           >  2:186128/1‑147 (MQ=255)                   gACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTgg                                                                                                                           <  1:186128/147‑1 (MQ=255)                    aCAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGa                                                                                                                       <  2:973529/150‑1 (MQ=255)                        aTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGaa                                                                                                                   >  1:1355887/1‑150 (MQ=255)                        aTGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGaa                                                                                                                   <  1:639336/150‑1 (MQ=255)                           caccacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGa                                                                                                                >  1:76366/1‑150 (MQ=255)                            accacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                               <  2:919185/150‑1 (MQ=255)                            accacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                               <  2:153633/150‑1 (MQ=255)                            accacGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                               <  1:366275/150‑1 (MQ=255)                                 gcgcAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAgggggt                                                                                                          >  1:766581/1‑150 (MQ=255)                                   gcAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCg                                                                                                        <  2:76366/150‑1 (MQ=255)                                     aaaTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgtg                                                                                                                                    <  2:42900/120‑1 (MQ=255)                                     aaaTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgtg                                                                                                                                    >  1:42900/1‑120 (MQ=255)                                       aTGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGTTGGTCGAGGc                                                                                                    >  2:400721/1‑150 (MQ=255)                                         gCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaa                                                                                                  >  2:962957/1‑150 (MQ=255)                                          cTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaaa                                                                                                 <  1:390874/150‑1 (MQ=255)                                          cTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCaaa                                                                                                 >  2:527256/1‑150 (MQ=255)                                               ccGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACa                                                                                            <  2:134575/150‑1 (MQ=255)                                                 gCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                               >  1:1132685/1‑129 (MQ=255)                                                 gCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAg                                                                                                               <  2:1132685/129‑1 (MQ=255)                                                 gCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                          <  2:1139701/150‑1 (MQ=255)                                                   cTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                        >  1:1012911/1‑150 (MQ=255)                                                          tAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGc                                                                                 >  1:589791/1‑150 (MQ=255)                                                          tAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGc                                                                                 >  1:160545/1‑150 (MQ=255)                                                             cgcTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGt                                                                              >  1:186483/1‑150 (MQ=255)                                                               cTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTg                                                                            <  1:1290263/150‑1 (MQ=255)                                                                     cAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAATCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAaca                                                                      <  1:522258/150‑1 (MQ=255)                                                                     cAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                          >  1:45991/1‑130 (MQ=255)                                                                     cAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGc                                                                                          <  2:45991/130‑1 (MQ=255)                                                                        ggTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgt                                                                                                                                     <  2:705085/84‑1 (MQ=255)                                                                        ggTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCgt                                                                                                                                     >  1:705085/1‑84 (MQ=255)                                                                        ggTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                      <  2:464163/115‑1 (MQ=255)                                                                        ggTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAg                                                                                                      >  1:464163/1‑115 (MQ=255)                                                                          ttttCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgc                                                                 <  2:836386/150‑1 (MQ=255)                                                                              cAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTgccgcc                                                             <  2:1228497/150‑1 (MQ=255)                                                                                gcgcCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCt                                                           >  1:574273/1‑150 (MQ=255)                                                                                  gccgcAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTcc                                                         <  1:962957/150‑1 (MQ=255)                                                                                        aCGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTcc                                                   <  2:160545/150‑1 (MQ=255)                                                                                          ggCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAGCTGGTTGAGGAACGAGGTGGTCGAg                                                                                                      <  2:1387882/97‑1 (MQ=255)                                                                                          ggCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAGCTGGTTGAGGAACGAGGTGGTCGAg                                                                                                      >  1:1387882/1‑97 (MQ=255)                                                                                          ggCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                        >  2:568550/1‑111 (MQ=255)                                                                                          ggCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGt                                                                                        <  1:568550/111‑1 (MQ=255)                                                                                          ggCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGACCGGGTCCGCGGCCGTGTCGGACTGGTTGACGAACGAGGTGGTCGAGACAAAGAACAGCGTGCCAGGCAGTTGCTAACACTTGCCGGACCCCGCGTCATc                                                 >  2:79674/1‑150 (MQ=255)                                                                                           gCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCa                                                >  2:460746/1‑150 (MQ=255)                                                                                             cGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                  <  2:1174431/66‑1 (MQ=255)                                                                                             cGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTc                                                                                                                                  >  1:1174431/1‑66 (MQ=255)                                                                                             cGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGa                                                                                                                <  2:1182601/84‑1 (MQ=255)                                                                                             cGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGa                                                                                                                >  1:1182601/1‑84 (MQ=255)                                                                                             cGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATg                                              >  1:255836/1‑150 (MQ=255)                                                                                                      gACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGACCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGGGTGACAAGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGc                                     >  1:1156275/1‑150 (MQ=255)                                                                                                       aCAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                     >  2:423955/1‑85 (MQ=255)                                                                                                       aCAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAgg                                                                                                     <  1:423955/85‑1 (MQ=255)                                                                                                       aCAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCa                                    >  2:159312/1‑150 (MQ=255)                                                                                                          aaGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATAc                                 <  2:788373/150‑1 (MQ=255)                                                                                                           aGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACa                                <  2:574273/150‑1 (MQ=255)                                                                                                              tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                             <  1:63748/134‑1 (MQ=255)                                                                                                              tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATgg                                             >  2:63748/1‑134 (MQ=255)                                                                                                              tCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGCTGTCGGCATACAGCt                             >  2:1348369/1‑150 (MQ=255)                                                                                                                            tCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCAtgctgc               >  1:732892/1‑150 (MQ=255)                                                                                                                              cAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAAATCTTTCCACTCACTGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCa             >  2:749236/1‑150 (MQ=255)                                                                                                                               aGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg            <  1:527256/150‑1 (MQ=255)                                                                                                                               aGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg            <  2:148786/150‑1 (MQ=255)                                                                                                                               aGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAg            <  2:186483/150‑1 (MQ=255)                                                                                                                               aGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGCTGCAGGCAGACAGCTCCCGCGCATGCTGCCAg            <  1:749236/150‑1 (MQ=255)                                                                                                                                  gTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCa         >  1:618658/1‑150 (MQ=255)                                                                                                                                     gccTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAg      >  2:1010940/2‑150 (MQ=255)                                                                                                                                        acggCCGGGGACCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTg   >  2:1213808/3‑150 (MQ=255)                                                                                                                                         gggCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGa  <  1:1010940/150‑1 (MQ=255)                                                                                                                                         gggCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTAAAGCTATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGa  <  1:337394/150‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                             AGCTTCTTGCACCGCTCCGACAAATGCACCACGCGCAAATGCTTGCCCGCCTTGGCATAGCGCTCGCGCAGGGTTTTCAGCGCCGCAACGGCCGAATAATCGACAAAGCTCAAATGCTGGCAATCCAGCGTGACCTGGGCCGGGTCCCCGGCCGTGTCGAACTGGTTGAGGAACGAGGTGGTCGAGGCAAAGAACAGCGTGCCATGCAGTTGATAACACTTGCCGCCCTCCGCGTCCTCATGGGTGTCGGCATACAGCTCCCGCGCATGCTGCCAGGCAAAGTTGA  >  NC_002947/4647752‑4648037

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