Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	346,119	T→A	15.4%	N576I (AAC→ATC) ‡	PP_0287 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	346,119	0	T	A	15.4%	153.5 / 18.6	65	N576I (AAC→ATC) ‡	PP_0287	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base A (5/5);  total (34/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.68e-01

GCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTT  >  NC_002947/345974‑346260                                                                                                                                                   |                                                                                                                                                gCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGCGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGa                                                                                                                                              >  3:123026/1‑149 (MQ=255)   cATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATg                                                                                                                                            >  5:124869/1‑149 (MQ=255)        cGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGGGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGAGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGaa                                                                                                                                       <  6:146551/149‑1 (MQ=255)             caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTg                                                                                                                                  >  6:199107/1‑149 (MQ=255)             caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGCgtgctgg                                                                                                                                   >  6:56328/1‑142 (MQ=255)             caccaGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGCgtgctgg                                                                                                                                   <  5:56328/148‑7 (MQ=255)                caGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGCgtgctg                                                                                                                                    >  4:223869/1‑139 (MQ=255)                caGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGCgtgctg                                                                                                                                    <  3:223869/144‑6 (MQ=255)                  ggCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTg                                                                                                                            >  4:71839/1‑149 (MQ=255)                     tGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCg                                                                                                                         >  4:217574/1‑149 (MQ=255)                          cgcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAg                                                                                                                    <  6:267474/149‑1 (MQ=255)                           gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAgg                                                                                                                   <  5:136197/149‑1 (MQ=255)                           gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGc                                                                                                                                          <  6:3863/127‑1 (MQ=255)                           gcAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCAACGATGGc                                                                                                                                          >  5:3863/1‑127 (MQ=255)                                    ggAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGAc                                                                                                           >  5:146703/1‑148 (MQ=255)                                     gAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCa                                                                                                         <  3:71839/149‑1 (MQ=255)                                       aGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                       <  6:124869/149‑1 (MQ=255)                                            gtgtCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCttt                                                                                                  >  2:312727/1‑149 (MQ=255)                                                  ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGtt                                                                                                                                                <  3:219760/98‑1 (MQ=255)                                                  ttGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGtt                                                                                                                                                >  4:219760/1‑98 (MQ=255)                                                    gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCt                                                                                            <  3:217574/147‑1 (MQ=255)                                                    gCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGt                                                                                          >  5:117525/1‑149 (MQ=255)                                                     ccGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGt                                                                                          <  3:226929/148‑1 (MQ=255)                                                     ccGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGt                                                                                          <  4:123026/148‑1 (MQ=255)                                                       gcgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTAc                                                                                       >  5:146560/1‑149 (MQ=255)                                                        cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGtt                                                                                                                                                >  1:3605/1‑92 (MQ=255)                                                        cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGtt                                                                                                                                                <  2:3605/92‑1 (MQ=255)                                                        cgTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTAcc                                                                                      >  5:34292/1‑149 (MQ=255)                                                             ttCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGAc                                                                                                           <  4:28420/123‑1 (MQ=255)                                                                    ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                       >  1:228457/1‑120 (MQ=255)                                                                    ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGt                                                                                                       <  2:228457/120‑1 (MQ=255)                                                                    ctCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAg                                                                          <  6:34292/149‑1 (MQ=255)                                                                       aGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGAt                                                                        >  1:13619/1‑148 (MQ=255)                                                                        gCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCa                                                                      <  5:199107/149‑1 (MQ=255)                                                                          gTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACt                                                                    <  6:146560/149‑1 (MQ=255)                                                                            cTTTCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                  >  6:11512/1‑149 (MQ=255)                                                                            cTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTgg                                                                  >  5:84945/1‑149 (MQ=255)                                                                             ttGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggt                                                                 >  2:193924/1‑149 (MQ=255)                                                                                cGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcag                                                              <  3:100868/149‑1 (MQ=255)                                                                                 gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                             >  2:2330/1‑149 (MQ=255)                                                                                 gATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTggtcagg                                                             <  6:146703/149‑1 (MQ=255)                                                                                                    ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGccc                                                 >  7:90364/1‑142 (MQ=255)                                                                                                    ttGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGccc                                                 <  8:90364/142‑1 (MQ=255)                                                                                                         acaACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTGACGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACcg‑gcgg                                     >  1:87038/1‑149 (MQ=255)                                                                                                         acaACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                             >  4:70368/1‑109 (MQ=255)                                                                                                         acaACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGt                                                                             <  3:70368/109‑1 (MQ=255)                                                                                                                ctgTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagc                                                       <  2:69052/124‑1 (MQ=255)                                                                                                                ctgTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTcagc                                                       >  1:69052/1‑124 (MQ=255)                                                                                                                ctgTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGtt                             <  6:84945/149‑1 (MQ=255)                                                                                                                   ttCCAGGTTGGCATTGACCAGCACGCGCATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTcc                           >  1:3666/1‑149 (MQ=255)                                                                                                                     ccAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCtt                         <  5:11512/148‑1 (MQ=255)                                                                                                                               aTTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTTACCG‑CCGGCGACAGGTCcttcgcttcgtgca              >  4:46056/1‑145 (MQ=255)                                                                                                                                   aCCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTc                   <  4:224629/140‑1 (MQ=255)                                                                                                                                   aCCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTc                   >  3:224629/1‑140 (MQ=255)                                                                                                                                   aCCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCg                  <  4:47752/141‑1 (MQ=255)                                                                                                                                   aCCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCg                  >  3:47752/1‑141 (MQ=255)                                                                                                                                   aCCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTtgat          >  5:36953/1‑149 (MQ=255)                                                                                                                                    ccAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCGAGCGGCGTCTGTTCCTTGCCTTCGGTCTtgat          <  1:193924/149‑1 (MQ=255)                                                                                                                                     cAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTtgatg         >  8:146937/1‑148 (MQ=255)                                                                                                                                      aGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa       <  1:2330/149‑1 (MQ=255)                                                                                                                                      aGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGGTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGaa       <  3:211192/149‑1 (MQ=255)                                                                                                                                       gCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAg      >  3:54878/1‑149 (MQ=255)                                                                                                                                        ttcGC‑CATCAACGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCtt                         >  2:213611/3‑129 (MQ=255)                                                                                                                                        ttcGC‑CATCAACGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCtt                         <  1:213611/127‑1 (MQ=255)                                                                                                                                        ttcGC‑CATCAACGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGCTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt     >  8:123211/3‑149 (MQ=255)                                                                                                                                        cACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgt     <  7:146937/149‑1 (MQ=255)                                                                                                                                          cGC‑CATCAACGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGCTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgtg    <  7:123211/148‑1 (MQ=255)                                                                                                                                           gC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAgtgt   >  3:211488/1‑148 (MQ=255)                                                                                                                                           gC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGtt  >  8:2130/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                GCCATTGCGTACCACCAGGCTGCCGCGCAGCTCCTGGAAGGGGGTGTCCTTGCCGCGTGGTTCGCTGCTCAGCGTCTTGCGATTGAGCGTGGCGATGGCTTGGCACAACTGCTGTTCCAGGTTGGCATTGACCAGCACGC‑CATCGTTGATGGCGAAGTTTG‑CGCTGCCGTTCAGGGTGTCGACCAGTGCTTTCTGGCTGTTACCGGTGGCTGTCAGATCACTGGTCAGGGTCAGCAGGCCCTTGACCG‑GCGGCGTCTGTTCCTTGCCTTCGGTCTTGATGAAGTGTT  >  NC_002947/345974‑346260

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