BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,553,198:1	+A	5.0%	coding (1772/1797 nt)	PP_4041 ←	glycoside hydrolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,553,198	1	.	A	5.0%	409.3 / 12.0	99	coding (1772/1797 nt)	PP_4041	glycoside hydrolase
Reads supporting (aligned to +/- strand): ref base . (51/43); new base A (3/2); total (54/45)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

TCGCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGG  >  NC_002947/4553051‑4553344

tCGCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑cc                                                                                                                                                  >  2:223351/1‑150 (MQ=255)

tCGCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACCCCGGTCCAGGCAGTCGCCCCACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑cc                                                                                                                                                  >  2:610078/1‑150 (MQ=255)

tCGCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGACCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑cc                                                                                                                                                  <  2:401949/150‑1 (MQ=255)

  gCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCgg                                                                                                                                                <  2:1867545/150‑1 (MQ=255)

    tGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGTGTTTGGCCTAGTCTACTCCTGTCATGGCCGTCGCCCTTCTCATTCGCTACCCACGACGCGCCCGG‑CGCGCGACCGGGc                                                                                                                                              <  2:325730/150‑1 (MQ=255)

       cATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCAcg                                                                                                                                           <  1:549637/150‑2 (MQ=255)

        aTTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTGGGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACtg                                                                                                                                          >  1:1470393/1‑150 (MQ=255)

        aTTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTCTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACtg                                                                                                                                          >  2:924610/1‑150 (MQ=255)

              cTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGcgcg                                                                                                                                    <  2:1367978/150‑1 (MQ=255)

                 gATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTTCTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACTGACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTg                                                                                                                                 >  1:1260940/1‑150 (MQ=255)

                     ttgtCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGa                                                                                                                             >  2:793773/1‑150 (MQ=255)

                     ttgtCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGa                                                                                                                             <  2:1826720/150‑1 (MQ=255)

                       gtCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGTCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGAtt                                                                                                                           <  2:427766/150‑1 (MQ=255)

                        tCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCGTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACGGTGCGCGGTGCAGATTc                                                                                                                          >  2:463682/1‑150 (MQ=255)

                            tCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAgcg                                                                                                                      >  2:1324579/1‑150 (MQ=255)

                              tGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGcgcg                                                                                                                    <  2:772784/150‑1 (MQ=255)

                               gCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCgg                                                                                                                   <  2:1750114/150‑1 (MQ=255)

                                  ggTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGgtgt                                                                                                                >  2:652156/1‑150 (MQ=255)

                                   ttCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACACCGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGtt                                                                                                               <  2:641968/149‑1 (MQ=255)

                                   gTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCTCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGtt                                                                                                               <  2:778503/150‑1 (MQ=255)

                                      aGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTtgat                                                                                                            >  2:134527/1‑150 (MQ=255)

                                      aGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTtgat                                                                                                            >  2:1141391/1‑150 (MQ=255)

                                      aGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTACTAAACTACACACAACGCGCCCTGTCGCGCG‑CAGGGCACTGTGCGCGGTTCAGATTCAGCGCGGTGTtgat                                                                                                            >  1:34411/1‑150 (MQ=255)

                                         ccatcgAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTgatgat                                                                                                         <  2:1847338/146‑1 (MQ=255)

                                         aCCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTgatgat                                                                                                         <  2:441539/150‑1 (MQ=255)

                                         aCCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGAATCATTTCGGTGTTgatgat                                                                                                         <  1:1005165/150‑1 (MQ=255)

                                         aCCGCGAATGCCCATGTATGCTCCGTGCAGAAAGGGTTTGGCCTAGTCTCCTCCGGAAAAGTCAGTCGACCGACTCAATCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTgatgat                                                                                                         <  1:610078/150‑1 (MQ=255)

                                           cgcgAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑cc                                                                                                                                                  <  1:74125/107‑1 (MQ=255)

                                           cgcgAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑cc                                                                                                                                                  >  2:74125/1‑107 (MQ=255)

                                           cgcgAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACCATTGTCATAGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGTTTCAGCTCGGGGTTgatgattc                                                                                                       >  2:739007/1‑148 (MQ=255)

                                                    ccATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGc                                                                                              <  1:739007/150‑1 (MQ=255)

                                                    ccATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCAAAACGCGCCTGGTCGCGCT‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGc                                                                                              >  1:698529/1‑150 (MQ=255)

                                                      aTGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAgc                                                                                                                       >  1:409218/1‑123 (MQ=255)

                                                      aTGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAgc                                                                                                                       <  2:409218/123‑1 (MQ=255)

                                                      aTGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTg                                                                                            >  1:75233/1‑150 (MQ=255)

                                                         tATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTtgat                                                                                                            <  2:1738123/131‑1 (MQ=255)

                                                         tATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTtgat                                                                                                            >  1:1738123/1‑131 (MQ=255)

                                                         tATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAg                                                                                         >  2:766674/1‑150 (MQ=255)

                                                         tATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAg                                                                                         <  1:1198049/150‑1 (MQ=255)

                                                          aTGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAgg                                                                                        <  1:652156/150‑1 (MQ=255)

                                                          aTGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAgg                                                                                        >  1:801509/1‑150 (MQ=255)

                                                            gCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGCTGATGCCTATCTGGCTGGAGGcc                                                                                      >  2:1542805/1‑150 (MQ=255)

                                                                cGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTggg                                                                                  >  1:1701719/1‑150 (MQ=255)

                                                                 gTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTgggg                                                                                 >  1:615762/1‑150 (MQ=255)

                                                                 gTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTgggg                                                                                 >  2:188835/1‑150 (MQ=255)

                                                                  tGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTggggg                                                                                >  2:1559850/1‑150 (MQ=255)

                                                                            gTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGa                                                                      >  1:254216/1‑150 (MQ=255)

                                                                                   cTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATgcg                                                               >  1:954365/1‑150 (MQ=255)

                                                                                   cTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATgcg                                                               >  1:52048/1‑150 (MQ=255)

                                                                                    tAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAg                                                                                         <  1:354207/123‑1 (MQ=255)

                                                                                    tAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAg                                                                                         >  2:354207/1‑123 (MQ=255)

                                                                                    tAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATgcgc                                                              >  2:64790/1‑150 (MQ=255)

                                                                                       tCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTgg                                                           >  1:687522/1‑150 (MQ=255)

                                                                                           cTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTggg                                                                                  <  1:1202745/123‑1 (MQ=255)

                                                                                           cTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTggg                                                                                  >  2:1202745/1‑123 (MQ=255)

                                                                                           cTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTACATGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTggccgg                                                       >  1:1452034/1‑150 (MQ=255)

                                                                                            tCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGa                                                                                                                             >  2:1315684/1‑79 (MQ=255)

                                                                                            tCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGa                                                                                                                             <  1:1315684/79‑1 (MQ=255)

                                                                                            tCCGGTCAAGGCAGTCGCCCTACTCATTAGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTggccggc                                                      >  1:816253/1‑150 (MQ=255)

                                                                                             ccGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTggccggcc                                                     <  1:1125012/150‑1 (MQ=255)

                                                                                              cGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCAGAGCATGTTCTGtcctgccg                                                    >  2:815194/1‑150 (MQ=255)

                                                                                               ggTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGccggccgg                                                   <  1:1056174/150‑1 (MQ=255)

                                                                                               ggTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGccggccgg                                                   <  1:924610/150‑1 (MQ=255)

                                                                                                 tCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTTCAGATGCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCGGCGGTTTCGCTGGCCGGCCGGGt                                                 >  1:575448/1‑150 (MQ=255)

                                                                                                   aaGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGCTGCCATCGATCTGGCGGAAGGCGTGGGGGTAGTTGCCGAGCATGCGCTGGCCGGCGGGGGCg                                               >  2:912648/1‑150 (MQ=255)

                                                                                                   aaGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGG‑CGCGCGACCGGGc                                                                                                                                              <  1:1185936/55‑1 (MQ=37)

                                                                                                   aaGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGG‑CGCGCGACCGGGc                                                                                                                                              >  2:1185936/1‑55 (MQ=37)

                                                                                                     ggCATTCGCCCTACTCATTCGCTACCCACCCCGCGCCAGGTCGCGCG‑GCGGGGACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGACTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTa                                             <  2:1260940/150‑1 (MQ=255)

                                                                                                       cAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCt                                                                                     >  2:1825416/1‑108 (MQ=255)

                                                                                                       cAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCt                                                                                     <  1:1825416/108‑1 (MQ=255)

                                                                                                       cAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCAGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACt                                           >  2:523933/1‑150 (MQ=255)

                                                                                                            gCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCg                                      >  2:343941/1‑150 (MQ=255)

                                                                                                             cccTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCgg                                     >  2:314891/1‑150 (MQ=255)

                                                                                                               cTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGc                                                                                                       >  2:1598858/1‑82 (MQ=255)

                                                                                                               cTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGc                                                                                                       <  1:1598858/82‑1 (MQ=255)

                                                                                                                 aCTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCcag                                 <  1:343941/150‑1 (MQ=255)

                                                                                                                   tCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCcagca                               <  2:687522/150‑1 (MQ=255)

                                                                                                                   tCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCcagca                               >  1:377671/1‑150 (MQ=255)

                                                                                                                   tCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCcagca                               >  2:1715370/1‑150 (MQ=255)

                                                                                                                     aTTCGCTACCCACAACGCGCCCGGTCGCTCG‑CCGGGCACTTTGCGCGGTGCATATTCAGCGCGGTGTTGATGATGCCTATCTTGCTGAAGGCCTGTGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGTTTCGTACTGTTCGGCCAGCAgg                             <  2:392933/150‑1 (MQ=255)

                                                                                                                      ttCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGc                                    <  1:1479012/142‑1 (MQ=255)

                                                                                                                      ttCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGc                                    >  2:1479012/1‑142 (MQ=255)

                                                                                                                          gatctcacaACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTGGATGATGCCTATCAGGCTGAAGGCCTGGGGCAGGGAGCTGGACATACACTGCCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAg                        <  1:975807/145‑1 (MQ=255)

                                                                                                                          cTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaa                        <  1:1559850/150‑2 (MQ=255)

                                                                                                                          cTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCAGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAg                        <  1:523933/150‑1 (MQ=255)

                                                                                                                           tACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaac                       <  2:1701719/150‑3 (MQ=255)

                                                                                                                           tACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCTGCTGGCCGaac                       <  2:1470393/150‑3 (MQ=255)

                                                                                                                             cccACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTgg                                                                                   <  1:1446499/88‑1 (MQ=255)

                                                                                                                             cccACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTgg                                                                                   >  2:1446499/1‑88 (MQ=255)

                                                                                                                              ccACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCg                    >  2:245816/1‑150 (MQ=255)

                                                                                                                               cacaACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGCTTCAACACTGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGt                   >  1:1294591/1‑150 (MQ=255)

                                                                                                                               cacaACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGt                   >  2:901268/1‑150 (MQ=255)

                                                                                                                                acaACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGtt                  <  1:64790/150‑1 (MQ=255)

                                                                                                                                acaACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGtt                  <  2:801509/150‑1 (MQ=255)

                                                                                                                                  aaCGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGc                >  2:897687/1‑150 (MQ=255)

                                                                                                                                    agtgcCCGG‑CGCGCGACCGGGCGCTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGcac              >  1:670309/4‑150 (MQ=255)

                                                                                                                                    agcgcCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGcac              >  1:1438646/2‑150 (MQ=255)

                                                                                                                                     gcgcCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGcaca             <  1:1537155/150‑1 (MQ=255)

                                                                                                                                      tgcCCGG‑CGCGCGACCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAg            >  2:672713/2‑150 (MQ=255)

                                                                                                                                       gcCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAgg           >  2:648168/1‑150 (MQ=255)

TCGCTGACATTGAACTCGATGTTGTCGATCTGCCGGTCAGTACCGCGAATGCCCATGTATGCTCCGTGCCGAAAGGGTTTGGCCTAGTCTACTCCGGTCAAGGCAGTCGCCCTACTCATTCGCTACCCACAACGCGCCCGGTCGCGCG‑CCGGGCACTGTGCGCGGTGCAGATTCAGCGCGGTGTTGATGATGCCTATCTGGCTGAAGGCCTGGGGGAAGTTGCCGAGCATGCGCTGGCCGGCCGGGTCGTACTGTTCGGCCAGCAGGCCGAGGTCGTTGCACAGGCCGGTCAGG  >  NC_002947/4553051‑4553344

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