Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,436,997	T→A	5.1%	S188C (AGC→TGC)	PP_4775 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,436,997	0	T	A	5.1%	307.5 / 4.7	98	S188C (AGC→TGC)	PP_4775	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (50/43);  new base A (2/3);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 6.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCCGTGCGGGCGGTA  >  NC_002947/5436849‑5437144                                                                                                                                                     |                                                                                                                                                    cGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCt                                                                                                                                                     >  1:6728/1‑149 (MQ=255)      gCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCt                                                                                                                                                 <  6:97878/148‑1 (MQ=255)        aCTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGc                                                                                                                                               >  8:277750/1‑148 (MQ=255)               gCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTAc                                                                                                                                        <  2:6728/148‑1 (MQ=255)                cGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTg                                                                                                                                      <  7:277750/149‑1 (MQ=255)                 ggTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCAGACAGCCAggaatatg                                                                                                                                     <  3:283424/149‑7 (MQ=255)                  gTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTgg                                                                                                                                     >  1:369961/1‑148 (MQ=255)                    cctgctcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                                                  >  3:22788/1‑149 (MQ=255)                     ctgctcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                                                  >  4:28254/1‑148 (MQ=255)                     ctgctcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                                                  <  2:61706/148‑1 (MQ=255)                       gctcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCCGGCTGTCTGCCAGGTACTGGTAATAg                                                                                                                               >  2:238961/1‑149 (MQ=255)                        ctcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAgg                                                                                                                              <  8:260172/149‑1 (MQ=255)                         tcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAg                                                                                                                               >  8:46451/1‑147 (MQ=255)                         tcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAg                                                                                                                               <  7:46451/147‑1 (MQ=255)                         tcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCAGACAGCCAggaatatg                                                                                                                                     <  7:125303/141‑7 (MQ=255)                         tcctgctGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCAGACAGCCAggaatatg                                                                                                                                     >  8:125303/1‑135 (MQ=255)                           ctgctgTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTAt                                                                                                                           >  2:20131/1‑149 (MQ=255)                           ctgctgTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTAt                                                                                                                           >  5:327645/1‑149 (MQ=255)                              ctgTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATc                                                                                                                          <  4:242405/147‑1 (MQ=255)                                 ttCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGcca                                                                                                                     >  4:163036/1‑149 (MQ=255)                                    cTGTCCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGGTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGccacc                                                                                                                   <  6:423987/148‑1 (MQ=255)                                    cTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCAGACAGCCAggaatat                                                                                                                                      <  6:420953/129‑6 (MQ=255)                                    cTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACAAATTCCTGGCAGACAGCCAggaatat                                                                                                                                      >  5:420953/1‑124 (MQ=255)                                     tGACCCTGTCCATACATCTGATAAGCCACCAGGTCATAGCCGCGATTTGGCGGGGCCAGCGGCAGGGGCGACGGCACGCCGAGCGGCGGGTTGAGCACCACATTTTCCTGGCTGTCTGCCAg                                                                                                                                            <  8:288787/122‑1 (MQ=255)                                     tGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAg                                                                                                                                            >  7:288787/1‑122 (MQ=255)                                     tGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGccac                                                                                                                    >  2:75089/1‑146 (MQ=255)                                       aCCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaa                                                                                                               <  2:238097/149‑1 (MQ=255)                                        cccTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaaa                                                                                                              <  3:163036/149‑1 (MQ=255)                                        cccTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaaa                                                                                                              <  3:153271/149‑1 (MQ=255)                                         ccTGGCCATACATCTGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCaga                                                                                                                >  1:224473/1‑146 (MQ=255)                                         ccTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCTACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaaag                                                                                                             <  7:283960/149‑1 (MQ=255)                                         ccTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCaga                                                                                                                >  7:31401/1‑146 (MQ=255)                                             gacATACATCGGATAAGCCACAACGTCATACCCGACATTTGCCGGGGCCACCGGCACGGGCGACGGCACGACGACCTGCTGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTTGTAATAGGTATCGGCCACCAGAAAGAACa                                                                                                         <  1:276747/147‑1 (MQ=255)                                               cATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaa                                                                                                               >  7:216685/1‑141 (MQ=255)                                               cATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCagaa                                                                                                               <  8:216685/141‑1 (MQ=255)                                                aTACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAgcg                                                                                                      <  4:22788/149‑1 (MQ=255)                                                  aCATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGca                                                                                                    >  7:96083/1‑149 (MQ=255)                                                    aTCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGATGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                                                  <  8:340841/117‑1 (MQ=255)                                                    aTCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTaa                                                                                                                                  >  7:340841/1‑117 (MQ=255)                                                    aTCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGcac                                                                                                   >  4:99586/1‑148 (MQ=255)                                                     tCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCTGCCACCAGAAAGAACAGCGcacca                                                                                                 >  8:196498/1‑149 (MQ=255)                                                          tAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCaa                                                                                            >  2:351494/1‑149 (MQ=255)                                                            aGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAAc                                                                                           <  1:20131/148‑1 (MQ=255)                                                             gCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACca                                                                                         >  8:64095/1‑149 (MQ=255)                                                               caccacGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACcaca                                                                                       >  8:169382/1‑149 (MQ=255)                                                                  cacGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACtt                                                                                    >  4:140892/1‑149 (MQ=255)                                                                      tCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCcg                                                                                >  7:345399/1‑149 (MQ=255)                                                                         tACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCcgcgc                                                                             <  6:327645/149‑1 (MQ=255)                                                                          aCCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAAc                                                                                           <  6:277117/134‑1 (MQ=255)                                                                          aCCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAAc                                                                                           >  5:271015/1‑134 (MQ=255)                                                                          aCCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAAc                                                                                           >  5:277117/1‑134 (MQ=255)                                                                          aCCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAAc                                                                                           <  6:271015/134‑1 (MQ=255)                                                                            ccGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGt                                                                           >  4:6816/1‑148 (MQ=255)                                                                             cGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTaa                                                                         <  7:322033/149‑1 (MQ=255)                                                                                 aTTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAggcg                                                                     <  7:169382/149‑1 (MQ=255)                                                                                    tGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAggcggc                                                                   <  1:75089/148‑1 (MQ=255)                                                                                     gCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTgg                                                                                                                                     >  5:349558/1‑81 (MQ=255)                                                                                     gCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTgg                                                                                                                                     <  6:349558/81‑1 (MQ=255)                                                                                     gCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAggcgg                                                                    >  4:176903/1‑146 (MQ=255)                                                                                      ccGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAggcgg                                                                    >  3:220178/1‑145 (MQ=255)                                                                                          ggCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAAc                                                            >  4:368921/1‑149 (MQ=255)                                                                                           gCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGt                                                                                                                             >  3:79530/1‑83 (MQ=255)                                                                                           gCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGt                                                                                                                             <  4:79530/83‑1 (MQ=255)                                                                                            ccACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACt                                                           >  4:318824/1‑148 (MQ=255)                                                                                                  gCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCa                                                                  <  8:361216/135‑1 (MQ=255)                                                                                                  gCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCa                                                                  >  7:361216/1‑135 (MQ=255)                                                                                                       ggCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGcc                                               >  6:252522/1‑149 (MQ=255)                                                                                                       ggCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGcc                                               <  7:196498/149‑1 (MQ=255)                                                                                                        gCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCg                                              >  3:306652/1‑149 (MQ=255)                                                                                                              gCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTAcggc                                        >  7:327779/1‑149 (MQ=255)                                                                                                              gCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTAcggc                                        <  3:368921/149‑1 (MQ=255)                                                                                                               cACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTAcggcg                                       <  3:140892/149‑1 (MQ=255)                                                                                                                aCGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTAcggcgg                                      >  6:258242/1‑149 (MQ=255)                                                                                                                 cGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCggg                                     <  1:351494/149‑1 (MQ=255)                                                                                                                        cggcggGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTcacca                              >  5:18056/1‑149 (MQ=255)                                                                                                                        cggcggGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTcacca                              <  2:369961/149‑1 (MQ=255)                                                                                                                         ggcggGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTcaccac                             >  5:180547/1‑149 (MQ=255)                                                                                                                            gggaggACCACCCCATATTCCTGGCTGTCTGCCAGGTACAGGTAATAGGTATCGGCCACCAGAAAGAACCGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGAt                          <  4:47161/144‑1 (MQ=255)                                                                                                                            gggTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGAt                          <  6:180547/149‑1 (MQ=255)                                                                                                                             ggTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATg                         <  3:6816/149‑1 (MQ=255)                                                                                                                             ggTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATg                         >  3:164612/1‑149 (MQ=255)                                                                                                                                 gACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTaa                      <  7:64095/148‑1 (MQ=255)                                                                                                                                   ccaccaCATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACgg                   <  5:252522/149‑1 (MQ=255)                                                                                                                                    caccacATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGc                  >  1:282201/1‑149 (MQ=255)                                                                                                                                    caccacATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGc                  >  4:24389/1‑149 (MQ=255)                                                                                                                                       cacaTATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACg                    <  5:348816/144‑1 (MQ=255)                                                                                                                                       cacaTATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACg                    >  6:348816/1‑144 (MQ=255)                                                                                                                                         caTATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCcg             >  6:18056‑M2/1‑147 (MQ=255)                                                                                                                                             ttCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGggggccgttacatcgtggtgaccccgccgtacggcgtgcg         >  8:133770‑M1/1‑109 (MQ=255)                                                                                                                                              tCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGggggccgttacatcgtggtgaccccgccgtacggcgtgcgg        >  1:104187‑M1/1‑108 (MQ=255)                                                                                                                                               ccTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCcgtgcgg        >  8:31401‑M2/1‑141 (MQ=255)                                                                                                                                                cTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGggggccgttacatcgtggtgaccccgccgtacggcgtgcgggt      <  2:104187‑M1/149‑44 (MQ=255)                                                                                                                                                cTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGt                        <  4:322902/131‑1 (MQ=255)                                                                                                                                                cTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGt                        >  3:322902/1‑131 (MQ=255)                                                                                                                                                cTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCcgtgcggg       >  8:345399‑M2/1‑140 (MQ=255)                                                                                                                                                 tGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCCGTgcgggcg     >  8:136084/1‑149 (MQ=255)                                                                                                                                                   gCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGggggccgttacatcgtggtgaccccgccgtacggcgtgcgggtcag   <  7:133770‑M1/149‑47 (MQ=255)                                                                                                                                                   gCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCCGTgcgggcgg    <  3:318824/148‑1 (MQ=255)                                                                                                                                                    cTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCCGTGCGGGCGGTa  <  5:258242/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGATGGCACTGGTAGCGGTCCTGCTCCTGCTGTTCCTGACCCTGGCCATACATCGGATAAGCCACCACGTCATACCCGCCATTTGCCGGGGCCACCGGCACGGGCGACGGCACGCCGACCGGCGGGTTGACCACCACATATTCCTGGCTGTCTGCCAGGTACTGGTAATAGGTATCGGCCACCAGAAAGAACAGCGCACCACCCAACCACACTTCCCGCGCGTAAGGCGGCAGATAACTGACCCGCACGCCGTACGGCGGGGTCACCACGATGTAACGGCCCCCGTGCGGGCGGTA  >  NC_002947/5436849‑5437144

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