Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,416,520	T→C	6.6%	I297T (ATC→ACC)	PP_1238 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,416,520	0	T	C	6.6%	171.6 / 2.7	60	I297T (ATC→ACC)	PP_1238	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base C (2/2);  total (30/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 = 9.79e-01

AGCCTGAGCGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAG  >  NC_002947/1416375‑1416665                                                                                                                                                  |                                                                                                                                                    aGCCTGAGCGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCa                                                                                                                                                   <  1:13896/148‑1 (MQ=255)  gCCTGAGCGAAGACGGCAGCGGCAACCCTGTGCTGTACCTTGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAAc                                                                                                                                                 <  2:277549/149‑1 (MQ=255)  gCCTGAGCGAAGACGGCAGCGGCAACCCTGTGCGGGACCTTGGCTCCGACCTGGCCCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGAGACTTGCGGGGTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAAc                                                                                                                                                 <  2:273508/149‑1 (MQ=255)   ccTGAGCGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAAcc                                                                                                                                                >  7:238176/1‑149 (MQ=255)    cTGAGCGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAAcca                                                                                                                                               <  5:191604/149‑2 (MQ=255)         cGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTcc                                                                                                                                           >  8:223865/1‑148 (MQ=255)           aaGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGa                                                                                                                                         >  1:86669/1‑148 (MQ=255)                ggCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑Gcc                                                                                                                                 <  3:206793/149‑1 (MQ=255)                         aaCCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGGGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTgacgac                                                                                                                          >  7:139400/1‑147 (MQ=255)                          aCCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAg                                                                                                                       <  7:88410/149‑1 (MQ=255)                           cccGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGc                                                                                                                      <  2:129347/149‑1 (MQ=255)                            ccGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAg                                                                                                                       >  7:100086/1‑147 (MQ=255)                             cGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAg                                                                                                                       >  7:188625/1‑146 (MQ=255)                                  cTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAAc                                                                                                                 >  2:110747/1‑147 (MQ=255)                                     tACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGcc                                                                                                            >  4:100604/1‑149 (MQ=255)                                            gCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGcttc                                                                                                     <  4:94589/149‑1 (MQ=255)                                            gCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGcttc                                                                                                     <  5:84915/149‑1 (MQ=255)                                             cTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGcttct                                                                                                    <  6:202171/149‑1 (MQ=255)                                                cGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTg                                                                                                                               >  6:193226/1‑119 (MQ=255)                                                cGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTg                                                                                                                               <  5:193226/119‑1 (MQ=255)                                                    cTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCa                                                                                                 <  6:338951/145‑1 (MQ=255)                                                    cTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCa                                                                                                 >  5:338951/1‑145 (MQ=255)                                                     tGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGc                                                                                                                      >  4:177953/1‑123 (MQ=255)                                                     tGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGc                                                                                                                      <  3:177953/123‑1 (MQ=255)                                                         ccGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGt                                                                                        <  8:188625/149‑1 (MQ=255)                                                          cgcgCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑Gc                                                                                                                                  >  7:258084/1‑106 (MQ=255)                                                          cgcgCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑Gc                                                                                                                                  <  8:258084/106‑1 (MQ=255)                                                              ccTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGc                                                                                   <  7:223865/149‑1 (MQ=255)                                                               cTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCa                                                                                  <  8:100086/149‑1 (MQ=255)                                                                        cGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGa                                                                         <  3:54618/149‑1 (MQ=255)                                                                        cGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGa                                                                         <  8:238176/149‑1 (MQ=255)                                                                         gTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAg                                                                                 >  4:301147/1‑140 (MQ=255)                                                                         gTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAg                                                                                 <  3:301147/140‑1 (MQ=255)                                                                            ggCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCaa                                                                     <  8:139400/149‑1 (MQ=255)                                                                                    ccTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGc                                                             >  5:293979/1‑149 (MQ=255)                                                                                    ccTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGc                                                             >  8:244562/1‑149 (MQ=255)                                                                                                tGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTg                                                 >  8:226262/1‑149 (MQ=255)                                                                                                    gAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGa                                             <  3:100604/149‑1 (MQ=255)                                                                                                     aaTGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAg                                            <  2:86669/149‑1 (MQ=255)                                                                                                              ttttgtaCATAAACCGAAGCTTGGGCTTGTAAAAAATCAACTTTTCCTAACACCCCCTGACGACTAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTAcc                                     <  1:110747/143‑1 (MQ=255)                                                                                                              ttGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAg                                   >  4:210410/1‑149 (MQ=255)                                                                                                              atGAGGACATCAACCGCAGCCTGGGCCTCTACTACATCAACCTCTCCGAAAA‑‑GCCTGACGACAAGCATAACCAGCAGGGCTTCTTCATCCCCATGTTCGGCAGCGAGCCGACCAAGCAAGAGCCTGAAGCACGTGCCGAGCGTAACCAg                                   >  2:249177/2‑149 (MQ=255)                                                                                                                 aGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGtt                                >  1:149827/1‑149 (MQ=255)                                                                                                                         aaCCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGa                         <  7:244562/148‑1 (MQ=255)                                                                                                                           ccGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGAGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGc                       >  4:209047/1‑148 (MQ=255)                                                                                                                             gCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGACCGTTACCAGGTTCGCCTGAGCAAg                    >  3:271035/1‑149 (MQ=255)                                                                                                                              cAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAgg                   >  2:188184/1‑149 (MQ=255)                                                                                                                                 tcTGGGCCTGTACTAAACCAAACTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACAAGCCGGGCTTCTTCATCCGCCTGTTAGGCAGCGAGCCGACCAATGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgg                <  4:88810/148‑1 (MQ=255)                                                                                                                                 ccTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgg                >  6:162883/1‑149 (MQ=255)                                                                                                                                 ccTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgg                >  3:120998/1‑149 (MQ=255)                                                                                                                                   tGGGCCTGTACTCCACCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgg                <  4:93504/147‑1 (MQ=255)                                                                                                                                   tGGGCCTGTACTACACCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGAAAggtggga              >  3:93504/1‑148 (MQ=255)                                                                                                                                   tGGGCCTGTACTACACAAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACAAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGACCTTGACGAGCGTTACCAGGTTCGCCTGAGCAAgatgt                >  3:88810/1‑146 (MQ=255)                                                                                                                                     ggCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTga            <  3:209047/149‑1 (MQ=255)                                                                                                                                          gTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCg       >  7:45255/1‑149 (MQ=255)                                                                                                                                              tacATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgggt              >  6:93488/1‑138 (MQ=255)                                                                                                                                              tacATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgggt              <  5:93488/138‑1 (MQ=255)                                                                                                                                               acATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAg  >  4:194314/1‑149 (MQ=255)                                                                                                                                                 aTCTACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgggt              <  6:243024/135‑1 (MQ=255)                                                                                                                                                 aTCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAggtgggt              >  5:243024/1‑135 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AGCCTGAGCGAAGACGGCAGCGGCAACCCGGTGCTGTACCTGGGCTCCGACCTGGACCGCGCCTGGTCTGGCGTGGGCCGTGCCCTGGAGCAGGGTGGCGAATGGCGGGTTGAGGACATCAACCGCAGCCTGGGCCTGTACTACATCAACCTGTCCGAAAA‑‑GCCTGACGACAAGCAGAACCAGCCGGGCTTCTTCAGCCGCCTGTTCGGCAGCGAGCCGACCAAGGAAGAGCGTGAAGCCCGTGCCGAGCGTTACCAGGTTCGCCTGAGCAAGGTGGGTGAGAGCGTGCAG  >  NC_002947/1416375‑1416665

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