BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,881,082	C→T	6.1%	G47D (GGC→GAC)	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,082	0	C	T	6.1%	470.1 / 11.8	147	G47D (GGC→GAC)	recF	DNA repair protein RecF
Reads supporting (aligned to +/- strand): ref base C (63/75); new base T (7/2); total (70/77)
Fisher's exact test for biased strand distribution p-value = 8.62e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.45e-01

Only 100 of 154 total aligned reads displayed.

TGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCA  >  NC_000913/3880933‑3881231

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

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

  cGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGCCGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCg                                                                                                                                                     >  1:1812588/1‑150 (MQ=255)

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

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

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

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

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

          gCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGtata                                                                                                                                             <  1:365516/150‑1 (MQ=255)

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

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

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

              tCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGGCGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATg                                                                                                                                         <  1:1704595/150‑1 (MQ=255)

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

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

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

                     gTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCa                                                                                                                                  >  2:325197/1‑150 (MQ=255)

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

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

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

                              tGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTgg                                                                                                                         <  1:2461301/150‑1 (MQ=255)

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

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

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

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

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

                                  tAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGtttt                                                                                                                     >  1:2404083/1‑150 (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)

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

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

                                                    gcgcTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCAcc                                                                                                   <  1:962744/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)

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

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

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

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

                                                                          ccGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGga                                                                             <  2:1014042/150‑1 (MQ=255)

                                                                           cGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGgag                                                                            <  2:50185/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:2235218/1‑150 (MQ=255)

                                                                                    cAAACGCATCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCga                                                                   >  1:27620/1‑150 (MQ=255)

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

                                                                                      aaCGCCTCCTGCTCATGCCGAAAGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCCACCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTAACAGAAAGTTAAAGCCGGGAGATAAGTCgaga                                                                 >  2:1898052/1‑150 (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)

                                                                                                   cATGGCGAATGACGCGACCAATCTGCAAACTCCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCTTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCaa                                                    >  2:620792/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:1779968/150‑1 (MQ=255)

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

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

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

                                                                                                                       aTCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATc                                <  2:1956741/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:1428788/1‑150 (MQ=255)

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

                                                                                                                                             gACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGAAACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGAGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTa          >  2:57607/1‑150 (MQ=255)

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

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

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

TGCGGACTTTGCTGTCGCCCTGTTTGTCTTTGGTTAAGCCAATCGCTGTCTCGCGCTCTTCGCCCTGTAATCGCCCGTGGAGAACAAACGCCTCCTGCTCATGGCGAATGACGCGACCAATCTGCAAACTGCGAAACGCCCGACCATGGCCGAGCGTATAGATGGCTTCCAGCACGCTGGTTTTGCCACTGCCGTTGGCACCTACCAGAAAGTTAAAGCCGGGAGATAAGGCGAGATCCGCGGTTTCAATGTTGCGGAAATCGCGGATCAACAAGCGGGTGAGGGACATTACAGTCTCA  >  NC_000913/3880933‑3881231

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