Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,189,689	G→T	6.9%	V171V (GTG→GTT)	PP_1937 →	helicase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,189,689	0	G	T	6.9%	172.9 / 2.1	59	V171V (GTG→GTT)	PP_1937	helicase
Reads supporting (aligned to +/- strand):  ref base G (26/27);  new base T (2/2);  total (29/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

AGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAG  >  NC_002947/2189541‑2189827                                                                                                                                                     |                                                                                                                                           aGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTgc                                                                                                                                           >  2:884442/1‑150 (MQ=255)     gATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCtt                                                                                                                                       <  2:329092/150‑1 (MQ=255)             ccGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTg                                                                                                                               <  2:1272113/150‑1 (MQ=255)                 tCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGa                                                                                                                           >  1:439032/1‑150 (MQ=255)                     gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATATTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       <  1:1356128/150‑1 (MQ=255)                     gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       >  2:1356128/1‑150 (MQ=255)                     gTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACCGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAaga                                                                                                                       <  1:489343/150‑1 (MQ=255)                     gTGCGAATTCGTTACCCGACAGGGTAAAGACGGCTACATGCAAAAATGTACTGCTCCGATGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACATACGCTTCACCTCTCGGGAAaga                                                                                                                       >  2:297556/1‑150 (MQ=255)                      tGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCatgctga                                                                                                                               <  1:635310/141‑7 (MQ=255)                      tGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCatgctga                                                                                                                               >  2:635310/1‑135 (MQ=255)                               gTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTttg                                                                                                             <  1:728997/150‑1 (MQ=255)                                ttACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACATGCGCTTCACCTCTGGGGAAAGAGCGCGCTttgt                                                                                                            <  1:346474/150‑1 (MQ=255)                                 tACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCACATAAACTAGAACAGGTTTACCAACTCAGCATAAAGCGAACGTAAGCTTCACCTGTGGGGAAAGCGCGCGCTttgcc                                                                                                           <  2:287789/150‑3 (MQ=255)                                 ctcccGACAGGGTACAGACGGCTACATGCAAAAATGTACCGCTCCGCTGGAGCAGTGACGCTTGGCTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGAGCATCCCCTCTGGGGAAAGAGCGCGCTttgtt                                                                                                           <  1:397765/148‑1 (MQ=255)                                    ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATCCTCAGGCAGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTgag                                                                                                        >  2:397765/1‑150 (MQ=255)                                    ccGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGTGCGCTCTGTTgag                                                                                                        >  2:807805/1‑150 (MQ=255)                                     cGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCa                                                                                                                                     >  2:146447/1‑120 (MQ=255)                                     cGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGCACGTTCGCTTCa                                                                                                                                     <  1:146447/120‑1 (MQ=255)                                      gACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGGAAGAGCGCGCTTTGTTGAGAc                                                                                                      >  1:966467/1‑150 (MQ=255)                                           ggTACAGACGGCTACATGCACAACTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGGAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtc                                                                                                 >  2:381921/1‑150 (MQ=255)                                                     gCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTGCTTACTCAGCATGAAGCTAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGa                                                                                       <  1:884442/150‑1 (MQ=255)                                                             cAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCtttt                                                                               <  1:743631/150‑1 (MQ=255)                                                             cAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCtttt                                                                               >  2:743631/1‑150 (MQ=255)                                                              aaaaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                              >  1:419714/1‑150 (MQ=255)                                                              aaaaaTGTACTGCTCCGCTGGAGCAGTGACGCTAGGGCCCCAAGGCTCAACGCATCGGGGGTTAGATACTCATCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGGGGTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTg                                                                              >  1:387698/1‑150 (MQ=255)                                                                aaaTGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGa                                                                            <  1:644017/150‑1 (MQ=255)                                                                   tGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTTGACCAGGTTTACATACTCAGCATGAATCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCt                                                                         <  1:297556/150‑1 (MQ=255)                                                                    tcaCTGCTCCAGCGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCtt                                                                        >  2:1204210/3‑150 (MQ=255)                                                                      aCTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTgg                                                                      >  2:1098564/1‑150 (MQ=255)                                                                          cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCttt                                                                                                              <  1:702380/106‑1 (MQ=255)                                                                          cTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCttt                                                                                                              >  2:702380/1‑106 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtctc                                                                                               >  1:1330392/1‑120 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAAtctc                                                                                               <  2:1330392/120‑1 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              <  1:1244443/121‑1 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              >  2:95143/1‑121 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              <  1:95143/121‑1 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTcc                                                                                              >  2:1244443/1‑121 (MQ=255)                                                                           tCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGgaga                                                                 >  1:1281310/1‑150 (MQ=255)                                                                                cggAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGc                                                            <  1:1204210/149‑1 (MQ=255)                                                                                    gCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGaa                                                        <  2:439032/150‑1 (MQ=255)                                                                                    gCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGaa                                                        >  2:1343104/1‑150 (MQ=255)                                                                                        tGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACg                                                    <  1:1098564/150‑1 (MQ=255)                                                                                            gCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCttt                                                >  2:520167/1‑150 (MQ=255)                                                                                                   tCCAGCAACTACACTGGACCAGGTTTACTTTCTCAGCATGAAGCGAACGTGCGCTTCATCTCTGGGGAAAGAGCGTGCTTTGTTGAGACCAATCTCCTAGCTTAGTGCTTTTTTATCTGTGGGATAAACGCCGAAGACGCTTTCATTTAg                                         <  2:648853/150‑1 (MQ=255)                                                                                                       gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAc                                                              >  1:704032/1‑125 (MQ=255)                                                                                                       gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAAc                                                              <  2:704032/125‑1 (MQ=255)                                                                                                       gCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAgctgc                                     <  2:1281310/150‑1 (MQ=255)                                                                                                               acTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGagcag                             >  2:1044849/1‑150 (MQ=255)                                                                                                                cTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGagcagc                            <  1:520167/150‑1 (MQ=255)                                                                                                                         gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTtctc                       >  2:620794/1‑146 (MQ=255)                                                                                                                         gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTtctc                       <  1:620794/146‑1 (MQ=255)                                                                                                                         gTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGt                   >  1:277015/1‑150 (MQ=255)                                                                                                                               aTACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCATTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             >  2:908883/1‑150 (MQ=255)                                                                                                                                 aCTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCg                                                          <  1:952492/103‑1 (MQ=255)                                                                                                                                 aCTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCg                                                          >  2:952492/1‑103 (MQ=255)                                                                                                                                  cTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCt          <  1:1044849/150‑1 (MQ=255)                                                                                                                                   tCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTg         >  2:1153983/1‑150 (MQ=255)                                                                                                                                        aTGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             >  2:1239061/1‑141 (MQ=255)                                                                                                                                        aTGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCa             <  1:1239061/141‑1 (MQ=255)                                                                                                                                          gAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAagag  <  1:1343104/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                           AGCTGATTCGCGCCGGTCACGTGCGAATTCGTTACCCGACAGGGTACAGACGGCTACATGCAAAAATGTACTGCTCCGCTGGAGCAGTGACGCTTGGGTCCAGCAACTACACTGGACCAGGTTTACATACTCAGCATGAAGCGAACGTGCGCTTCACCTCTGGGGAAAGAGCGCGCTTTGTTGAGACCAATCTCCTAGCTGAGTGCTTTTGGAGCTTGGGGAGAAACGCCGAAGACGCTTTCATTGAGCTGCTAGAGCAGCTTCTCAAGTTCGTCACCTGGCAAGAG  >  NC_002947/2189541‑2189827

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