Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,706,587	T→A	5.2%	L202Q (CTG→CAG)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,706,587	0	T	A	5.2%	240.2 / 2.1	76	L202Q (CTG→CAG)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (37/35);  new base A (2/2);  total (39/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.81e-01

GGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCA  >  NC_002947/5706440‑5706728                                                                                                                                                    |                                                                                                                                              ggCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTGGACGGGGTACACGACAGCCTg                                                                                                                                              >  1:315495/1‑149 (MQ=255) ggCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACTACAGCCTg                                                                                                                                              <  8:165796/149‑1 (MQ=255)  tcTCGACGAAGACAGCCATGTCTCGCGCATTCGCTGGCTGGAGCACACCATGAAGGGTGCTTTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCTATAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                             <  4:22741/148‑1 (MQ=255)  gCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTgg                                                                                                                                             >  2:344213/1‑149 (MQ=255)  gCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAgg                                                                                                                                             >  8:338424/1‑149 (MQ=255)    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCt                                                                                                                                           >  3:241075/1‑149 (MQ=255)    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCt                                                                                                                                           >  3:236442/1‑149 (MQ=255)    tCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCt                                                                                                                                           >  8:3994/1‑149 (MQ=255)        cgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCg                                                                                                                                       <  4:236442/149‑1 (MQ=255)        cgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCg                                                                                                                                       <  4:241075/149‑1 (MQ=255)         gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgt                                                                                                                                      <  7:338424/149‑2 (MQ=255)         gaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCAGGCTGTCgt                                                                                                                                      <  1:96333/149‑2 (MQ=255)          aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                     <  4:80163/149‑1 (MQ=255)           aGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGACCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACCGCCTGGCTGTCGAgg                                                                                                                                    >  4:214669/1‑149 (MQ=255)             aCAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                   >  8:6290/1‑148 (MQ=255)              gacccAAGTCTCGCGGATGCTCTGGCTGGAGCACACCAGCAAGTATGCATTGCTCTAGAACCTACCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGa                                                                                                                                 <  8:326794/146‑1 (MQ=255)                        tcGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCGTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCAc                                                                                                                       >  6:89735/1‑149 (MQ=255)                          gcgcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGt                                                                                                                                         <  4:348246/129‑1 (MQ=255)                          gcgcATGCTCTGGCTGGAGCACACCATAAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGt                                                                                                                                         >  3:348246/1‑129 (MQ=255)                             cATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCg                                                                                                                  >  5:109276/1‑149 (MQ=255)                              aTGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCg                                                                                                                  >  7:305460/1‑148 (MQ=255)                                 ctctGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAAc                                                                                                              <  4:159252/149‑1 (MQ=255)                                    tggctggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACg                                                                                                             <  1:344213/147‑1 (MQ=255)                                        tggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGc                                                                                                       >  5:200249/1‑149 (MQ=255)                                         ggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                      >  2:131327/1‑149 (MQ=255)                                          gAGCAGACCATCAAGGATGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCATATGATGAAAAAGCTGCAATGTACGCCCTGCTCTACGTGGTACACGACAGCATGGCTGGCGAGGCGATACAAGTCACGATCGGAAAGGATGGCgg                                                                                                     >  8:289890/1‑149 (MQ=255)                                          gAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTTCACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgg                                                                                                     <  7:3994/149‑1 (MQ=255)                                          gAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                      <  7:6290/148‑1 (MQ=255)                                           aGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGGTCGGAACGGATGGCggg                                                                                                    <  6:200249/149‑1 (MQ=255)                                            gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                   <  7:107443/149‑1 (MQ=255)                                                acCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcaca                                                                                               >  4:15000/1‑149 (MQ=255)                                                acCATCAAGGCTCCATTGCTCGAGAACCTCCCTCGCAAATGCCTGGCAGATTATGCAAAAGCTGAAATGTACGCCCTGCTCTACGGGGTACACGACAGCCTGGCTGTCGAGTCGATACAAGTCACGCTCGGAACTGATGGCGGGGcaca                                                                                               <  7:289890/149‑1 (MQ=255)                                                      aaGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAACTCACGCTCGGAACGGATGGCGGGGCACACACCGt                                                                                         >  5:88340/1‑149 (MQ=255)                                                        ggCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTa                                                                                       <  5:55431/149‑1 (MQ=255)                                                              aTTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                   <  6:358500/99‑1 (MQ=255)                                                              aTTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGc                                                                                                                                   >  5:358500/1‑99 (MQ=255)                                                                 gCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCg                                                                               >  5:135978/1‑148 (MQ=255)                                                                     gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGAt                                                                            <  6:135978/147‑1 (MQ=255)                                                                       gaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCAc                                                                                                                       >  5:195763/1‑102 (MQ=255)                                                                       gaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCAc                                                                                                                       <  6:195763/102‑1 (MQ=255)                                                                         aCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTg                                                                      <  6:109276/149‑1 (MQ=255)                                                                            tccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                    >  8:121334/1‑148 (MQ=255)                                                                              cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                            <  1:38656/90‑1 (MQ=255)                                                                              cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                            >  2:38656/1‑90 (MQ=255)                                                                              cctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtc                                                                   <  2:82915/147‑1 (MQ=255)                                                                               ctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcatc                                                                <  3:15000/149‑1 (MQ=255)                                                                               ctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcat                                                                 >  2:79568/1‑148 (MQ=255)                                                                                  gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcac                                                              >  7:242887/1‑148 (MQ=255)                                                                                    cAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcaccac                                                           <  1:131327/149‑1 (MQ=255)                                                                                    cAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGACACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                            >  7:338328/1‑148 (MQ=255)                                                                                      aGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGGCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACaa                                                         <  6:248063/149‑1 (MQ=255)                                                                                      aGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACaa                                                         <  8:242887/149‑1 (MQ=255)                                                                                        gCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGt                                                       <  7:121334/149‑1 (MQ=255)                                                                                            ggCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                    >  7:96807/1‑132 (MQ=255)                                                                                            ggCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                    <  8:96807/132‑1 (MQ=255)                                                                                                 atgatgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGTATGGCGGGGcac                                                                                                <  8:171724/99‑1 (MQ=255)                                                                                                 atgatgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGcac                                                                                                >  7:171724/1‑99 (MQ=255)                                                                                                     tgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaacc                                          >  6:104618/1‑149 (MQ=255)                                                                                                     tgCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaacc                                          >  6:65984/1‑149 (MQ=255)                                                                                                        aaaaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGt                                                       >  1:214/1‑133 (MQ=255)                                                                                                        aaaaGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGt                                                       <  2:214/133‑1 (MQ=255)                                                                                                              tGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtg                                  >  2:205062/1‑148 (MQ=255)                                                                                                              tGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtg                                  >  5:220007/1‑148 (MQ=255)                                                                                                                     tACGCCCTTCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgcaa                          <  6:88340/149‑1 (MQ=255)                                                                                                                     tACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgc                            >  7:277971/1‑147 (MQ=255)                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                             <  6:263870/93‑1 (MQ=255)                                                                                                                          ccTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                             >  5:263870/1‑93 (MQ=255)                                                                                                                           cTGCTCGACGGGGTACACGACATCCTGGCTTGCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTTGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTg                    <  8:338328/149‑1 (MQ=255)                                                                                                                               tCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggt                >  8:149051/1‑149 (MQ=255)                                                                                                                                   cGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcgg            >  7:203839/1‑149 (MQ=255)                                                                                                                                      ggTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCggttcggttc         <  5:104618/149‑1 (MQ=255)                                                                                                                                          ctcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtcg     >  8:2176/3‑149 (MQ=255)                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgcaa                          >  5:46419/1‑127 (MQ=255)                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgcaa                          <  6:46419/127‑1 (MQ=255)                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc    >  5:267858/1‑149 (MQ=255)                                                                                                                                             gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCa  <  3:347016/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                              GGCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCA  >  NC_002947/5706440‑5706728

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