Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,382,537	T→A	9.8%	D255V (GAC→GTC)	PP_1204 ←	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,382,537	0	T	A	9.8%	163.4 / 8.2	61	D255V (GAC→GTC)	PP_1204	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (32/23);  new base A (3/3);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.46e-01

GAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAC  >  NC_002947/1382398‑1382683                                                                                                                                              |                                                                                                                                                   gAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGATGACCTCGAAGa                                                                                                                                           <  6:192247/149‑1 (MQ=255)    cTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGATGACCTCGAAGAtga                                                                                                                                        <  4:315877/149‑3 (MQ=255)     tCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAg                                                                                                                                       >  1:226271/1‑149 (MQ=255)      ccTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGaa                                                                                                                                             <  7:148362/142‑1 (MQ=255)      ccTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGaa                                                                                                                                             >  8:148362/1‑142 (MQ=255)       ctgctgCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCt                                                                                                                                     >  2:59590/1‑149 (MQ=255)          ctgcAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCAt                                                                                                                                  >  8:262651/1‑149 (MQ=255)             cAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATgg                                                                                                                                <  8:36452/148‑1 (MQ=255)                    ggATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCATCGAGGTCATCGAAGATCAGCTCATGGGCTTCGag                                                                                                                        <  3:261945/149‑1 (MQ=255)                         gTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagcagc                                                                                                                    <  2:184024/148‑1 (MQ=255)                              gggTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagc                                                                                                                       >  7:180708/1‑140 (MQ=255)                              gggTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGagc                                                                                                                       <  8:180708/140‑1 (MQ=255)                                gTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTGGTGCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTg                                                                                                              <  1:267627/149‑1 (MQ=255)                                gTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCtt                                                                                                                            >  1:256824/1‑133 (MQ=255)                                gTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCtt                                                                                                                            <  2:256824/133‑1 (MQ=255)                                        gATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                               <  6:9355/106‑1 (MQ=255)                                        gATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCg                                                                                                                                               >  5:9355/1‑106 (MQ=255)                                              gCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGa                                                                                              >  5:83733/1‑149 (MQ=255)                                              gCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGa                                                                                              >  2:33955/1‑149 (MQ=255)                                                  tGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCaa                                                                                          >  7:141927/1‑149 (MQ=255)                                                    ttCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACa                                                                                        <  1:33955/149‑1 (MQ=255)                                                         tttCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTg                                                                                   >  6:270445/1‑149 (MQ=255)                                                                 gCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCg                                                                             <  6:269474/147‑1 (MQ=255)                                                                 gCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGcc                                                                           >  3:164305/1‑149 (MQ=255)                                                                   ttCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTc                                                                         >  2:2311/1‑149 (MQ=255)                                                                       tctTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAc                                                                     >  3:194045/1‑149 (MQ=255)                                                                       tctTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTAc                                                                     >  6:126835/1‑149 (MQ=255)                                                                         ttCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCg                                                                   >  3:260937/1‑149 (MQ=255)                                                                         ttCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCg                                                                   >  7:12507/1‑149 (MQ=255)                                                                           cAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCa                                                                 >  4:135181/1‑149 (MQ=255)                                                                             aGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCAta                                                               <  6:83733/149‑2 (MQ=255)                                                                              ggCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGATGACCTCGAAGAtga                                                                                                                                        >  7:12258/1‑73 (MQ=38)                                                                              ggCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGATGACCTCGAAGAtga                                                                                                                                        <  8:12258/75‑3 (MQ=38)                                                                               gCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGt                                                             >  1:329638/1‑149 (MQ=255)                                                                                cccAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGTGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTc                                                            <  1:80914/149‑1 (MQ=255)                                                                                cccAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTc                                                            >  2:34982/1‑149 (MQ=255)                                                                                 ccAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACTGGCTGGCATCGCCTCGTACCGCATTGTcc                                                           >  1:58863/1‑149 (MQ=255)                                                                                      tCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGt                                                       >  5:58938/1‑148 (MQ=255)                                                                                       cGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTc                                                      <  4:164305/148‑1 (MQ=255)                                                                                              gcgcTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACtt                                              >  3:287054/1‑149 (MQ=255)                                                                                               cgcTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtg                                             <  7:262651/149‑1 (MQ=255)                                                                                                gcTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgc                                            <  1:2311/149‑1 (MQ=255)                                                                                                gcTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtg                                             >  6:277205/1‑148 (MQ=255)                                                                                                    gACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggt                                        >  8:275354/1‑149 (MQ=255)                                                                                                      cccGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtg                                       >  5:236784/1‑148 (MQ=255)                                                                                                        cGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcag                                    <  5:270445/149‑1 (MQ=255)                                                                                                             cACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGc                               >  1:68941/1‑149 (MQ=255)                                                                                                             cACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGc                               <  2:329638/149‑1 (MQ=255)                                                                                                                ccACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGa                            >  1:221447/1‑149 (MQ=255)                                                                                                                 cACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGAc                           <  6:236784/149‑1 (MQ=255)                                                                                                                   cTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGt                         <  2:221447/149‑1 (MQ=255)                                                                                                                            aCATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   <  7:205966/132‑1 (MQ=255)                                                                                                                            aCATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTtgcaggtgcagg                                   >  8:205966/1‑132 (MQ=255)                                                                                                                            aCATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAg                  >  3:52287/1‑149 (MQ=255)                                                                                                                                 ctgatcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCt             >  1:96908/4‑149 (MQ=255)                                                                                                                                  tcatcTTCGAGGTCATCGAAGATCAGCTCATGGTCTACGAGCAGCTTAATGCGCACGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCCCATTGACCTCGTCAACTACTTGCAGGTGCAGGGTGCCGACGTGGTGTAGGTTCTg            >  6:344360/1‑149 (MQ=255)                                                                                                                                   gatcTTCGATGACCTCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGc           >  8:55535/2‑149 (MQ=255)                                                                                                                                    atcTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCt          <  2:68941/149‑1 (MQ=255)                                                                                                                                         cGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGga      <  7:275354/148‑1 (MQ=255)                                                                                                                                            ggTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAc  <  6:58938/149‑1 (MQ=255)                                                                                                                                            ggTCATCGAAGATCAGCTCATGGGCTTCCAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAc  >  1:191098/1‑149 (MQ=255)                                                                                                                                              |                                                                                                                                                   GAACTCCTGCTGCAGCTCCGGATCGTTGCGGGTCAGCTCGATCAGGCTCTGTTCCATTTCACTGGCTTCCTCTTCAAGGCCCAGTTCGGACAGGCGCTTGACCCGGTGCACCCACTG‑‑GCCGACATCATCATCTTCGAGGTCATCGAAGATCAGCTCATGGGCTTCGAGCAGCTTATTGCGCAGGGTGGCGCTGACCAACAGGCTGGCATCGCCTCGTACCGCATTGTCCTCGTCAAGTACTTGCAGGTGCAGGGTGCCGACGTGGTTGAGGTTCTGCTCGGAGAAC  >  NC_002947/1382398‑1382683

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