BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,881,085	A→G	6.1%	L46P (CTC→CCC)	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,085	0	A	G	6.1%	440.4 / 13.5	147	L46P (CTC→CCC)	recF	DNA repair protein RecF
Reads supporting (aligned to +/- strand): ref base A (64/74); new base G (7/2); total (71/76)
Fisher's exact test for biased strand distribution p-value = 8.94e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.59e-01

Only 100 of 151 total aligned reads displayed.

GGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCATTG  >  NC_000913/3880936‑3881234

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

ggACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGa                                                                                                                                                       >  2:1946943/1‑150 (MQ=255)

   cTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCg                                                                                                                                                    <  1:2485562/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)

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

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

          gTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGAt                                                                                                                                             >  1:2358644/1‑150 (MQ=255)

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

                 tGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTcc                                                                                                                                      >  2:2092904/1‑150 (MQ=255)

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

                  gTTTGTCTTTGGTTAAGCCAATCGCTGTCGCGCGCGCTTCGCCCTGTAATCGCCCGGGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCa                                                                                                                                     <  2:485885/150‑1 (MQ=255)

                    ttGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGc                                                                                                                                   >  1:1424036/1‑150 (MQ=255)

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

                       tCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACg                                                                                                                                >  2:1838694/1‑150 (MQ=255)

                           tGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGCCCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTgg                                                                                                                            >  2:1726767/1‑150 (MQ=255)

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

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

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

                                aaGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTg                                                                                                                       <  2:2044696/150‑1 (MQ=255)

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

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

                                  gCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGACAGTGGAGAACAAACGCCTCCTGCCCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGcc                                                                                                                     <  2:859117/150‑1 (MQ=255)

                                   ccAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCa                                                                                                                    >  1:1927854/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

                                                    ctctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAc                                                                                                   >  2:2163489/1‑150 (MQ=255)

                                                    ctctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGCTTTGCCACTGCCGTTGGCACCTAc                                                                                                   >  1:325705/1‑150 (MQ=255)

                                                     tctTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAcc                                                                                                  <  2:1351311/150‑1 (MQ=255)

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

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

                                                        tCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGACAGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGCTTTTGCCACTGCCGTTGGCACCTACCAGa                                                                                               >  1:759680/1‑150 (MQ=255)

                                                           cccTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAg                                                                                            <  1:2042864/150‑1 (MQ=255)

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

                                                             cTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGtt                                                                                          <  1:20960/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

                                                                                   aaCGCCTCCTGCTCATGCCGAAAGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCCACCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAACAGAAAGTTAAAGCCGGGAGATAAGTCgaga                                                                    >  2:1898052/1‑150 (MQ=255)

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

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

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

                                                                                              cTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTc                                                         >  2:2123402/1‑150 (MQ=255)

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

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

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

                                                                                                         tGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGa                                              <  2:1424036/150‑1 (MQ=255)

                                                                                                              cgACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATcg                                         <  2:568244/150‑1 (MQ=255)

                                                                                                                aCCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATcgcg                                       <  2:27620/150‑1 (MQ=255)

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

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

                                                                                                                  cAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGa                                     >  1:1972904/1‑150 (MQ=255)

                                                                                                                     tCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATca                                  <  2:1036319/150‑1 (MQ=255)

                                                                                                                     tCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATca                                  <  2:2485528/150‑1 (MQ=255)

                                                                                                                      cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaa                                 >  1:782740/1‑150 (MQ=255)

                                                                                                                      cTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaa                                 <  1:1876568/150‑1 (MQ=255)

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

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

                                                                                                                        gCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATcaaca                               <  1:1573724/150‑1 (MQ=255)

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

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

                                                                                                                            aCTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCg                           <  1:1412918/150‑1 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                      gCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGTCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGAc                 <  2:1655995/150‑1 (MQ=255)

                                                                                                                                       cccGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACa                <  2:1781875/150‑1 (MQ=255)

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

                                                                                                                                           aCCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTAc            <  1:2060809/150‑1 (MQ=255)

                                                                                                                                            ccATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACa           <  2:603485/150‑1 (MQ=255)

                                                                                                                                               tGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGtc        >  2:1385754/1‑150 (MQ=255)

                                                                                                                                                ggCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGtct       <  1:317365/150‑1 (MQ=255)

                                                                                                                                                  ccGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGACACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCa     >  2:936282/1‑150 (MQ=255)

                                                                                                                                                   cGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCAt    >  2:100219/1‑150 (MQ=255)

                                                                                                                                                     aGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCATTg  <  1:2092904/150‑1 (MQ=255)

GGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCATTG  >  NC_000913/3880936‑3881234

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