BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,260,675	C→G	5.4%	A711P (GCG→CCG)	PP_3734 ←	ABC transporter permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,260,675	0	C	G	5.4%	415.4 / 10.1	129	A711P (GCG→CCG)	PP_3734	ABC transporter permease
Reads supporting (aligned to +/- strand): ref base C (58/64); new base G (4/3); total (62/67)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 134 total aligned reads displayed.

GACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCTGTCA  >  NC_002947/4260526‑4260823

gACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCgcgc                                                                                                                                                       <  2:2298698/150‑1 (MQ=255)

 accaggcaccaggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCGGGGGGCCCAGCCGCGcc                                                                                                                                                      <  1:1694370/150‑1 (MQ=255)

     ggcacaaggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑Gc                                                                                                                                                 <  1:735410/150‑1 (MQ=255)

        accaggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCa                                                                                                                                              >  2:1076213/1‑150 (MQ=255)

          caggcCAGCAACAGCCTAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAACCGCGCCTG‑GCCCAgg                                                                                                                                            <  2:1217048/150‑1 (MQ=255)

          caggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAgg                                                                                                                                            <  2:1629550/150‑1 (MQ=255)

            ggcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGt                                                                                                                                          >  2:1427184/1‑150 (MQ=255)

             gcCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTc                                                                                                                                         >  1:202890/1‑150 (MQ=255)

             gcCAGCAACAGCCCAAGCGGAATCGCAAGCAGTACGGTGAAGAAGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCTCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTc                                                                                                                                         <  2:1593466/150‑1 (MQ=255)

               cAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTcag                                                                                                                                       >  2:2016361/1‑150 (MQ=255)

                aGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTcagc                                                                                                                                      <  1:1134848/150‑1 (MQ=255)

                aGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTcagc                                                                                                                                      <  1:1621416/150‑1 (MQ=255)

                   aaCAGCCCAAGCGGATTTGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGGGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAgg                                                                                                                                   <  2:1877755/150‑1 (MQ=255)

                   aaCAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAgg                                                                                                                                   >  1:347182/1‑150 (MQ=255)

                     cAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTTGGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGtt                                                                                                                                 <  1:2367663/150‑1 (MQ=255)

                     cAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGtt                                                                                                                                 <  1:990238/150‑1 (MQ=255)

                       gCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGATTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCAACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtga                                                                                                                               <  1:882047/150‑1 (MQ=255)

                        cccAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgat                                                                                                                              <  1:1076213/150‑1 (MQ=255)

                        cccAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgat                                                                                                                              >  2:2178922/1‑150 (MQ=255)

                         ccAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatg                                                                                                                             >  2:703130/1‑150 (MQ=255)

                          cAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCCTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatga                                                                                                                            >  2:19727/1‑150 (MQ=255)

                          cAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatga                                                                                                                            <  1:1533497/150‑1 (MQ=255)

                           aaGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGGCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtga                                                                                                                               >  2:1531896/1‑146 (MQ=255)

                           aaGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGaa                                                                                                                           <  1:218457/150‑1 (MQ=255)

                           aaGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGaa                                                                                                                           <  1:2262540/150‑1 (MQ=255)

                            aGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCA‑CCCCAGTGCCCACAGCGGTGTCAAGGGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGaaca                                                                                                                         >  1:2090228/1‑150 (MQ=255)

                                   tCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTcagca                                                                                                                                     >  1:1674503/1‑132 (MQ=255)

                                   tCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTcagca                                                                                                                                     <  2:1674503/132‑1 (MQ=255)

                                       cagcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACAcc                                                                                                               >  1:126085/1‑150 (MQ=255)

                                        agcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatga                                                                                                                            >  1:1357970/1‑136 (MQ=255)

                                        agcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatga                                                                                                                            <  2:1357970/136‑1 (MQ=255)

                                        agcagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCg                                                                                                              <  1:703130/150‑1 (MQ=255)

                                          cagTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑CCCCAGGGTCTTCAGGTTGATGAACAACGCCACACCGGc                                                                                                            <  1:63155/150‑1 (MQ=255)

                                           agTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgat                                                                                                                              >  2:2316156/1‑131 (MQ=255)

                                           agTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgat                                                                                                                              <  1:2316156/131‑1 (MQ=255)

                                             tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAg                                                                                                                                             <  2:2191520/114‑1 (MQ=255)

                                             tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAg                                                                                                                                             >  1:2191520/1‑114 (MQ=255)

                                             tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGaac                                                                                                                          >  2:2048647/1‑133 (MQ=255)

                                             tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGaac                                                                                                                          <  1:2048647/133‑1 (MQ=255)

                                             tACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCAc                                                                                                         >  1:506283/1‑150 (MQ=255)

                                                ggTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTca                                                                                                                                        >  1:1466686/1‑116 (MQ=255)

                                                ggTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGcc                                                                                                      >  1:2162726/1‑150 (MQ=255)

                                                  tGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCaa                                                                                                    <  2:506283/150‑1 (MQ=255)

                                                  tGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCaa                                                                                                    <  2:2264895/150‑1 (MQ=255)

                                                    aaGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAgg                                                                                                  <  1:938697/150‑1 (MQ=255)

                                                     aGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgat                                                                                                                              >  1:2040757/1‑121 (MQ=255)

                                                       ccGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTa                                                                                               >  1:1036220/1‑150 (MQ=255)

                                                        cGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCcaca                                                                                                                 <  2:328425/131‑1 (MQ=255)

                                                              gCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGt                                                                                        <  1:998190/150‑1 (MQ=255)

                                                                aTGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCCTGCGCCGCAAACCCACTGCACACAGCGGTCACAAC‑GGTCCAACCGCCCGTGAGTCCAGTGGCATGAGGTTTATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTg                                                                                      >  2:1463750/1‑150 (MQ=255)

                                                                aTGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCaa                                                                                                    <  2:1913100/136‑1 (MQ=255)

                                                                aTGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCaa                                                                                                    >  1:1913100/1‑136 (MQ=255)

                                                                 tgcgtgtCGGCCCCAGGCTCAACCCGCCCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCACCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgc                                                                                     <  2:1150173/147‑2 (MQ=255)

                                                                 tGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTTGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTga                                                                                     <  1:841776/150‑1 (MQ=255)

                                                                  gAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTaa                                                                                              >  2:2218677/1‑140 (MQ=255)

                                                                   aGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTgaga                                                                                   >  1:717012/1‑150 (MQ=255)

                                                                        cTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGt                                                                                        >  2:726096/1‑140 (MQ=255)

                                                                        cTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGt                                                                                        <  1:726096/140‑1 (MQ=255)

                                                                         tGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatg                                                                                                                             <  1:2202939/102‑1 (MQ=255)

                                                                         tGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTtgatg                                                                                                                             >  2:2202939/1‑102 (MQ=255)

                                                                          gCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCaa                                                                                                    >  2:2040474/1‑126 (MQ=255)

                                                                             ccAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTgg                                                                         <  1:1736829/150‑1 (MQ=255)

                                                                             ccAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGAGTTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTgg                                                                         <  1:2367499/150‑1 (MQ=255)

                                                                                 gCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTggcggc                                                                     >  1:2406232/1‑150 (MQ=255)

                                                                                     aaCCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAg                                                                 <  2:826990/150‑1 (MQ=255)

                                                                                      accagaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAgg                                                                                                                                            >  2:1449114/1‑74 (MQ=255)

                                                                                      accagaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAgg                                                                                                                                            <  1:1449114/74‑1 (MQ=255)

                                                                                      accagaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGG‑CCAGGCGCGGCTG‑GGCCa                                                                                                                                              <  1:350782/71‑1 (MQ=37)

                                                                                       ccagaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGt                                                               >  1:1938950/1‑150 (MQ=255)

                                                                                          gaccagTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAATGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCg                                                            <  1:2016361/150‑1 (MQ=255)

                                                                                                ttGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                      <  2:1157699/150‑1 (MQ=255)

                                                                                                ttGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                      >  1:1157699/1‑150 (MQ=255)

                                                                                                  gCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGa                                                    <  2:905339/150‑1 (MQ=255)

                                                                                                    aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACg                                                  >  1:1526032/1‑150 (MQ=255)

                                                                                                    aTGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACg                                                  >  1:2304566/1‑150 (MQ=255)

                                                                                                        gccgcACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGaa                                                      <  2:2121939/142‑1 (MQ=255)

                                                                                                        gccgcACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACgtcgg                                              >  1:314452/1‑150 (MQ=255)

                                                                                                                cccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAg                                      >  2:485582/1‑150 (MQ=255)

                                                                                                                cccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAg                                      >  1:928874/1‑150 (MQ=255)

                                                                                                                 ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGc                                     >  1:1806529/1‑150 (MQ=255)

                                                                                                                 ccAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGc                                     >  1:2012665/1‑150 (MQ=255)

                                                                                                                      gCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTc                                <  2:928874/150‑1 (MQ=255)

                                                                                                                        ccACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAg                              >  1:1236566/1‑150 (MQ=255)

                                                                                                                         cacaGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAgg                             >  1:2158862/1‑150 (MQ=255)

                                                                                                                              cGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGgcgcgc                        <  2:126085/150‑1 (MQ=255)

                                                                                                                                  gCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGc                    <  1:1427184/150‑1 (MQ=255)

                                                                                                                                    cAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc                  <  2:1526032/150‑1 (MQ=255)

                                                                                                                                     aaCTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCg                 >  2:1619017/1‑150 (MQ=255)

                                                                                                                                      ccTGGGCCAGGCGCGGCTG‑GGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGtt                                                                                       <  1:244377/78‑1 (MQ=37)

                                                                                                                                      ccTGGGCCAGGCGCGGCTG‑GGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGtt                                                                                       >  2:244377/2‑79 (MQ=37)

                                                                                                                                      ccTGGGCCAGGCGCGGCTG‑GGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGa                >  1:606671/2‑150 (MQ=255)

                                                                                                                                       cTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACgtcggt                                             >  1:1388702/1‑120 (MQ=255)

                                                                                                                                        tgggccAGGCGCGGCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTgg                                                                         <  2:593648/87‑1 (MQ=255)

                                                                                                                                           gccAGGCGCGGCTG‑GGCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGaccacc           >  1:798488/2‑150 (MQ=255)

                                                                                                                                            ccAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCAccc          >  2:1838977/1‑150 (MQ=255)

                                                                                                                                             cAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCg         >  1:1607659/1‑150 (MQ=255)

                                                                                                                                              aGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc                  <  2:282445/140‑1 (MQ=255)

                                                                                                                                              aGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTc                  >  1:282445/1‑140 (MQ=255)

                                                                                                                                                 cgcgCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCTg     <  1:614077/150‑1 (MQ=255)

                                                                                                                                                    gcCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCTGTCa  <  1:2148728/150‑1 (MQ=255)

GACCAGGCACCAGGCCAGCAACAGCCCAAGCGGAATCGCCAGCAGTACGGTGAAGCCGCTGAGCATGAGTGTCTGCCCCAGGCTCAACCAGACCAGTTGCATGCGCCGCACCCCCAGTGCCCACAGCGGTGCCAACTGGCCCAGCCGCGCCTG‑GCCCAGGGTCAGCAGGTTGATGAACAACGCCACACCGGCCACGCCAAGGGTAAGGCTGTTGAGAGCTGCGGTGGCGGCAAAGGTGCGGTTGAAGACGTCGGTCGACCAGCGTTTCAGGCGCGCCTGCTCGACCACCCGGCTGTCA  >  NC_002947/4260526‑4260823

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