BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,171,851	T→A	6.1%	T106S (ACG→TCG)	PP_3670 ←	DMT superfamily permease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,171,851	0	T	A	6.1%	336.3 / 9.1	111	T106S (ACG→TCG)	PP_3670	DMT superfamily permease
Reads supporting (aligned to +/- strand): ref base T (49/55); new base A (4/3); total (53/58)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 113 total aligned reads displayed.

GGTAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCA  >  NC_002947/4171702‑4171997

ggTAAACATCGCCCCGCCAGGGTGCATCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGt                                                                                                                                                    <  1:2346545/150‑1 (MQ=255)

 gTAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTg                                                                                                                                                   >  1:905822/1‑150 (MQ=255)

  tAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTgg                                                                                                                                                  <  1:1221277/150‑1 (MQ=255)

  tAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTgg                                                                                                                                                  >  2:2146443/1‑150 (MQ=255)

  tAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTgg                                                                                                                                                  <  2:905822/150‑1 (MQ=255)

          gCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGcc                                                                                                                                          >  1:1815176/1‑150 (MQ=255)

            cccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAg                                                                                                                                        >  2:350300/1‑150 (MQ=255)

             ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGt                                                                                                                                       <  1:611096/150‑1 (MQ=255)

             ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGt                                                                                                                                       <  2:1662566/150‑1 (MQ=255)

             ccGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                       <  2:110644/150‑4 (MQ=255)

                 cAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCa                                                                                                                                   <  1:765851/150‑1 (MQ=255)

                  aGGGTGCACCCCCACCAGCACGCCGATACATCCCATCAACACCGCCACTCACTGACCCACTGTCACTCGCTGTGTGAGCAGCGGCGCCGACAGCGCGGTGCCCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAg                                                                                                                                  <  2:216068/150‑1 (MQ=255)

                  aGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAg                                                                                                                                  <  1:1921459/150‑1 (MQ=255)

                      tGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTAt                                                                                                                              >  2:1936643/1‑150 (MQ=255)

                          ccaccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCa                                                                                                                          <  1:350300/150‑1 (MQ=255)

                           caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                         >  1:2165986/1‑150 (MQ=255)

                           caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                         <  1:1155001/150‑1 (MQ=255)

                           caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                         <  1:1846583/150‑1 (MQ=255)

                           caccaccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCaa                                                                                                                         >  1:2171102/1‑150 (MQ=255)

                               accaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCg                                                                                                                     <  1:2146443/150‑1 (MQ=255)

                                ccaGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTa                                                                                                                               >  1:1285151/1‑139 (MQ=255)

                                     aCGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACGGTCACTCGCTCTTTGAGCAGCGGCGCCGAGAGCGCGGTGACAAGCAAAGGGGCAAGGAACTTGACGCCCGAGTCCTCGGCCAGTGGCAGGTATTGAAAGCCGGTCGTg                                                                                                               >  2:1331354/1‑150 (MQ=255)

                                      cGCCGATAAAGCCCATCACCACCGACACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGa                                                                                                              >  1:290941/1‑150 (MQ=255)

                                       gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGaa                                                                                                             >  2:1992110/1‑150 (MQ=255)

                                       gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                       >  2:823079/1‑121 (MQ=255)

                                       gCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGAGGCCACGGccgtc                                                                                                                                       <  1:823079/124‑4 (MQ=255)

                                        ccGATAAAGCCCATCACCACCGCCAACCACTGCACGACTTGCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTCTTGCAAGCCGGTCGTGAAc                                                                                                            >  1:1145733/1‑150 (MQ=255)

                                         cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAca                                                                                                           <  2:1687714/150‑1 (MQ=255)

                                         cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAca                                                                                                           <  1:221297/150‑1 (MQ=255)

                                         cGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAAc                                                                                                            <  2:2192195/149‑1 (MQ=255)

                                           aTAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGAGGTGAAcagc                                                                                                         >  1:100930/1‑150 (MQ=255)

                                               aGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGc                                                                                                     <  2:2165986/150‑1 (MQ=255)

                                               aGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGc                                                                                                     <  2:2171102/150‑1 (MQ=255)

                                                gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGTGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                    >  1:361894/1‑150 (MQ=255)

                                                gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                    >  1:1405362/1‑150 (MQ=255)

                                                gCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCACCCGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCt                                                                                                    >  1:1240328/1‑150 (MQ=255)

                                                 cccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                   >  1:337795/1‑150 (MQ=255)

                                                 cccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACAGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTg                                                                                                   <  1:1477122/150‑1 (MQ=255)

                                                  ccATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTgg                                                                                                  >  1:2041960/1‑150 (MQ=255)

                                                   cATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCATTGCCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGt                                                                                                 >  1:122835/1‑150 (MQ=255)

                                                    aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTTACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACACCATGCTGGTg                                                                                                >  1:702276/1‑150 (MQ=255)

                                                    aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGc                                                                                                                       <  1:930880/127‑1 (MQ=255)

                                                    aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGc                                                                                                                       >  2:930880/1‑127 (MQ=255)

                                                    aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTg                                                                                                >  2:291201/1‑150 (MQ=255)

                                                    aTCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTg                                                                                                <  1:1539136/150‑1 (MQ=255)

                                                        ccaccGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTca                                                                                            >  2:756179/1‑150 (MQ=255)

                                                        ccaccGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTca                                                                                            >  2:1423647/1‑150 (MQ=255)

                                                            cGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTcagcag                                                                                        >  1:1106655/1‑150 (MQ=255)

                                                              ccacccacTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCGGCAggg                                                                                      >  2:1707655/1‑149 (MQ=255)

                                                                           cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGGTCGGCAGGCGCGCGGCCCGCaa                                                                         >  2:1880022/1‑150 (MQ=255)

                                                                           cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCaa                                                                         >  2:1789640/1‑150 (MQ=255)

                                                                           cACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCaa                                                                         >  2:1416934/1‑150 (MQ=255)

                                                                              tGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTgg                                                                                                                                                  >  1:1657834/1‑74 (MQ=255)

                                                                              tGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTgg                                                                                                                                                  <  2:1657834/74‑1 (MQ=255)

                                                                              tGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGt                                                                      <  1:850351/150‑1 (MQ=255)

                                                                                  aCTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGc                                                                  >  1:1457694/1‑150 (MQ=255)

                                                                                  aCTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGc                                                                  >  2:1982669/1‑150 (MQ=255)

                                                                                   cTCNCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCa                                                                 >  2:1498012/1‑150 (MQ=255)

                                                                                       ctctTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAc                                                             >  2:597554/1‑150 (MQ=255)

                                                                                         ctTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAccg                                                           <  2:290941/150‑1 (MQ=255)

                                                                                          tttGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGTAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGccc                                                                             >  2:1860729/1‑131 (MQ=255)

                                                                                          tttGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGTAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGccc                                                                             <  1:1860729/131‑1 (MQ=255)

                                                                                          tttGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATAccgg                                                          >  2:1685590/1‑150 (MQ=255)

                                                                                               gcagcGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACcggccggc                                                     <  2:478514/150‑1 (MQ=255)

                                                                                                   cggcgCCGACCGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCGCGCAAGGTCTGCAATACCGGCCGGCGGGt                                                 <  2:718099/150‑1 (MQ=255)

                                                                                                    ggcgCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCa                                                                                   <  2:1929721/115‑1 (MQ=255)

                                                                                                    ggcgCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCa                                                                                   >  1:1929721/1‑115 (MQ=255)

                                                                                                         cGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCaa                                           <  1:1498012/150‑1 (MQ=255)

                                                                                                             aGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACg                                       <  2:1683537/150‑1 (MQ=255)

                                                                                                              gcgcGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGt                                      <  2:1457694/150‑1 (MQ=255)

                                                                                                               cgcgGTGACCAGCACCGGTGCAAGGAAGTTGACGGCCGTGCCCTCTTCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGtt                                     <  1:959372/150‑1 (MQ=255)

                                                                                                                 cgGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGa                                   <  1:1789640/150‑1 (MQ=255)

                                                                                                                 cgGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGa                                   <  2:1614800/150‑1 (MQ=255)

                                                                                                                  ggTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCCATGTCTTCCGGCGGGTGCGCAACACGTTGAg                                  >  1:603413/1‑150 (MQ=255)

                                                                                                                   gTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACCCCAGGCTGGTGCTCAGCAGGCTCAGTGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGc                                 <  1:1992110/150‑1 (MQ=255)

                                                                                                                   gTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGc                                 <  1:609748/150‑1 (MQ=255)

                                                                                                                      aCCAGCACCGGGGCACGGAAGTTGCCGGACGTGGCGTCGGCCAGTGGCACGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCa                              <  2:122835/150‑1 (MQ=255)

                                                                                                                      aCCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCa                              <  2:2203472/150‑1 (MQ=255)

                                                                                                                        cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGt                                                 >  2:1314311/1‑129 (MQ=255)

                                                                                                                        cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGa                            >  2:773983/1‑150 (MQ=255)

                                                                                                                        cAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGa                            >  2:2379863/1‑150 (MQ=255)

                                                                                                                               gtggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGc                     <  2:361894/148‑1 (MQ=255)

                                                                                                                                gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGTGCaa                                           <  2:351453/127‑1 (MQ=255)

                                                                                                                                gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCa                    <  2:2041960/150‑1 (MQ=255)

                                                                                                                                gggCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCa                    <  2:337795/150‑1 (MQ=255)

                                                                                                                                    aaGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGaa                >  1:1409107/1‑150 (MQ=255)

                                                                                                                                    aaGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGaa                <  1:243537/150‑1 (MQ=255)

                                                                                                                                     aGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTgg                       >  2:698524/1‑142 (MQ=255)

                                                                                                                                     aGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTgg                       <  1:698524/142‑1 (MQ=255)

                                                                                                                                       gAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCATCACGTTGAGCCCAGACTTGGGCAGGAAGAt             <  2:1405362/150‑1 (MQ=255)

                                                                                                                                       gAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGAt             <  1:1685590/150‑1 (MQ=255)

                                                                                                                                           ttGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCg         >  1:1672951/1‑150 (MQ=255)

                                                                                                                                            cactggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGATCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCTCAGACTTGGGCAg                   >  1:1417266/5‑139 (MQ=255)

                                                                                                                                            cactggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGATCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCTCAGACTTGGGCAg                   <  2:1417266/135‑1 (MQ=255)

                                                                                                                                             actggCCGAGGCCACGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGc       >  1:1430130/4‑150 (MQ=255)

                                                                                                                                              aCGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGACCAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGcc      <  1:2379863/150‑1 (MQ=255)

                                                                                                                                                ggCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCAt    >  1:359405/1‑150 (MQ=255)

                                                                                                                                                 gCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACCGAAATACGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATc   <  1:1331354/150‑1 (MQ=255)

                                                                                                                                                 ccactGGCCGAGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATc   >  2:223627/5‑150 (MQ=255)

                                                                                                                                                  ccGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGACCGCAAGGTCTGCAATGCCGGCCTGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCa  <  1:291201/150‑1 (MQ=255)

GGTAAACATCGCCCCGCCAGGGTGCACCACCACCAGCACGCCGATAAAGCCCATCACCACCGCCACCCACTGCCCCACTGTCACTCGCTCTTTGAGCAGCGGCGCCGACAGCGCGGTGACCAGCACCGGGGCAAGGAAGTTGACGGCCGTGGCCTCGGCCAGTGGCAGGTATTGCAAGCCGGTCGTGAACAGCAGGCTGGTGCTCAGCAGGCTCAGGGCCCGCAAGGTCTGCAATACCGGCCGGCGGGTGCGCAACACGTTGAGCCCAGACTTGGGCAGGAAGATGCCGGCCATCA  >  NC_002947/4171702‑4171997

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