Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,431,270	T→G	6.2%	L347W (TTG→TGG) ‡	PP_3925 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,431,270	0	T	G	6.2%	253.5 / 7.0	81	L347W (TTG→TGG) ‡	PP_3925	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (43/33);  new base G (2/3);  total (45/36)
Fisher's exact test for biased strand distribution p-value = 6.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.56e-01

ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCAT  >  NC_002947/4431123‑4431414                                                                                                                                                    |                                                                                                                                                 aCAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTg                                                                                                                                                 >  2:13060/1‑149 (MQ=255)   aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                >  3:311178/1‑148 (MQ=255)   aaCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGc                                                                                                                                                <  4:311178/148‑1 (MQ=255)    aCTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAg                                                                                                                                              <  8:194560/149‑1 (MQ=255)      ttCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGt                                                                                                                                            >  8:220993/1‑149 (MQ=255)         gATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCa                                                                                                                                          <  3:170372/148‑1 (MQ=255)          aTTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtt                                                                                                                                        <  1:243805/149‑1 (MQ=255)           ttACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAttt                                                                                                                                       <  8:17066/149‑1 (MQ=255)               tactaTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTGCAAGGTCAtttgttc                                                                                                                                   <  8:61754/149‑5 (MQ=255)                  taTTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                >  7:148974/1‑149 (MQ=255)                    ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTc                                                                                                                                >  1:333626/1‑147 (MQ=255)                    ttGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCa                                                                                                                              <  8:271931/149‑1 (MQ=255)                          tGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCctct                                                                                                                        <  1:261994/149‑1 (MQ=255)                            gCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTg                                                                                                                       <  1:13060/148‑1 (MQ=255)                              cTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAg                                                                                                                     >  6:334167/1‑148 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                               <  2:48022/137‑1 (MQ=255)                               tGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTcc                                                                                                                               >  1:48022/1‑137 (MQ=255)                                    cGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                              >  2:250163/1‑149 (MQ=255)                                     gCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGcc                                                                                                              >  4:258563/1‑148 (MQ=255)                                        tGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCac                                                                                                           >  8:249882/1‑148 (MQ=255)                                          gcaggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacg                                                                                                          <  4:203202/147‑1 (MQ=255)                                             ggcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgac                                                                                                     <  8:106776/149‑1 (MQ=255)                                              gcaTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                    <  1:256759/149‑1 (MQ=255)                                               caTCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCacgacgacg                                                                                                    >  2:205572/1‑148 (MQ=255)                                                   gTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                               <  5:238678/149‑1 (MQ=255)                                                    ttCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATaa                                                                                               >  2:159507/1‑148 (MQ=255)                                                      cagcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACa                                                                                             <  7:220993/148‑1 (MQ=255)                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                             <  3:119931/98‑1 (MQ=255)                                                        gcagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAgg                                                                                                                                             >  4:119931/1‑98 (MQ=255)                                                         cagATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAgg                                                                                          >  7:95287/1‑148 (MQ=255)                                                          agATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCTTTTCGCACAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCAAGCATCTGAGCACAGCCCCCGACGACTGATAACAAGGa                                                                                         >  3:237004/1‑148 (MQ=255)                                                               tactacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAg                                                                                   >  8:296169/1‑149 (MQ=255)                                                                 ctacCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGc                                                                                  <  8:95287/148‑1 (MQ=255)                                                                    ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                     <  3:90965/142‑1 (MQ=255)                                                                    ccGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAAc                                                                                     >  4:90965/1‑142 (MQ=255)                                                                     cGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTc                                                                             <  8:148974/149‑1 (MQ=255)                                                                          ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGAGCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                        <  7:105934/149‑1 (MQ=255)                                                                          ctacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGa                                                                        <  1:250163/149‑1 (MQ=255)                                                                           tacCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAACCAGCATGTCCAAGAc                                                                       >  4:29252/1‑149 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                 <  8:212156/121‑1 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGAcc                                                                      >  6:21001/1‑148 (MQ=255)                                                                             ccATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCAGTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGAt                                                                                                 >  7:212156/1‑121 (MQ=255)                                                                               aTGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCAc                                                                   <  5:71965/149‑1 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                            >  2:207631/1‑88 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGc                                                                                                                            <  1:207631/88‑1 (MQ=255)                                                                                   ccAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGt                                                               >  8:305583/1‑149 (MQ=255)                                                                                    cAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGtt                                                              <  5:21001/149‑1 (MQ=255)                                                                                        cccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                          >  2:268988/1‑149 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAcc                                                         >  8:98853/1‑149 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                          >  8:94097/1‑148 (MQ=255)                                                                                         ccAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                          >  5:228918/1‑148 (MQ=255)                                                                                           aGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTc                                                       >  3:247072/1‑149 (MQ=255)                                                                                             gACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                     <  1:205572/149‑1 (MQ=255)                                                                                                 aGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                 <  1:268988/149‑1 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGCCGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGc                                                 >  6:182000/1‑148 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                          <  4:161219/139‑1 (MQ=255)                                                                                                  gCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGAc                                                          >  3:161219/1‑139 (MQ=255)                                                                                                    gTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCaa                                               <  6:228918/148‑1 (MQ=255)                                                                                                          cccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGcc                                        >  6:320147/1‑149 (MQ=255)                                                                                                           ccAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCt                                       >  3:217661/1‑149 (MQ=255)                                                                                                             aGTTCATTCACATGAACCGCTTGCTTGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCttt                                     <  5:334167/149‑1 (MQ=255)                                                                                                                  attcaCATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTc                                <  1:159507/149‑1 (MQ=255)                                                                                                                   ttgacaTGAACGGCTTGCTGGAACAATTGACCTTGCAGGTCATTTCCAAGTCCTGCCTCTGAGCACAGCCCAGGACGACGGATAACAAGGAAATGAGCATGTGCAAGACGCTCCTGTTCGTCCTCGACGGGAAGATTGCCTTTGTTTcc                               <  6:187312/146‑1 (MQ=255)                                                                                                                        aTGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGgcgc                          >  3:134793/1‑149 (MQ=255)                                                                                                                           aaCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCa                        <  7:249882/148‑1 (MQ=255)                                                                                                                            aCCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGc                      >  8:233403/1‑149 (MQ=255)                                                                                                                              cGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCg                    >  4:33991/1‑149 (MQ=255)                                                                                                                                cTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTg                  <  5:320147/149‑1 (MQ=255)                                                                                                                                   gCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGCGCACAGCCAACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCACGACGGCAAGAATGCCTAAGATTCCGGCGCCAGCCGTGGc                >  7:296205/1‑148 (MQ=255)                                                                                                                                     tGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATc             >  3:195323/1‑149 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg           >  1:143470/5‑149 (MQ=255)                                                                                                                                       ttggaaaTGACCTGCAAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCgg           >  6:255332/5‑149 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                     >  2:16002/1‑107 (MQ=255)                                                                                                                                       gAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGa                                                     <  1:16002/107‑1 (MQ=255)                                                                                                                                           aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg       <  4:223562/149‑1 (MQ=255)                                                                                                                                           aaaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAg       >  4:257488/1‑149 (MQ=255)                                                                                                                                            aaTGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAgg      >  6:32449/1‑149 (MQ=255)                                                                                                                                              tGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGc     >  4:153813/1‑148 (MQ=255)                                                                                                                                               gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGcc    >  8:265609/1‑148 (MQ=255)                                                                                                                                               gACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCa   >  2:55271/1‑149 (MQ=255)                                                                                                                                                aCCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCAt  >  4:17478/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                 ACAACTTCGATTACTACTATTGCTATGGCCTGTTCCGCCTGGCAGGCATCGTTCAGCAGATCTACTACCGCTACTACCATGGCCAGACCCAGGACAAGCGTTTCGCCCAGTTCATTCACATGAACCGCTTGCTGGAACAAATGACCTTGCAGGTCATTTCCAAGTCCAGCCTCTGAGCACAGCCCACGACGACGGATAACAAGGAAAACAGCATGTCCAAGACCCACCTGTTCGACCTCGACGGCAAGATTGCCTTTGTTTCCGGCGCCAGCCGTGGCATCGGCGAGGCCAT  >  NC_002947/4431123‑4431414

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