BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,881,072	C→T	6.0%	R50R (CGG→CGA)	recF ←	DNA repair protein RecF

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,881,072	0	C	T	6.0%	477.4 / 13.7	150	R50R (CGG→CGA)	recF	DNA repair protein RecF
Reads supporting (aligned to +/- strand): ref base C (65/76); new base T (7/2); total (72/78)
Fisher's exact test for biased strand distribution p-value = 8.82e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.33e-01

Only 100 of 157 total aligned reads displayed.

GTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACA  >  NC_000913/3880923‑3881220

gTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGcc                                                                                                                                                      <  2:1478616/150‑1 (MQ=255)

gTACAGTCGATGCTCACTTTGAAGTCCCCCTCTTTGTCTTTCGTAAAGCCAATCGATGTCTAGCGCTCTACGCCCTGAAATCGCCCGAGCAGAACAAACGCCTCCTGCTCATGGCGAATCATGCCACCAATCTGCAAACTGCTAAAcccc                                                                                                                                                      >  2:1864303/1‑150 (MQ=255)

 tACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGccc                                                                                                                                                     <  2:2447458/150‑1 (MQ=255)

    cGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGAc                                                                                                                                                  >  2:20960/1‑150 (MQ=255)

     gTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGAcc                                                                                                                                                 >  1:650205/1‑150 (MQ=255)

      tCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACTCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCa                                                                                                                                                <  2:1020672/150‑1 (MQ=255)

      tCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCa                                                                                                                                                <  1:235796/150‑1 (MQ=255)

        gATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATg                                                                                                                                              <  2:1990513/150‑1 (MQ=255)

        gATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATg                                                                                                                                              <  1:1258333/150‑1 (MQ=255)

         aTGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGCCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATgg                                                                                                                                             >  2:36917/1‑150 (MQ=255)

         aTGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGACAACAAACGCCTCCTGCTCATGGTGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATgg                                                                                                                                             >  2:1651961/1‑150 (MQ=255)

          tGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                            <  1:1692618/150‑1 (MQ=255)

          tGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                            >  2:317365/1‑150 (MQ=255)

          tGCGGACTTTGCTGTCGCCCTGTTTCTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGc                                                                                                                                            >  1:2076561/1‑150 (MQ=255)

           gCGGACTTTGCGGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCATCCTGATCATGGCGAATGACGCGACCAATGTGCAAACTGCGAAACGCCCGACCATGGcc                                                                                                                                           <  2:783788/150‑1 (MQ=255)

             ggACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGa                                                                                                                                         <  1:1039244/150‑1 (MQ=255)

                cTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCg                                                                                                                                      <  1:2485562/150‑1 (MQ=255)

                cTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCg                                                                                                                                      <  1:2079958/150‑1 (MQ=255)

                 tttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGt                                                                                                                                     <  1:704144/150‑1 (MQ=255)

                  ttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                    >  1:59898/1‑150 (MQ=255)

                  ttGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                    >  2:1942184/1‑150 (MQ=255)

                   tGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGtat                                                                                                                                   <  2:807695/150‑1 (MQ=255)

                     cTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCTCGGCCATGGTCGGGCGtttcg                                                                                                                                 <  1:1112538/150‑5 (MQ=255)

                     cTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCTCGGCCATGGTCGGGCGtttcg                                                                                                                                 <  1:1121035/150‑5 (MQ=255)

                        tCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATg                                                                                                                              <  1:2204978/150‑1 (MQ=255)

                          gCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGc                                                                                                                            <  1:1748522/150‑1 (MQ=255)

                               gTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGATGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCa                                                                                                                       <  1:1864303/150‑1 (MQ=255)

                                  tGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCa                                                                                                                    >  2:1526983/1‑150 (MQ=255)

                                     cTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                    >  2:2290071/1‑131 (MQ=255)

                                     cTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGta                                                                                                                                    <  1:2290071/131‑1 (MQ=255)

                                         ggTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGt                                                                                                             <  2:1021339/150‑1 (MQ=255)

                                          gTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGtt                                                                                                            >  1:50185/1‑150 (MQ=255)

                                           ttAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGttt                                                                                                           >  2:1729619/1‑150 (MQ=255)

                                              aGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGc                                                                                                        >  1:568244/1‑150 (MQ=255)

                                               gCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGcc                                                                                                       <  1:380711/150‑1 (MQ=255)

                                                 cAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCAc                                                                                                     <  2:2128975/150‑1 (MQ=255)

                                                    tCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAAAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGc                                                                                                  >  2:1590207/1‑150 (MQ=255)

                                                         gTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTg                                                                                             <  1:1802641/150‑1 (MQ=255)

                                                           ctcGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGc                                                                                           >  2:769143/1‑150 (MQ=255)

                                                           ctcGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGc                                                                                           >  2:505833/1‑150 (MQ=255)

                                                           ctcGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGc                                                                                           <  1:999015/150‑1 (MQ=255)

                                                             cgcgCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCAc                                                                                         <  1:766263/150‑1 (MQ=255)

                                                              gcgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCAcc                                                                                        <  1:2470468/150‑1 (MQ=255)

                                                              gcgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCAcc                                                                                        <  1:367952/150‑1 (MQ=255)

                                                               cgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCt                                                                                       >  2:2333861/1‑150 (MQ=255)

                                                                  tctTCGCCCTGTAATCGCCCGTGGACAACAAACGCCTCCTGCTCATGGTGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAcc                                                                                    <  1:1651961/150‑1 (MQ=255)

                                                                   ctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCa                                                                                   >  2:2060809/1‑150 (MQ=255)

                                                                         ccTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGt                                                                             <  2:324879/150‑1 (MQ=255)

                                                                           tGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTa                                                                           >  1:2409505/1‑150 (MQ=255)

                                                                           tGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTa                                                                           <  2:636418/150‑1 (MQ=255)

                                                                             tAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCGGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTaaa                                                                         <  2:2076561/150‑1 (MQ=255)

                                                                             tAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGCCCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTaaa                                                                         >  2:1165795/1‑150 (MQ=255)

                                                                             tAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTaaa                                                                         <  2:397070/150‑1 (MQ=255)

                                                                                 cGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCg                                                                     <  2:1927854/150‑1 (MQ=255)

                                                                                     cGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGgag                                                                 <  2:50185/150‑1 (MQ=255)

                                                                                        ggAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATa                                                              >  1:1956741/1‑150 (MQ=255)

                                                                                         gagaACAAACGCCTCCTGCTCATGGCGAATGACGCGCCCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATaa                                                             >  2:958851/1‑150 (MQ=255)

                                                                                          agaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAg                                                            >  1:5941/1‑150 (MQ=255)

                                                                                           gaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAgg                                                           >  1:2235218/1‑150 (MQ=255)

                                                                                           gaACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAgg                                                           >  1:781153/1‑150 (MQ=255)

                                                                                                aaCGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCgaga                                                      >  1:1894613/1‑150 (MQ=255)

                                                                                                aaCGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCgaga                                                      >  1:1846392/1‑150 (MQ=255)

                                                                                                 aCGCCTCATGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGAt                                                     <  1:2387708/150‑1 (MQ=255)

                                                                                                     ctcctGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCcgc                                                 <  1:1942184/150‑1 (MQ=255)

                                                                                                      tcctGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCcgcg                                                <  1:1925885/150‑1 (MQ=255)

                                                                                                         ggCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCGGGGTTTGCCACTGCGGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGCCGAGATCCGCGGtt                                             <  2:1852772/149‑1 (MQ=255)

                                                                                                             cATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCaa                                         >  2:2487348/1‑150 (MQ=255)

                                                                                                             cATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTCGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGACATAAGGCGAGATCCGCGCTTTCaa                                         >  2:612040/1‑150 (MQ=255)

                                                                                                                ggCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGt                                      >  1:895231/1‑150 (MQ=255)

                                                                                                                  cGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTg                                    <  2:1765502/150‑1 (MQ=255)

                                                                                                                    aaTGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCg                                  <  2:650205/150‑1 (MQ=255)

                                                                                                                         cgcgACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAAt                             >  1:1864202/1‑150 (MQ=255)

                                                                                                                              ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                        <  2:1894613/150‑1 (MQ=255)

                                                                                                                              ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                        <  1:1533237/150‑1 (MQ=255)

                                                                                                                              ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                        <  2:1812588/150‑1 (MQ=255)

                                                                                                                              ccAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCgg                        <  2:1779968/150‑1 (MQ=255)

                                                                                                                               cAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGa                       <  2:1689381/150‑1 (MQ=255)

                                                                                                                                 aTCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATc                     <  2:1956741/150‑1 (MQ=255)

                                                                                                                                   cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAACGCGAGATCCGCGGGTTCAATGTTGCGGAAATCGCGGATcaa                   <  2:2155602/150‑1 (MQ=255)

                                                                                                                                    tGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaac                  >  1:1045021/1‑150 (MQ=255)

                                                                                                                                     gCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaaca                 >  1:1428788/1‑150 (MQ=255)

                                                                                                                                       aaaCTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAg               >  1:37008/1‑150 (MQ=255)

                                                                                                                                          catctataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCgg            >  1:1998047/8‑150 (MQ=255)

                                                                                                                                          cTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCgg            >  1:1196638/1‑150 (MQ=255)

                                                                                                                                           tGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCggg           <  1:325197/150‑1 (MQ=255)

                                                                                                                                            gCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGt          >  1:242842/1‑150 (MQ=255)

                                                                                                                                             cGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTg         <  2:1050101/150‑1 (MQ=255)

                                                                                                                                              tataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGa        >  1:1413221/4‑150 (MQ=255)

                                                                                                                                              tataCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGAATTTGGCACATACAAGAAAGTTAAAGCCGGGAGATAATACGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGa        >  2:1831554/4‑150 (MQ=255)

                                                                                                                                               aaaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAg       <  2:2404083/150‑1 (MQ=255)

                                                                                                                                               aaaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAg       >  1:332582/1‑150 (MQ=255)

                                                                                                                                               aaaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAg       >  2:1038264/1‑150 (MQ=255)

                                                                                                                                                aaCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGAGGAAATCGCGGATCAACAAGCGGGTGAgg      <  2:14200/150‑1 (MQ=255)

                                                                                                                                                 aCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg     >  1:290971/1‑150 (MQ=255)

                                                                                                                                                 aCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg     >  1:948566/1‑150 (MQ=255)

                                                                                                                                                 aCGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg     >  2:84800/1‑150 (MQ=255)

                                                                                                                                                 aCGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAggg     >  1:1146024/1‑150 (MQ=255)

                                                                                                                                                  cGCTCGGCCATGGTCGGGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGa    >  2:171150/1‑150 (MQ=255)

                                                                                                                                                   gCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGAc   <  2:1262344/150‑1 (MQ=255)

                                                                                                                                                    cccGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACa  >  1:1359791/1‑150 (MQ=255)

GTACCGTCGATGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACA  >  NC_000913/3880923‑3881220

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