Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,474,635	G→T	6.3%	F315L (TTC→TTA)	PP_3964 ←	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,474,635	0	G	T	6.3%	251.4 / 3.3	79	F315L (TTC→TTA)	PP_3964	transposase
Reads supporting (aligned to +/- strand):  ref base G (34/40);  new base T (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

GACCTCGTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTGGGTTT  >  NC_002947/4474488‑4474783                                                                                                                                                    |                                                                                                                                                     gACCTCGTACCAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCTGGGCGTGAGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGaa                                                                                                                                                    <  2:1416986/150‑1 (MQ=255) gACCTCGTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGaa                                                                                                                                                    <  2:1776898/150‑1 (MQ=255) gACCTCGTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGaa                                                                                                                                                    <  1:771808/150‑1 (MQ=255)  aCCTCGTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAg                                                                                                                                                   <  1:1193454/150‑1 (MQ=255)       gTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAgtcgtc                                                                                                                                              <  1:1299121/150‑1 (MQ=255)       gTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAgtcgtc                                                                                                                                              <  1:477569/150‑1 (MQ=255)                gCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCTAAGTCGTCGCAGACCaa                                                                                                                                     >  2:633820/1‑150 (MQ=255)                gCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCTAAGTCGTCGCAGACCaa                                                                                                                                     >  2:633955/1‑150 (MQ=255)                gCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGCCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCaa                                                                                                                                     >  2:999492/1‑150 (MQ=255)                 cTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAAc                                                                                                                                    >  2:1272143/1‑150 (MQ=255)                   gATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCTACCCGTTTAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACtt                                                                                                                                  <  1:407425/150‑1 (MQ=17)                    aTCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTg                                                                                                                                 >  1:361988/1‑150 (MQ=255)                        ctggaacggCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGt                                                                                                                             >  2:454779/8‑150 (MQ=255)                        aGTGCAGGGCCTGTTCAGCCAGCAGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGt                                                                                                                             >  2:944103/1‑150 (MQ=255)                          tGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCTAAGTCGTCGCAGACCAACTTGCCGTTc                                                                                                                           <  1:633820/150‑1 (MQ=255)                          tGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCTAAGTCGTCGCAGACCAACTTGCCGTTc                                                                                                                           <  1:633955/150‑1 (MQ=255)                           gCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAATTTGCCGTTcc                                                                                                                          >  2:1770647/1‑150 (MQ=255)                             acggCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAg                                                                                                                        <  1:454779/148‑1 (MQ=255)                             aGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTcc                                                                                                                          >  2:771008/1‑148 (MQ=255)                             aGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTcc                                                                                                                          <  1:771008/148‑1 (MQ=255)                                     ttCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAg                                                                                                                >  1:1010438/1‑150 (MQ=255)                                      tCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCTCTTACCAACTTGCCGTTCCAGTCACCCAgg                                                                                                               <  1:989009/150‑1 (MQ=255)                                          ccagcTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGt                                                                                                           >  2:1176819/1‑150 (MQ=255)                                           cagcTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGtt                                                                                                          <  1:1658980/150‑1 (MQ=255)                                           cagcTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGtt                                                                                                          <  1:1260053/150‑1 (MQ=255)                                           cagcTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCACCGTACCTTGCCGTTCCAGTCACCCAGGAAGtt                                                                                                          <  1:50801/150‑1 (MQ=18)                                               tGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTgcgc                                                                                                      >  2:1774598/1‑150 (MQ=255)                                                   ttttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTc                                                                                                                           <  1:1455581/125‑1 (MQ=255)                                                   ttttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTc                                                                                                                           >  2:1455581/1‑125 (MQ=255)                                                   ttttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCAt                                                                                                  >  1:436419/1‑150 (MQ=255)                                                   ttttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCAt                                                                                                  <  1:481349/150‑1 (MQ=255)                                                   ttttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTAGTAAAACGGTTCGTCGATTCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCaa                                                                                                  >  2:1483186/1‑149 (MQ=17)                                                     ttGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGt                                                                                                <  1:1828448/150‑1 (MQ=255)                                                       gTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTc                                                                                              <  1:944103/150‑1 (MQ=255)                                                       gTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTc                                                                                              >  2:1866906/1‑150 (MQ=255)                                                           gCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCa                                                                                          <  2:15393/150‑1 (MQ=255)                                                           gCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCa                                                                                          >  1:262971/1‑150 (MQ=255)                                                           gCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCa                                                                                          <  1:242889/150‑1 (MQ=255)                                                           gCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCa                                                                                          >  1:15393/1‑150 (MQ=255)                                                                 tGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACgg                                                                                    >  1:238986/1‑150 (MQ=255)                                                                      aTCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTgg                                                                               >  1:1143660/1‑150 (MQ=255)                                                                      aTCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCATTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTgg                                                                               <  1:1163107/150‑1 (MQ=255)                                                                         aaaGAACTTGCGCCGGGCGTGGGCCATGGAGCCGATTTCAGTGATGCCTTGTTCGAGACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCATTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCt                                                                            >  2:1000143/1‑150 (MQ=255)                                                                         aaaGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCt                                                                            >  1:549187/1‑150 (MQ=255)                                                                                ttGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAACCAGCCTTTGTAGCCAGCGGAGTCTTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCg                                                                     >  1:1673041/1‑150 (MQ=14)                                                                                 tGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGATTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGt                                                                    <  2:660158/150‑1 (MQ=255)                                                                                    gcCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTaaa                                                                 >  1:1159111/1‑150 (MQ=255)                                                                                    gcCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAACACGGCTTGGGCTGAAGTCGTaaa                                                                 <  2:1636495/150‑1 (MQ=255)                                                                                     ccGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAac                                                                <  1:1176819/150‑1 (MQ=255)                                                                                              ggCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCttt                                                       >  2:604909/1‑150 (MQ=255)                                                                                              ggCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCttt                                                       <  2:1310562/150‑1 (MQ=255)                                                                                                  aTGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGc                                                   >  1:1917508/1‑150 (MQ=255)                                                                                                     cAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCt                                                <  2:1143660/150‑1 (MQ=255)                                                                                                         cGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGaaa                                            >  2:1466988/1‑150 (MQ=255)                                                                                                          gATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGaaaa                                           >  2:1281087/1‑150 (MQ=255)                                                                                                           aTTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAAc                                          <  1:350695/150‑1 (MQ=255)                                                                                                                  tGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGcgtcg                                   >  2:1718646/1‑150 (MQ=255)                                                                                                                      gCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCt                               <  2:1010438/150‑1 (MQ=255)                                                                                                                        cTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGt                             >  1:459075/1‑150 (MQ=255)                                                                                                                         ttGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCCTTGAGCCCTGAAAACGGCGTCGTGCTGTa                            >  2:743139/1‑150 (MQ=255)                                                                                                                          tGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTAc                           <  1:604909/150‑1 (MQ=255)                                                                                                                               gAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCa                                                                                       <  2:1445860/85‑1 (MQ=255)                                                                                                                               gAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCa                                                                                       >  1:1445860/1‑85 (MQ=255)                                                                                                                                     ggCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATa                >  2:646870/1‑150 (MQ=255)                                                                                                                                      gCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGt                                                                    >  2:755699/1‑97 (MQ=255)                                                                                                                                      gCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGt                                                                    <  1:755699/97‑1 (MQ=255)                                                                                                                                      gCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAg               <  2:1673041/150‑1 (MQ=255)                                                                                                                                             aGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGt        <  1:999492/150‑1 (MQ=255)                                                                                                                                              gCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCAt                                                                                                  >  2:146800/1‑59 (MQ=255)                                                                                                                                              gCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCAt                                                                                                  <  1:146800/59‑1 (MQ=255)                                                                                                                                              gCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTg       <  2:361988/150‑1 (MQ=255)                                                                                                                                               ccAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTgg      >  1:884051/1‑150 (MQ=255)                                                                                                                                               ccAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCCTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTgg      <  1:743139/150‑1 (MQ=255)                                                                                                                                               ccAGCGAAGTCGTCGCAGACCAACTCGCCGTTCCCGACGTCCGTGAGCTGGCGCGCATGTTCGCCAGCACGGCGTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTgg      <  1:1866906/150‑1 (MQ=17)                                                                                                                                                cAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTggg     <  2:238986/150‑1 (MQ=255)                                                                                                                                                cAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAAACACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTggg     >  1:911994/1‑150 (MQ=255)                                                                                                                                                  gCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCTTGAAGGCGTCGTCGTGCTGTACGCCCAGACATAGGCCCGGTGGGtt   <  1:1774598/150‑1 (MQ=255)                                                                                                                                                  gCGAAGTCGTCGCAGACCAACTTGCCGTTCAAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTGGGtt   <  2:911994/150‑1 (MQ=255)                                                                                                                                                   cGAAGTCGTCGCAGACCAATTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGTTTTGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTGGGttt  <  1:1770647/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GACCTCGTACAAGCCGCTGATCGAGTGCAGGGCCTGTTCAGCCAGCTGACTTTTGTTCGCCACGTGCAAATCAAAGAACTTGCGCCGGGCGTGGGCCATGCAGCCGATTTCAGTGATGCCTTGTTCGAAACCGGCTTTGTAGCCAGCGAAGTCGTCGCAGACCAACTTGCCGTTCCAGTCACCCAGGAAGTTGCGCGCATGTTCGCCAGCACGGCTTGGGCTGAAGTCGTAAACCACCGCTTTGAGCCCTGAAAACGGCGTCGTGCTGTACGCCCAGACATAGGCCCGGTGGGTTT  >  NC_002947/4474488‑4474783

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