Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	202,062	T→A	16.4%	L2523Q (CTG→CAG)	PP_0168 →	putative surface adhesion protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	202,062	0	T	A	16.4%	160.8 / 25.4	67	L2523Q (CTG→CAG)	PP_0168	putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (25/31);  new base A (4/7);  total (29/38)
Fisher's exact test for biased strand distribution p-value = 7.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

TCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAC  >  NC_002947/201915‑202207                                                                                                                                                    |                                                                                                                                                  tCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                 >  2:1325859/1‑150 (MQ=255)  catcACTACCATCACTGACCAGCCAGGAACTACAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGt                                                                                                                                                <  2:1366870/147‑1 (MQ=14)       cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           <  2:1056394/150‑1 (MQ=255)       cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           <  2:831363/150‑1 (MQ=255)       cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           <  1:32598/150‑1 (MQ=255)       cTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           <  2:2751/150‑1 (MQ=255)        tACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGAc                                                                                                                                          >  1:1196644/1‑150 (MQ=255)            gTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGagcca                                                                                                                                       >  1:1140324/1‑145 (MQ=255)            gTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGAgcca                                                                                                                                       <  2:1140324/149‑2 (MQ=255)             gcACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGa                                                                                                                                     <  2:511923/149‑1 (MQ=40)                     gAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           <  2:833896/136‑1 (MQ=255)                     gAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCAGGTCGTGa                                                                                                                                           >  1:833896/1‑136 (MQ=255)                      aGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCCCCATCACGGCCGCGCAGACTTCGGTAGCCGATAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACg                                                                                                                            <  2:1228728/150‑1 (MQ=255)                        ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                               >  2:136883/1‑129 (MQ=39)                        ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                               <  1:136883/129‑1 (MQ=40)                           ggCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                       <  2:1196644/150‑1 (MQ=255)                           ggCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                       >  2:282907/1‑150 (MQ=40)                                tCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAc                                                                                                                  >  2:545481/1‑150 (MQ=40)                                    ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                               <  2:81769/117‑1 (MQ=40)                                    ggCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                               >  1:81769/1‑117 (MQ=40)                                     gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                                             >  2:778464/1‑150 (MQ=255)                                      cAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACTCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                                            <  2:511108/150‑1 (MQ=255)                                      cAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                                            <  1:545481/150‑1 (MQ=255)                                       aaCGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                           <  1:778464/150‑1 (MQ=255)                                             ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGTCCgaa                                                                                                     <  2:955347/150‑3 (MQ=255)                                                  cATCTTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCCGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGAc                                                                                                <  2:341181/150‑1 (MQ=255)                                                    tCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCa                                                                                              <  2:1268102/150‑1 (MQ=255)                                                                ccATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                               >  2:490255/1‑121 (MQ=255)                                                                ccATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                               <  1:490255/121‑1 (MQ=255)                                                                 cATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                           <  1:1233519/92‑1 (MQ=255)                                                                 cATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGa                                                                                                                                           >  2:1233519/1‑92 (MQ=255)                                                                           gACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCCCGAcgcg                                                                       >  1:183031/1‑150 (MQ=255)                                                                            aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACgcgc                                                                      >  2:602886/1‑150 (MQ=255)                                                                                gACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTga                                                                                                   <  1:1259721/117‑1 (MQ=255)                                                                                gACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTga                                                                                                   >  2:1259721/1‑117 (MQ=255)                                                                                   ttCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGa                                                               <  1:282907/150‑1 (MQ=255)                                                                                       gTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                             <  1:76234/84‑1 (MQ=39)                                                                                       gTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                             >  2:76234/1‑84 (MQ=39)                                                                                       gTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGAc                                                           <  2:130777/150‑1 (MQ=255)                                                                                         aGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGt                                                         <  1:1325859/150‑1 (MQ=255)                                                                                            cGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGAt                                                              >  1:417883/1‑142 (MQ=255)                                                                                            cGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGAt                                                              <  2:417883/142‑1 (MQ=255)                                                                                                   atCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                               >  1:1095344/2‑150 (MQ=255)                                                                                                    tCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTa                                                                                    >  2:753141/1‑112 (MQ=255)                                                                                                    tCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTa                                                                                    <  1:753141/112‑1 (MQ=255)                                                                                                       aGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTg                                           >  2:368527/1‑150 (MQ=255)                                                                                                       aGCCAACCTCCCCCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGCCCAGTGCGCCCTACCCCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTg                                           <  1:568806/150‑1 (MQ=14)                                                                                                        ccagaCCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACAccc                                                                               >  2:602204/5‑113 (MQ=255)                                                                                                        ccagaCCTTCACCGTGCACGTCAACCAGCCGCAAGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACAccc                                                                               <  1:602204/109‑1 (MQ=255)                                                                                                          cAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGcc                                        >  1:1172720/1‑150 (MQ=255)                                                                                                            aCCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAg                                      <  2:183031/150‑1 (MQ=255)                                                                                                                    cGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCt                              >  1:1150622/1‑150 (MQ=255)                                                                                                                       gCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTggg                           <  2:1095344/150‑1 (MQ=255)                                                                                                                       gCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTggg                           <  2:778368/150‑1 (MQ=255)                                                                                                                        cACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACTCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGc                          >  2:1020303/1‑150 (MQ=255)                                                                                                                          cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTg                             <  1:1370121/145‑1 (MQ=255)                                                                                                                          cGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTg                             >  2:1370121/1‑145 (MQ=255)                                                                                                                            tCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCa                      <  2:1172720/150‑1 (MQ=255)                                                                                                                                           gtCACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACAACCACGACGCGCAAGGCGATGACGTCTACCAGGACGATGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATg       >  1:803294/2‑150 (MQ=255)                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGcc                                                                                       <  1:236599/69‑1 (MQ=255)                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGcc                                                                                       >  2:236599/1‑69 (MQ=255)                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt      >  1:254235/1‑150 (MQ=255)                                                                                                                                            tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt      >  1:494452/1‑150 (MQ=255)                                                                                                                                            tCACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt      >  1:504383/1‑150 (MQ=255)                                                                                                                                             cACGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTa     <  2:504383/150‑1 (MQ=255)                                                                                                                                                gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAc  <  2:494452/150‑1 (MQ=255)                                                                                                                                                gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTAACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAc  <  1:166078/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  TCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAC  >  NC_002947/201915‑202207

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